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Hypersonic Flight and the Re-Entry Problem 


The Twenty-First Wright Brothers Lecture 


H. JULIAN ALLEN* 
Ames Aeronautical Laboratory, N. 1CA 


INTRODUCTION 

P* \GRESS IN TRANSPORTATION has been brought 

about more by revolutionary than by evolutionary 
changes in methods of propulsion. Over many centur- 
ies of use, sailing ships were greatly improved, yet they 
gave way to the steamship less than 100 years after 
the invention of this new means of propulsion. In less 
than 20 years the horse-drawn carriage was replaced by 
the automobile, and more recently, in this country, the 
steam locomotive has submitted to the Diesel engine. 

Although these changes have often resulted in more 
economical transportation, they have not always done 
so. Speed improvement has also been an impetus. 
True, speed and economy often go hand in hand, but 
speed alone may constitute the sole reason for change, 
particularly when the vehicle has military usefulness. 
In the history of development of the airplane, speed has 
played a particularly prominent role. 

At its inception, the airplane could hardly have been 
considered useful, militarily or otherwise. With its en 
gine, which weighed 17 Ibs. per hp., the first Wright 
brothers machine could barely fly. The initial success, 
in fact, is attributable to the knowledge and ingenuity 
of the inventors. They knew the vital need for pro- 
viding their craft with a sufficiently great wing span to 
minimize the induced drag penalty imposed by the 
heavy engine. They then devised a biplane configu- 
tation which permitted the required long span to be 
fealized with a minimum of structural weight. In 
spite of the marginal performance of the first Wright 
machine, the situation was clearly correctable. Charles 
Manly had already constructed, for the ill-fated Langley 
Aerodrome, an engine no heavier than the 12-hp. 


Presented before the Institute of the Aeronautical Sciences in 
the Natural History Building Auditorium, Smithsonian Institu 
tion, Washington, D.C 

Chief, High-Speed Research Division 


Wright engine which had produced, in three 10-hour 
tests, more than 50 hp.—a most remarkable achieve- 
ment. 

To the sagacious. no doubt, while the general future 
of the airplane seemed bright, the military future seemed 
brilliant. Although the first Wright airplane had only 
reached 31 m.p.h., it was evident that high speed would 
be a particular virtue of the aircraft to come. This 
high speed would be unusually significant for the mili 
tary airplane since it would provide an invulnerability 
to possible counterweapons. Progress in improvement 
was sufficient in the next 10 years to enable the airplane 
to play a supporting role in World War I, and, in the 
following 20 years to give it a decisive role in World 
War II. In spite of the fact that over this whole period 
the performance gains had been truly amazing, it then 
appeared that further significant increase in speed was 
unlikely. One cause was traceable to the propeller. 
As the aircraft speed had increased, the tip speed of the 
propeller had approached too close to the speed of pres 
sure propagation in air. Thus compression shock losses 
occurred which seriously decreased the attainable lift 
drag ratio of the blades, and the propulsive efficiency 
dropped accordingly. The use of thinner blade sec- 
tions helped to delay the onset of the difficulty. How 
ever, this cure clearly was short-lived at best. Even 
presuming that the propeller difficulties could be circum- 
vented, the increase in airplane drag to be anticipated 
at transonic and higher speeds would require such a large 
increase in power as to make the propulsion system 
again excessively heavy. The propeller-driven airplane 
powered by the piston engine had reached an impasse. 

A drastic change was needed and it came. Before 
World War II was over, Sir Frank Whittle in England 
and engineers at the Junkers Company in Germany had 
developed turbojet engines that gave promise of pro- 
viding the required high power at much lower weight 
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than could be attained with a piston engine and pro- 
peller combination. In the years since that war the 
promised performance of the turbojet engine has been 
realized. For the turbojet, higher speeds are not only 
desirable but necessary if high propulsive efficiency is 
to be achieved. It fell to the aerodynamicist to make 
the attainment of efficient transonic and, later, super- 
sonic flight a reality. Thin wings of low aspect ratio 
provided one solution for reducing drag; swept wings, 
another; and means for promoting favorable interfer- 
ence between aircraft components, a third. The aero- 
dynamic improvements, in general, brought new and 
severe structural problems. Nevertheless, in time the 
revolution was complete. The new breed of airplane 
did not provide the increased speed without penalty to 
range, for the aerodynamic efficiency had suffered a 
permanent setback from the wave drag which had been 
incurred on entering the supersonic regime. By careful 
aerodynamic design, however, the penalty incurred 
due to wave drag could be kept within acceptable 
bounds. The effect on range of the unavoidable de- 
crease in lift-drag ratio could be minimized by the use 
of more efficient structures to provide a more favorable 
ratio of gross weight to empty weight. 

Now that the sonic speed hurdle has been passed and 
truly supersonic flight has become commonplace, the 
quest for more speed continues, and it devolves upon the 
power plant to produce the increased thrust required 
for higher speed with a minimum increase in engine 
weight. Two engine types have received consideration 
in this regard. 
progression from the turbojet or turbojet with after- 


One is the ram-jet, which is a natural 


burner. The other is the rocket motor. 

Fig. 1 shows the weights of the several power plants 
per unit thrust horsepower as a function of flight speed. 
The attractiveness of the rocket motor from this stand- 
point is obvious. It has a second advantage in having 
no maximum limit on speed as do the air-breathing 
engines. The limit of usefulness for the piston engine 
and propeller drive occurs when the speed is sufficiently 
high that the turbojet engine performs about as effi- 
ciently but with less weight. A similar limit occurs for 
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the turbojet since at a sufficiently high speed (Mag 
the 
serves only to decrease the efficiency and increase th, 


Number of order of 4) the compressor-turbjp, 


weight in comparison with the ram-jet. For the ram. 
jet the limit of usefulness as a heat engine occurs whep 
due to the heat of compression, the temperature of the 
air entering the combustor reaches the temperature fo 
chemical equilibrium through combustion so no heat 
can be added. The problems of air-breathing engines 
have been! and are now the subject of much research 
The limit speeds will certainly be increased but prob 
ably not indefinitely. A third advantage of the rocket 
motor is that since it does not require oxygen from 
the atmosphere it is not altitude limited. It performs 
best, in fact, 7m vacuo. 

In spite of these advantages the rocket motor has th 








important disadvantage that it requires all of th 
chemicals needed for its operation to be carried aboard 
the vehicle it powers, while the other engines requir 
only fuel. Fig. 2 shows the specific impulses in pounds 
thrust per pound per second of propellant for the several 


engines as a function of speed, and indicates the ver 


inferior position of chemical rocket motors in this 
regard. 
Clearly for high-speed short-range flight the rocket’s 


advantage of light engine weight far outweighs its dis 


advantage of low specific impulse. However, for high 


speed longer range flight the situation is not so obvious 


It is the purpose in this paper to discuss such questions 
as: Can rocket vehicles compete with supersonic air 
planes on an efficiency basis for long-range flight? 
What types of rocket vehicles, if any, appear attractive 
and under what circumstances? What new problems 
occur with these vehicles and do they appear surmount 


able? 


PERFORMANCE 


There are three types of long-range, high-speed ve- 
hicles which appear to be of particular interest—the 
ballistic, the glide, and the skip rocket. The typical 
flight trajectories of these rockets are shown in Fig. 5 
The ballistic trajectory is so well known that no dis 
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HYPERSONIC FLIGHT 





cussion of it is needed, but it is perhaps desirable to dis 
cuss the other two briefly. 

The glide rocket is boosted by the rocket motor to an 
iltitude and speed such that at the end of boost the 
dynamic pressure is that required for the vehicle to 
This 


iltitude, which gradually decreases as the vehicle loses 


fly without power at some given lift coefficient. 


momentum due to drag, has been termed by Sanger 
the “equilibrium altitude.’ The aerodynamic lift re 
quired for flight is the weight less the centrifugal force 
resulting from the curved flight around the earth. 
Thus the aerodynamic lift, which equals the weight for 
low speed, becomes zero as satellite speed is approached. 
The equilibrium altitude as a function of speed, there 
fore, varies as shown in Fig. 4 for a range of wing load- 
ings from 10 to 100 Ibs. per sq.ft., and for a lift coeffi 
cient of 1 10. It is at first, I think, a little surprising 
that even for speeds closely approaching satellite speed 
The 


glide path relative to the earth below it, therefore, re 


the equilibrium altitudes remain below 250,000 ft. 


mains very nearly flat at all supersonic speeds. 

rhe skip-rocket trajectory is composed of a succes 
sion of ballistic paths each connected to the next by 
a ‘skipping phase’’ during which the vehicle enters the 
itmosphere, negotiates a turn at some given lift coeffi- 
cient, and is then ejected from the atmosphere. In 
each skip the minimum altitude must be less, of course, 
than the equilibrium altitude at which the glide rocket 
would fly at the same lift coefficient and speed, since 
in increased lift is required to execute the turn. 

For all three rocket vehicles we will be concerned 
only with those trajectories which yield the maximum 
range for a given energy input.’ For the glide rocket 
and for the skip rocket, when in the atmosphere, the lift 
coefficient for maximum lift-drag ratio must be main- 
tained to achieve this end. For the glide rocket, as 
noted earlier, the flight path angle during glide is the 
small angle required to maintain equilibrium altitude 
as speed decreases during flight. For the ballistic 
rocket the least-energy trajectory requires, as a function 
of range, the flight path angles at end of boost shown 
in Fig. 5. For the skip rocket having a very low lift- 
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drag ratio the optimum angle at end of boost approaches 
that for the ballistic rocket. As the lift-drag ratio in 
creases, the optimum angle will be less than that for 
the ballistic vehicle except at the longest ranges. This 
case is indicated in Fig. 5 for a skip rocket having a lift 
drag ratio of 6. 

It is a purpose of this paper not only to compare 
optimum performance of these three rocket vehicles 
with one another, but also with that of a typical super- 
The 


efficiency of flight is perhaps best measured by the cost 


sonic airplane powered by air-breathing engines. 
of delivering a given payload a given distance—the 
higher the cost the lower the efficiency. It is clearly 
beyond the scope of this paper to actually compute this 
cost. Rather I will use as a measure of the efficiency 
the ratio of initial weight to payload. All the com 
ponents which go to make up the initial weight (fuel, 
structure, engines, etc.) do not have the same unit cost 
nor, for a nonexpendable vehicle, is all this material 
wasted. Nevertheless, this ratio should provide a 
fair estimate of efficiency, particularly if the vehicle 
is considered expendable. 

In the following, I will first discuss the ratio of initial 


weight to final (“‘empty’’) weight and then the ratio of 
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final weight to payload so that, lastly, some reasonably 
intelligent estimates, I hope, can be made of the rela- 
tive flight efficiency. 

For the aeronautical engineer the ratio of initial 
weight to final weight, and the range, are related by the 
Breguet equation. One form of this equation is 


R = (L/D)IV In (W;,/W;,) (1) 
where 


L/D = lift-drag ratio 


W, = take-off (initial) weight 
IW, = landing (final) weight 

5 = specific impulse of the fuel 
V = flight speed 


The product /1 is equal, of course, to the thermal pro- 
pulsive efficiency times the heat value of a unit weight 
of fuel. 

For the comparison of the several vehicles it would 
clearly be desirable to develop for rockets an equation 
corresponding to the Breguet equation. This was done 
in reference 4 in the following way. First, it should be 
noted that the speed at end of rocket boost can always 


be written 
V, = 1g 1n (W;,/W,) (2) 


where g = acceleration of gravity and /, = effective 
specific impulse. The effective specific impulse would 
be the actual specific impulse of the propellant if the 
thrust were extremely large compared to the weight, 
and aerodynamic drag during boost were negligible. 
These conditions are not met in practice so the effective 
impulse must always be less than the actual, but for 
efficient designs the difference is not great. 

Next we define an “effective drag’’ force in unpow- 
ered flight which, when multiplied by the range, equals 
the kinetic energy of the rocket vehicle—that is, 


D.R = (W,/29) V2 (3) 


and an “‘effective lift’’ equal to the final weight 


L. = W,; 
From Eqs. (3) and (4) it follows that 


R = (L/D), (V,7/2g) 


From Eqs. (2) and (5), then, a rocket equation simila 


in form to Breguet’s equation is obtained, 
R = (L/D) [.V.1n (W,/Wy) 
where the effective velocity is 
Ve = V,/2 7 
With Eqs. (1) and (6) we will be able to make th 
first comparisons we desire, for it is the product 
(L/D)IV 


which constitutes our standard of excellence. 
The power-plant characteristics of Fig. 2 may be re 
plotted as /,1’, (which for airplane propulsion is /| 


as a function of flight speed. These products are shown 
in Fig. 6. It is clear that this ‘propulsion product” js 


very low generally for the rocketcraft, for only at the 
highest flight speeds are they at all comparable with the 
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air-breathing engines. This, I think, is the dismal 
picture that has too often brought on the usual rejec- 
tion of rocketcraft as reasonably efficient vehicles. 
However, it is not alone this propulsion product which 
determines the efficiency since it is also necessary to 
consider the effect of lift-drag ratio on the range-weight 
relation. 

The effective lift-drag ratios for rockets have been 
determined in reference 4. In Fig. 7(a) are shown the 
effective lift-drag ratios for the ballistic rocket, and for 
the glide and skip rockets when the maximum aerody- 
namic lift-drag ratio has the low value of 2. The skip 
vehicle is best and, oddly enough, the glider is barely 
better than the ballistic rocket. In Fig. 7(b) is shown 
the effective lift-drag ratio for the ballistic rocket, 
again, and for the glide and skip rockets having the 
moderately high aerodynamic lift-drag ratio of 6. Here, 
the glide and skip vehicles are very nearly equal and, 
The 
most striking and important feature to be noted in both 
these Figures is that for all the rocketcraft the effective 
lift-drag ratio continuously increases with increasing 


of course, the ballistic vehicle is much inferior. 


One physical explanation for this increase for 
When the range is 


range. 
the ballistic rocket is the following. 
half the circumference of the earth, the speed at end of 
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boost is required to be just satellite speed and for greater 
range no increase in energy input is required. Thus 
the effective drag, by definition, continuously decreases 
with increasing range. For the glide rocket a physical 
explanation that may perhaps be more apparent is 
that, as speed is increased to obtain longer range, the 
centrifugal force increases. Thus, less of the weight 
must be supported by aerodynamic lift, so that the 
aerodynamic drag is less. 

If one now combines the results of Figs. 6 and 7, the 
ratio of initial-to-final weight can be obtained. These 
ratios are shown in Figs. S(a) and S(b) for the rockets 
and the airplane. The assumptions here are that the 
effective specific impulse for the rockets is a presently 
obtainable value (300) in Fig. S(a) and twice this value 
(600), which might be obtainable in the future, in Fig. 
S(b). It is also assumed that the /1° product for the 
air-breathing, hydrocarbon-burning engine is slightly 
more than SOO nautical miles (see Fig. 6) and the aero 
dynamic lift-drag ratio for the airplane, the glide rocket, 
and the skip rocket is 6. 

It is seen that when the range is sufficiently great, 
the improvement in effective lift-drag ratio offsets 
the disadvantageous propulsion characteristics of the 
rocket so that the rocket vehicles on the basis of weight 
ratio become competitive with the airplane. 

It is next in order to consider the component weights, 
other than payload, that go to make up the final weight. 
As regards the propulsion system weight, the advantage 
For the 
structural weight the situation is not so clear. The 


ballistic vehicle would appear to have the advantage 


here is, as we have seen, with the rocketcraft. 


since it has no wings; however, it will tend to have the 
largest tankage weight. Consideration must also be 
given to the fact that it may experience large aerody- 
namic loads on entering the atmosphere which may 
adversely affect structural weight. The skip rocket 
will be similarly affected. A factor of primary signifi- 
cance for all of the long-range rocketcraft is that they 
are subject to intense aerodynamic heating as a result 
of the high speeds attained. Even if the heat can be 
radiated away, the high surface temperature must ad- 
versely affect structural weight. If the heat cannot 
all be radiated, then the final weight must be increased 
by the required weight of coolant needed to protect 
the vehicle. It is fitting, then, to discuss in some detail 
several of the more important factors which influence 
final weight in an effort to gain insight into the relative 


ratios of final weight to payload. 


AERODYNAMIC HEATING 


Before discussing the detailed heating problems asso- 
ciated with each of the rocket types it is well to review 
the nature of the problem from a general point of view. 
First, as indicated in Fig. 9, for long-range rockets 
the speeds at rocket burnout are generally 10,000 ft. 
per sec. or greater. In the usual case, the speed at 
landing will be very small compared to burnout speed 
so that virtually all of the kinetic energy imparted to 
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3 | ae SKIP The ballistic vehicle requires the first method aa 
2 | Or i solution for the heating problem since the convectiy time 
£ 10 | J AERODYNAMIC “/ = 6 anon heat input rates will usually be several orders of mag skip 
3 —— nitude greater than the rate at which surfaces at ten ers 
: peratures near the melting point of metals can radiat, lari 
8 " —. hone + mel ies = Fortunately,” the total eens heat ape can b equi 
: Range, thousands of nautical miles kept satisfactorily low since the fraction of the total tem] 
— kinetic energy change which must be accepted as cor wal 
7 vected heat to the vehicle is approximately an 
S/2CpA | ore’ 
HEAT EQUIVALENT TO KINETIC ENERGY si laa 
lar where C, = the frictional drag coefficient per unit of | 
| the wetted area S and Cy, = the total drag coefficient | 
- 2 J based on the reference area A. | TI 
8 10} : That is to say, by making the friction drag small | mor 
as 8} HEAT ENERGY PER POUND compared to the total drag, a large fraction of the heat | craft 
Sz | developed is given to the atmosphere (wasted in shock | ysya 
as 6r waves, etc.), and the heat convected to the vehicle is | tg hj 
35 4p SPEED kept small. Thus for ballistic vehicles the solution to | tane 
~7 e the heating problem is to employ blunt shapes which | res] 
= 2b ICE TO 1000°F STEAM have high pressure-drag coefficients.* In this way the the s 
= ee er heat convected can be kept to a per cent or so of the tres 
. ’ af . ad 25 kinetic energy change so that the weight of coolant re yr 
ne ee quired may be correspondingly small. failu 
ree. 10 For the glide vehicle, this solution of the heating — mete 
problem is not possible since, for the glider to be su By 
these craft must appear as heat. The heat equivalent perior to the ballistic vehicle, it must develop a hig! rang 
per pound of weight as a function of flight speed is aerodynamic lift-to-drag ratio which is incompatibk ta 
shown in Fig. 10. Also shown is the heat required to with a small ratio of frictional to total drag. However, | yal 
convert one pound of ice to steam at 1,000°F. as a the glider gradually converts its kinetic energy into | tyra] 
“measuring stick” of what a fairly good coolant can range so that, unlike the ballistic vehicle, the time rate | prob 
handle. It is evident that if all of the kinetic energy of convective heating is not large, as has been shown sale 
appeared as heat within the vehicle and if the time rate Hence for the glider there is a possibility of radiating | ttai 
of heat addition were so great that little of this heat all or nearly all the convected heat with surface tem plan 
could be radiated away, the problem might well be in- peratures which are structurally permissible. Eggers | of th 
surmountable even for some cooling system far better has shown, as indicated in Fig. 11 taken from reference | gree 
than this ice-to-steam system. To make this point 6, that the maximum temperature for radiation equilib- | paid 
clear, consider a ballistic vehicle having a flight range rium of an average surface element for a conical glide the « 
of 1,500 nautical miles which requires (see Fig. 9) vehicle could be about 1,600°F., which is within that | tor ; 
speed at end of boost of something over 15,000 ft. sec. allowable for some presently available materials. The ibou 
If the coolant to absorb the kinetic energy were three radiation equilibrium temperature is adversely affected Si 
times as effective as the ice, all of the final weight would by increasing the wing loading and accordingly satis to tl 
be in coolant (see Fig. 10) so that no payload could factory values of loading will generally be low ones iS We 
be carried. A practical vehicle could obviously not be Moreover, while the use of coolants at such “hot pher 
built if the coolant weight required became even a spots” as the bow of the body or wing leading edges shai 
large fraction of the final weight. when present, would probably be required,” 7 the weight nent 
How can this situation be avoided? One of two solu- of the coolant should not necessarily be excessive. 
tions to the problem may be applied. If the rate of 
heat input is extremely high so that only a small part of * This idea has no doubt occurred to many. It has be wher 
the heat convected to the vehicle can be radiated away, brought to my attention, for example, that Dr. 1 ” — altit 
ae ig e es Director of the American Meteorite Museum, had suggested t 
then it is mandatory that a minimum fraction of the possible advantage of blunt shapes for missile re-entry bodies # TI 
total kinetic energy appear as heat within the vehicle. a result of examination of the shape and surface conditiot singt 


On the other hand, if an excessively large fraction of many metallic meteorites. With 
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HYPERSONIC FLIGHT 


For the skip rocket, the heating problem appears to 
be much more formidable. As with the glider, high 
lift-drag ratios must be developed if it is to be a useful 
type, sO that a low ratio of frictional to total drag can- 
not be realized. A large fraction of the kinetic energy 
change in the first skip, which is, in itself, a fairly large 
fraction of the total kinetic energy at burnout, must 
be convected to the vehicle. On the other hand, the 
time spent in the atmosphere is so small during this 
skip that the rate of convective heating is high. Eg- 
un, has indicated that a conical skip rocket 


gers,’ ag 

during the first skip would reach the very high radiation 

equilibrium temperatures shown in Fig. 11. These 
would require the extensive use of 


temperatures 
coolant. Thus the weight penalty would probably be 
so excessive as to rule out this vehicle as impractical 


or even impossible except for short-range flights. 


AERODYNAMIC LOADS 


[he effects of aerodynamic loads on rocketcraft are 
more complex than they are for more conventional air- 
craft since when loads are high aerodynamic heating is 
usually intense. Thus rocket structures when subjected 
to high stress due, say, to bending moments may simul- 
taneously be subjected to additional localized stress 
resulting from severe temperature differences within 
the structure. For hypersonic vehicles, in fact, thermal 
stress may easily be a principal cause for structural fail- 
ure as is evidenced by the explosive spalling types of 
failures which, fortunately for us, occur with many 
meteors.” 

By astute choice of materials and ingenuity in ar- 
rangement, however, the designer of a rocketcraft can 
keep these thermal stresses from reaching untenable 
values. Much work is, of course, in progress by struc- 
tural engineers to find solutions to these thermal stress 
problems’ and much more will be required before the 
same confidence in design can be provided that has been 
attained in the conventional structural design of air- 
planes. The time well may come when our knowledge 
of these thermal aspects will have advanced to the de- 
gree that only small penalty in weight will have to be 
paid for these complicating effects. On the other hand, 
the direct aerodynamic load will always be a vital fac- 
tor in determining the weight of structure, and it is 
about these loads that we will be concerned here. 

Since the aerodynamic loads are directly proportional 
to the dynamic pressure and hence the air density, it 
is Well to note that the density variation in the atmos- 
phere, over the range which is important from the load 
the expo- 


standpoint, # can be approximated by 


nential relation 


where p, is sea-level density, 8 a constant, and y the 
altitude 

his particular functional relationship leads to some 
smgularly interesting simplifications.’ The variations 


with altitude of the deceleration experienced by a ballis- 


{ND 


s| 
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tic missile in its flight down through the atmosphere 
is a case in point. The velocity can be expressed as 


3 n@ 


V = Ve 
where 
A; trajectory angle below the horizontal 
Ve = velocity on entrance to the atmosphere 
Cy = drag coefficient 
g acceleration of gravity 
1” weight 
A = reference area on which the drag coeflicient is 


based 
and the deceleration in terms of g is 
CppA VV; IW x 


If the speed at impact is less than about 61 per cent 


of entrance speed, which will usually be the case, then 


it can be shown that the maximum deceleration is 
simply 
(dl dt Z|», r= 7 Bl’; sin Oy 2 g¢ 


where ¢ is the Naperian logarithm base. That is to say 


the maximum deceleration that can be experienced is 
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EFFECT OF SIZE ON DECELERATION OF IRON SPHERES 
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MAXIMUM DECELERATION FOR MINIMUM ENERGY 
FLIGHT OF BALLISTIC VEHICLE 
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independent of the drag coefficient and, hence, of the 
vehicle’s shape. This deceleration is also independent 
of the weight. 

As an arbitrary example consider a series of solid iron 
spheres which have diameters ranging from 1,100 ft. 
to 10 ft. and for which the entrance velocity is 20,000 
ft. per sec. and flight path is vertical (@g = 90°). These 
spheres range in weight from less than 2 grains to over 
100 tons. 
12. It is seen that not only is the maximum deceler- 
ation the same for all, but the functional relation of the 
deceleration with altitude is also identical. The de- 
celeration curve is simply shifted to higher altitudes 


The decelerations for each is shown in Fig. 


the lighter the sphere. Moreover, even for the heaviest 
sphere the maximum deceleration is reached prior to 
impact, so it is to be expected that for any ballistic 
vehicle which employs a high-drag shape to minimize 
the heating problem the maximum deceleration will 
similarly be reached before impact. It is also impor- 
tant to note that for each of the iron spheres the decel- 
erations are large compared to gravity over a range of 
altitude of about 100,000 ft. 
over which the heating rates are high. 

For a ballistic vehicle which for a given range has 
the least kinetic energy at rocket burnout, both the 
velocity and trajectory are determined so that the de- 
This deceleration is given 


This is also the range 


celeration can be calculated. 
in Fig. 13. Itis seen that the worst deceleration occurs 
for a range of 4,000 nautical miles when it is something 
less than 60 times gravitational acceleration. The 
maximum deceleration falls off with increasing range in 
spite of the fact that the burnout speed is increasing 
because the angle of approach to the earth becomes 
more flattened and so the time rate of density change 
is decreasing enough to more than offset the speed in- 
crease. 

In any event, for the nose portion of a ballistic vehicle 
which is of bluff shape, necessitated by aerodynamic 
heating considerations, the decelerations shown in 
Fig. 13 are not so high as to seriously increase the struc- 
tural weight. In fact for ranges of half the earth cir- 
cumference and greater, the decelerations are within 
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human tolerance so that return of a manned satellit 
in the form of a ballistic-type vehicle seems reasonab 
even when aerodynamic lift is not employed to mip; 
mize the re-entry forces experienced. !! 

For the glide rocket, the aerodynamic load problen 
is, in the main, the same as that for the conventiona| 
airplane with the exception that we must deal, in whok 
or in part, with a hot structure. Moreover, the gys 
load problem is as yet not well defined. 

For the skip rocket, the aerodynamic loads, in coy 
trast with aerodynamic heating, do not appear to pr 
sent too severe a problem. The normal accelerations 
which occur in the first skip are shown as a function of 
total range in Fig. 14. The aerodynamic lift-drag ratj 
assumed is 6. It is seen that these normal accelerations 
generally are of the order we have become accustomed 
to in fighter airplane design. Even aside from aer 
dynamic heating consideration, however, the structural 
problem is far more serious than for the glide rocket. 


STABILITY 


The aerodynamic stability of rocketcraft presents 
some unusual problems which can influence both struc 
ture and guidance importantly. 

For the ballistic missile, Friedrich and Dore hay¢ 
given a general method for the analysis of the stability 
In reference 13, this method was used with certait 
simplifying assumptions to describe the oscillatory be 
havior of a ballistic body at supersonic speeds. The 
general solution for the angle of attack, a, is 


) 


} a Bia J ay /9) 
a=ce [CiJo(2V ke + C2 Yo(2V eo « 2 
In this equation /,( ) and J( ) are the zero-order Bessel 
functions of the first and second kinds, respectively 


The “dynamic stability’ factor is 


k, = [gooA (4610 sin Oz) | [Cp — Cr, + 

(Crg + (Cmre l a) 
while the “‘static stability” factor is 
| [Cnall o)*] 


ky = —[gpA (2871V1/ sin? 6; 





MAXIMUM ACCELERATION OF SKIP ROCKETS 
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and 
C. the drag coefficient 
c.. the rate of change of lift coefficient with 
angle of attack 
“e the rate of change of moment coefficient 
. with angular velocity 
C. the rate of change of moment coefficient 
with time rate of change of angle of 
attack 
‘eae rate of change of moment coefficient with 
angle of attack 
g ratio of the characteristic length to the 
radius of gyration 
ee & constants of integration 


The other symbols are as previously defined. 

If the body on entering the atmosphere has its axis 
misaligned by an angle ag with respect to the flight path 
but has no angular velocity, then C, is zero and the 


solution becomes 
kie ~ PY » By/2 
a/Aaptg = ¢ Jo(2v ke e : ) 


Then to illustrate the typical oscillatory behavior, 
suppose k» has a value of 10°, which is a likely one, and 
b, has in turn the values —10 (damped), 0, +10 (un- 
The angle of attack variations with altitude 
If we look first 


damped). 
would then be those shown in Fig. 15. 
it the oscillations when the dynamic stability factor 
is zero, it will be seen that during the descent as the 
air density increases the missile responds by pitching 
about The 


motion, thus, is a damped one, as seen by the envelope 


zero angle with decreasing amplitude. 
curve, but not damped in the usual sense that energy 
has been removed from the system. Rather, this be- 
havior is akin to the motion that would occur with an 
oscillating mass on a spring if the spring constant were 
to increase continually. The effect of the dynamic 
stability factor on the amplitude, in fact, only becomes 
significant in decreasing (k; = — 10) or increasing (k; = 
+10) the amplitude at altitudes below about 100,000 





OSCILLATORY MOTION OF A BALLISTIC VEHICLE 
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OSCILLATORY MOTION OF A SKIP ROCKET 
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ft. where the amplitude has already been reduced to 
one-tenth the original value. 

The effect of the that 
allowed for in the design of whatever cooling system is 
used to protect the vehicle from the aerodynamic heat 
ing that is experienced, since the local heat input rate 
will be changed by the motion from what they would 


oscillations occur must be 


be for the case when a is zero. 

The frequency can also be determined by the method 
of reference 13. It increases with decreasing altitude 
until the altitude for maximum deceleration is reached 
i.e., where the velocity is 61 per cent of the entrance 
speed—and_ then The 
typical case, may easily become as high as 10 cycles 


decreases. frequency, in a 
per sec. and hence may introduce important stresses 
due to inertial loads which must be considered in the 
structural design. Proper orientation of the vehicle by 
reaction controls prior to the entrance to the atmosphere 
can, of course, prevent the occurrence of such additional 
problems. 

For the long-range glide rocket which properly fol- 
lows the “equilibrium” trajectory, the stability prob- 
lems are generally those for conventional airplanes 
except that, at the highest speeds, when a sizable part 
of the vehicle’s weight may be offset by the centrifugal 
force due to the curved path around the earth, the dy- 
namic pressure is less than at slower speeds so the 
frequency of all oscillatory motions is correspondingly 
lower. 

The angular motion of a skip rocket during the first 
part of the skip when the vehicle is approaching the 
noted 


resembles the behavior for ballistic ve- 


During the second part of a skip when the 


earth 
hicles. 
vehicle leaves the atmosphere some interesting motions 
are possible and the behavior during a complete skip 
for a typical vehicle is thus worth some discussion. As 
an example, consider a skip rocket which consists essen- 
tially of a triangular plan-form wing with large leading- 
edge-sweep angle and with a root-chord length of about 
50 ft. For an assumed wing loading of about 20 Ibs. 
per sq. ft., a lift-drag ratio of 6, and a flight range of 
about 4,000 nautical miles, the speed on entering the 
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atmosphere will be 14,500 ft. per sec. before the first skip 
at a flight path angle, 6, of about 12° below the hori- 
The path during the first skip is shown in the 
Suppose on the approach to the 


zontal 
lower part of Fig. 16. 
earth the vehicle has no angular velocity but is pitched 
If the 


static margin is 5 per cent of root chord, then the initial 


an angle ag away from the correct trim angle. 


angular response to the increasing dynamic pressure, 
which, as with the ballistic vehicle, is mathematically 
expressible as a Bessel function of the first kind (zero 
order), is that shown on the left in the upper part of 
Fig. 16. The amplitude of the oscillation diminishes 
and frequeney increases as the bottom of the skip is 
approached. The 
slightly less than one cycle per sec. 
bility were zero the motion would diverge as the vehicle 
Now, curiously enough, while 


maximum frequency reached is 


If the dynamic sta- 


left the atmosphere. 
the motion during the outgoing flight could be repre- 
sented by the reverse of the same Bessel function, it 
might not be the same since the general solution is of 


the form 


and it is possible that C; may be zero rather than C». 
In the event that this is the case, the Bessel function 
of the second kind (zero order) gives the motion shown 
on the right in the upper part of the Figure—that is, 
the vehicle is left with an angular velocity as it leaves 
the atmosphere. This tumbling would not present a 
serious problem since the angular rate is only 1, 10 
r.p.m. which could be checked readily by reaction con- 
trols. Fortunately, it appears that the dynamic sta- 
bility would actually be so stabilizing for the case as- 
sumed as to well damp the motion as the bottom of the 
skip was approached, as is also indicated in the Figure. 

In the flight of a glide rocket, the flight path at rocket 
burnout may accidentally be higher or lower than the 
equilibrium altitude or the flight path angle may acci- 
dentally be higher or lower than that required to follow 
the equilibrium trajectory. In either event, the vehicle 
will follow a skip path before recovery to the equilibrium 
path may be effected. In this event it follows, from 
what has been said about the skip-rocket stability 
problem, that if the glide rocket oscillates, it will do so 
at lower frequency than it would at equilibrium alti- 
tude when above this altitude, and at higher frequency 
when below it. The amplitude variation would, of 
course, be the opposite of this behavior. 


COMPARISON OF FLIGHT EFFICIENCIES 


In the section on performance it was pointed out 
that the measure of efficiency for long-range flight of 
the three types of rockets and of the supersonic air- 
plane would be chosen as the ratio of the initial weight 
to payload. It is in order, now, in the light of the dis- 
cussion of aerodynamic heating, loading, and stability, 
to adjudge the relative efficiencies of these vehicles. 
To review, on the basis of the ratios of initial to final 
weight given in Fig. S(a), it was noted that for ranges 
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less than half the circumference of the earth, th super 
sonic airplane was the most attractive and the ballist) 
vehicle the least. The glide and skip rockets are no 
only intermediate in this regard but have approy; 
mately equal weight ratios. 

On the other hand, for the ballistic vehicle the pay 
load will generally be a larger fraction of the final weight 
than for the other vehicles since (1) by using a high 
drag shape for the re-entry body the coolant weight t 
protect it from aerodynamic heating will generally by 
small; (2) the re-entry body is relatively small and ro 
bust so that its structural weight should be low in spite 
of the large drag and side forces that may be exper 
enced; (3) if, then, the payload is a large fraction of the 
re-entry body weight, the propellant tankage weight 
which is directly geared to this weight, can be kept 
small through clever design; and (4+) the motor weight 
as for all the rocketcraft, will be low (see Fig. | 

For the glide rocket 1n comparison with the airplane 
the ratio of payload to final weight has a chance of be- 
ing superior by virtue of the engine weight advantag 
However, it was noted that low wing loading appears 
to be a necessity for these vehicles if little or no coolant 
is to be required, and thus this rocketcraft does not 
have the natural robustness of the ballistic re-entry 
body. Rather, the structural problem is more nearly 
that for the conventional airplane with the exception 
that the surface, at least, must be able to withstand 
high temperatures. Success or failure of glide rockets 
is clearly to be determined by the state of the metallur 
gical and structural arts and the ingenuity of the 
designer. 

The skip rocket is certainly unattractive for the 
present at least. Its coolant requirements are severe 
Moreover, the structural loads are fairly large and occur 
when aerodynamic heating is intense. Thus, it is t 
be expected that the disadvantage of high coolant and 
structural weight will far more than offset the ver 
marginal advantage it has as regards initial-to-final 
weight ratio when compared with the glide rocket. 

In short, it is my opinion that the ballistic vehicle 
can compare very favorably with the supersonic air 
plane for long range as well as short range flight and the 
glide rocket may also prove to be attractive, but not 


so the skip rocket. 


NEW PROBLEMS OF HYPERVELOCITY VEHICLES 


Up to this point the discussion of the problems oi 
rocket vehicles has been confined to the effects of phe- 
nomena which have in the past been important ones for 
lower speed aircraft and will continue to be important 
for aircraft of all speeds. Now with considerable ex- 
tension of both speed and altitude, other phenomena 
also become important. The nature of some problems 
will be altered, as a result, and new problems, of courst 
will be encountered. 

First, it is well to note that our interest in bluff bodies 
for ballistic vehicles in particular, and in rounded-nosed 
bodies generally, has changed our emphasis in aerody- 
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namics. The detached bow waves which occur with 
such bodies at high supersonic speeds complicate the 
culations of the flow-field In the 
present period, much attention is being given to such 
In addition, at the very high altitudes at- 


characteristics. 


studies. 
tained by most of the rocketcraft, the mean free path 
fair molecules can be of the same order, or long, com- 
Thus, slip- 


pared to the dimensions of the vehicles. 
dow and free-molecule-flow studies are of interest, par- 
ticularly for satellite vehicles.'? The aerodynamicist 
must deal with air having unfamiliar states and proper- 
ties. 

Second, at hypersonic speeds where, for example, air 
is greatly decelerated, it may undergo considerable 
change in composition,” the degree of change depend- 
ing upon many factors. Dissociation of oxygen and 
nitrogen molecules can occur and, in addition, thermal 
ionization of many of the constituents. It is naturally 
to be expected that the convective heat transfer will, 
asa result, be altered from what it was for the ‘‘perfect”’ 
gas, and this has been the subject of much recent re- 
search effort.!® the be- 


comes capable of radiating energy and the radiative 


decelerated gas 


Moreover, 


heat transfer must generally be considered for hyper- 


sonic vehicles, particularly for long-range ballistic 
rockets. It is not only the aerodynamic heating prob- 


lems that are affected. The fact that at very high air 
temperature the gas becomes electrically conductive’® 
introduces new problems in radio wave transmission 
In addition, a conducting gas flow can, 
The study 


and reception. 
of course, be influenced by a magnetic field. 


“magneto gas 


of such flows, which has been termed ge 
dynamics,” is still in too primitive a state to indicate 
how important a role it can play, but many interesting 
possibilities suggest themselves.” 

Third, our experience with airplanes powered by air- 
breathing engines has naturally been restricted to the 
stratosphere, or lower. Our ignorance increases with 
altitude. For rockets, literally, “‘the sky's the limit,” 
and it is not surprising that a great emphasis has now 
been placed on obtaining a more thorough understand- 
ing of the whole atmosphere.'"’ These studies are not 
ained at an understanding of the chemical and physical 
characteristics alone, but also of the occurrence of high- 
energy particles, from meteors to cosmic rays, and the 
nature of the problems they will promote. 

This discussion of new problems has only touched 
upon a few of the known phenomena which become of 
Cer- 


tainly numerous others will appear as the conquest of 


interest in consideration of high-speed rockets. 


Space progresses. Faced with such a nebulous state of 


affairs it is not surprising that our approach to these 


r 
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new problems is a cautious one. It is well to note, 
however, that in this regard the present situation is 
certainly analogous to that which the Wright brothers 
faced at the turn of the century. If we give the same 
painstaking and intelligent treatment to our problems 
as they gave to theirs a half century ago, our success 


seems assured. 
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Discussion of the Lecture 


Antonio Ferri, Director, Aerodynamics Laboratory, 
Polytechnic Institute of Brooklyn: 


to be invited to comment on this excellent 


I was very pleased 
and stim- 
ulating paper which treats the subject of long-range 
high-speed flight since I am extremely interested in 
this field. I agree completely with the speaker on his 
concluding remark that the appears 
attractive for long-range flight, and I concur in general 


rocket glider 


with many of the technical analyses presented. 

I shall limit my comments to two topics. First, I 
believe that the comparison presented between the 
heating problems encountered by using a skip trajec- 
tory and a glider trajectory is not necessarily indicative 
of the relative merits of the two solutions because the 


comparison has been made on the basis of maximum 
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temperature for radiation equilibrium. The skip ve- 
hicle can take advantage of transient effects in order to 
reduce the surface temperature by utilizing the heat 
capacity of the structure. effects if 
used to their best advantage can be important and 
In reference a skip 


vehicle These 


therefore must be considered. 
trajectory has been analyzed for a different application 
and it was shown that the transient effects due to the 
presence of the small heat sink represented by the 
structure of the vehicle can keep the surface tempera- 
ture much below the value corresponding to the maxi- 
mum temperature of radiation equilibrium. 

The second comment is related to the comparison be- 
tween ballistic vehicles, gliders, and airplanes with air 
breathing engines. I have a somewhat different opinion 
on the relative possibilities of rocket and jet engines and 
of airplanes with continuous propulsion and rocket 
glider vehicles when compared only from the point of 
view of long-range flight possibilities. 

I am much more optimistic on the possibilities of air 
breathing engines and airplanes for the application con- 
As an example, I do not find any physical 
in this analysis, we should not consider 


sidered. 
reason why, 
high energy fuels as described in the technical literature 
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(see, for example, reference 2). Then much larger im 


pulses can be obtained in ram-jets than assumed by the 
speaker (see Fig. 1). These improvements alone affect 
strongly the terms of the comparisons between air 
planes and gliders, as shown in Fig. 2. 

In connection with the high-speed long-range air- 
plane configurations, I firmly believe that for fuels con 


sidered today, a multistage airplane has much better 


possibilities than the single-stage airplane. In this 
scheme the first stage is used for take-off and accel 


eration to supersonic speed, and uses turbojets. The 
second and possibly third stages are used for the long- 
range high-speed flight. We have investigated thes 
solutions in some detail.* The two-stage solution ap 
pears very attractive for high-speed transport flight 
The first stage is the largest part of the system, and ts 
it can be used 


attached to a given airport. Therefore 


* This work has been performed with the cooperation of I 
and W. Daskin of the Gruen Applied Science Labora 


tories, Inc. 


Feldman 
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HYPFERSUNIC FLIGCReT 
for several flights scheduled at very small intervals of 
time. The second stage which travels the distance is 
much smaller than the first and therefore is less ex- 
pensive. The three-stage solution requires dropping 
the second stage; therefore in this respect it corre- 
sponds qualitatively to the solution of a rocket glider 
where a new booster is required for each flight. A com- 
parison of the possibilities of such configurations with 
the others mentioned before for a range of 10,000 miles 
is presented in Fig. 3. Here as parameter the ratio of 
initial weight of the system to the payload weight is 
A probable consistent value of this ratio 


presented. 
It can be seen that the 


fora rocket glider 1s also shown. 
airplane compares very favorably with the glider. 
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Wallace D. Hayes, Professor of Aeronautical En- 
gineering, Princeton University: I want to extend my 
congratulations to Mr. Allen for his excellent review of 
hyperballistic flight, with his clear presentation of the 
most important factors involved. Of particular inter- 
est to me is the final demise of the skip vehicle, whose 
attractiveness I have felt to have a nostalgic rather 
than a logical basis. It appears likely that space or 
security has cramped Mr. Allen's style somewhat, in 
that he has not discussed particularly novel or uncon- 
ventional ideas for solving some of the serious problems 
in the field. We should like, but of course are not 
permitted, to see what we might be using in hyperballis- 
tic flight ten years from now. Perhaps an idea must 
have been published in Life magazine before it can be 
declassified for a Wright Brothers Lecture. 

The considerations which have been presented are for 
relatively inanimate, suitably engineered payloads for 
which the large re-entry decelerations of the ballistic 
missile may be tolerated. With a human or canine 
payload the considerations would be different and would 
appear to preclude the use of a purely ballistic missile. 
The question of maneuverability also lies outside the 
scope which has been chosen by Mr. Allen for his lec- 
ture. In this respect we may note that the conventional 
concept of maneuverability is inapplicable at hyper- 
sonic speeds. Busemann noted in 1944 that with a 
fixed LD and a speed of 1/+/2 of satellite speed the 
same amount of energy must be expended to fly an air- 
plane around a pylon as to fly it around the world. 
And, of course, a true satellite outside the atmosphere is 
completely unmaneuverable by aerodynamic means. 

Mr. Allen's use of the Brequet formula is perhaps the 
most logical compromise to make in comparing the per- 
formance of completely dissimilar vehicles. The re- 
sults for the range of the ballistic missile and the glide 
vehicle have had to be forced into the Brequet formula, 
which is a natural formula only for the airplane. Some 
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such compromise was necessary in order to obtain a 
comparison. However, it is essential for us to keep in 
mind that the application to the ballistic missile is ar- 
bitrary with respect to the definitions of (L/D), and of 
/,V., and that the only valid comparison which may be 
made must be on the product of these two factors. 
Hence, neither Fig. 6 nor either of Figs. 7(a) and 7(b) 
has any meaning taken alone, and they acquire meaning 
only when taken together. 

The question of possible multiple staging of airplane 
configurations has already been raised by Professor 
Ferri. With regard to the comparison between the 
ballistic missile and the glide vehicle, the performance of 
both would presumably be equally improved by equiv- 
alent suitable staging. Thus, with such equivalent 
staging Mr. Allen’s comparison would not be changed. 
However, the design problems associated with multiple 
staging may be considerably different for the two types 
of vehicles, and it may be appreciably easier to design 
multistage ballistic missiles than multistage glide vehi- 
cles, or vice versa. Different staging for the two types 
of vehicles would alter the comparison. 

I am in complete agreement with Mr. Allen on the 
very great importance of rarefied gas flows in the design 
of hyperballistic vehicles. However, we should note 
that it may be very difficult to develop appreciable lift 
on a body which is in a flow so rarefied that it may be 
considered a free-molecule flow. In the limiting case of 
free-molecule flow at infinite Mach Number on a body 
at zero absolute temperature with unit accommodation 
coefficients, the lift given by the theory is strictly zero, 
independent of the shape. The drag, of course, is finite, 
so that the lift-drag ratio is zero. 

Again, I want to thank Mr. Allen for a most interest- 


ing and stimulating lecture. 


H. H. Kurzweg, Associate Technical Director for 
Aeroballistic Research, U.S. Naval Ordnance Labora 
I suppose you are all, as I am, confronted with 
to find a good way to keep ade- 


tory: 
the same problem 
quately informed on the rapid and continuous progress 
in the thinking about fast flying missiles and aircraft. 
Every day piles of papers and reports appear and dis 
appear on our desks and it becomes more and more 
difficult timewise to find the grains of gold in the sand. 
Our system of communication seems to become out- 
dated. We must find new ways to increase the amount 
of thoughts per unit time to be transplanted from one 
brain to the other. One step in the direction is the 
condensation of the results on latest investigations in a 
The 
lecture we have just heard is an excellent example of 
The wide scatter of opinions 


brief, logical, and easily understandable system. 


how it should be done. 
on the usefulness of various types of long-range vehicles 
has been focused to a few clear lines which makes all 
further discussions on the subject so much easier and 
more profitable. 

Years ago, before the ICBM 's and IRBM’s entered 
the newspaper stage, when we measured heat transfer 
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On Radial Deflection of a Cylinder of Finite 
Length With Various End Conditions 


L. TING* ano S. W. YUAN** 
Polytechnic Institute of Brooklyn 


SUMMARY 


By a method different from that given in references 1 and 2, 
the radial deformations of a thin-walled circular cylindrical shell 
subjected to a pair of equal and opposite concentrated radial 
loads are obtained in the present paper for the cases of cylinders 
of finite length with built-in, simply supported, and free ends. 
When the length of the cylinder is not too large compared to its 
radius, the solutions for different end conditions differ from each 


other considerably 


free agree with Rayleigh’s inextensional solution,’ while the solu- 
tions for cylinders with both ends simply supported agree with 
the 


Odqvist’s solution? and also that in reference 1. When 


length of the cylinder is large compared to its radius, the solu- 
tions corresponding to different end conditions approach that for 


the cylinder of infinite length, as expected 


tensional solution, Odqvist’s solution, and the solution in reference 


1 are no longer valid in this range 


SYMBOLS 


a = radius of cylinder 

h = wall thickness 

k = [3(1 — v?) (a/h)?] "4 

n = integer 

g = radially distributed load 

x = axial coordinate 

w = radial displacement 

Lin, Ax = real part of complex roots of char 
acteristic equation 

Bin, Bo = imaginary part of complex roots of 


characteristic equation 
Cins Cons Giny Ge = amplitudes defined under Eq. (7) 


Eh*/(12(1 — v?)] = 


D= bending rigidity of wall of shell per 
unit length 

E = Young’s modulus of elasticity 

L, fa, t4 = length of shell 

0 = of the order of 

Pp = concentrated radial load 

w = (w/h)/(P/Eh*? = nondimensional radial deflection 

Ajn = complex roots of characteristic 
equation 

v = Poisson’s ratio 

Y = circumferential coordinate, rad 

p = u,v, or w (see Fig. | 


INTRODUCTION 


ie REFERENCES | AND 2, the solutions for cylinders of 
finite length are constructed from the corresponding 
solution for cylinders of infinite length by the method 
Received May 15, 1957. 
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The solutions for cylinders with both ends 


Rayleigh’s inex- 


230 


of images. The solutions so obtained are correct on] 
if both ends of the cylinder are simply supported. F 


cylinders with different end conditions, such as thos 


with free or built-in ends, the method of images cannot 


be applied. It is necessary to use a different method 


In the present paper, the solution for cylinders of j1 


finite length is employed as the particular integral 


while the end conditions are fulfilled by finding th, 
proper solutions to the homogeneous differential equa 
tion. The additional work in finding the proper soly 
tion is the onus inherent in solving a wider class oj 
problems than the method of images affords. 

In reference 2, it was pointed out that the leading 
terms” = 2,4,..., in the Fourier series in ¢ obtained 
in reference | for the radial deformation of cylinders oj 
infinite length are inaccurate and that therefore the ex 
pressions based upon the Fliigge equation? should be 
used. If a combination of expressions from the Fliigge 
equation and from that in reference | are used as the 
particular integral, the task of finding the proper solu 
tion to the homogeneous Fliigge equation is extremel 
tedious. It is necessary to simplify the Fliigge equa 
tion. 

It was shown by a careful analysis of order of magni 
tude that the Fliigge equation can be simplified to the 
complete Donnell equation,* 


Sw + [12(1 — v?) a*h?] (O%w Ox) + 


(2/a*) (O8w/O¢*) (+1/a5) (O’w/ O¢') — (1/D)V 4g = 


It is called the complete Donnell equation in order t 
the 
and quoted in references | and 2. 


differentiate it from one derived in reference | 


Based upon Eq. (1), simplified expressions for radial 
deformation of cylinders of infinite length are obtained 
The those 
based upon the Fliigge equation.” 
the 
equation in addition to the analytic proof in reference 5 


in the present paper. results agree with 


This is also a con 


firmation of accuracy of the complete Donnell 
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Fic. 1. 
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evlinder. 











Ehw,.. (x,0)/P 


(x,0) 


won 


ite 


rect on] 
ed For 
iS thos, 
‘S Cannot 
iethod 
ers of ir 
integral 
ding the 


ial equa 
per s ly 


class ol 


leading 
/btained 
nders of 
> the ex 
ould be 





: Fliigge 
1 as the 
er solu 
tremel\ 


e equa 


magni 
| to the 


rder t 


ence ( 


; radial 





ytained 
those 
acon- | 
Jonnell 


ence 9 


With the simplified expressions as the particular SIMPLIFIED SOLUTIONS FOR THE RADIAL DEFORMATION 


ntegral, the corresponding solutions to the homogene- OF CYLINDERS OF INFINITE LENGTH 


us equauon of Eq. (1) are obtained for cylinders of ' ; 
odin ai Aig é' P When the complete Donnell equation,’ Eq. (1), is 
nite length with built-in, simply supported, and free a ; 
. : used as the governing differential equation for the radial 
displacement, the solution can be obtained by the same 


method as in reference 1, 


k -k x ki ~ cosne 
exp ( cos (kx a) + sin (ki x a)| + x 
7 a T 94 nN » 
> Asa” + 1)° exp N\X| Njn/t 
284 X \ + | 2(k n)'r 
vhere A, = —Asze = (Ar, + 1B) /m and A —Nip = | + 7B,,) 





re the four complex roots with positive imaginary part of the characteristic equation, 


(A> + 1)4 + 4(R/m)*A4A — 2n-? + 1-4 = O (3 

The last two terms in Ec. (3), —2”~* + n~4, are the extra terms as compared to the corresponding equation in refer- 

nce 1. It is obvious that they are negligible for n > 1 as compared to the term (A* + 1)', therefore, the terms in 
reference | for >> 1 are correct. This agrees with the conclusion in reference 2. 

Since the terms in reference 1 for nm = 2, 4,6..., ie., for n? < k*, are incorrect, it is desirable to evaluate the 


terms in Eq. (2) for n? < k? only. 
[he roots can be expressed as power series of (m &)*, and the coefficients of the power series are determined by 


equating the coefficients of like powers of (7 k)* in Eq. (3), 


nd \ —dy = (4.2, + B:,)/n 
k ) ', 1 /n 1 /n ) ay 4 1 /n\’? l ( n ) } 5 
n/] \ 2 \k 1} \k ris ee > ie | k “—7 

Since Rk iV 3(1 — v2) (a h)!/?] > 1, there exists a set of values of m such that m? > 1 and also (n k)? <1. This 


implies that when n? > 1 but (7 &)* is still small enough for the power series in Eqs. (4) and (5) to converge with 
suflicient rapidity, the results given by the expressions (4) and (5) match with those given in reference 1, with accu 


racy not less than that of the approximate differential equation, Eq. (1 
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After substituting the expressions for \,, Eqs. (4) and (5), into Eq. (2), the expressions for @,, can be simplified 


k —k| x| Rix ‘ ki x ; — —By, x : Ayn|X 
al exp cos + sin a oe > exp C1, cos ~ 
2a a a a n=2, 4 a a 
 « Pb Aalhe —B,,|\x As,l: A | 
Gj, sin ( ) Joos ne + > exp ( ) | Cs cos ( ) + G,, sin ( = ) cosnge (| 
a n=2,4 a a a J 


2 





where, for (n/k)? < 1, 
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for n > 1; the corresponding results in reference 1 can be used. 

For practical purposes, the infinite series can be approximated by the first few terms, 7 = 2, 4,..., m with m « 
and the constants C),, Gi,, and 1;,, Bi, can be approximated, respectively, by the leading term in their power series 
of (n/k)*, [see Eqs. (4), (5), and (7)]. Hence Eq. (6) is further simplified to 


m «KI ° 9 
; ” ws Dp k\? l ' ; —Uyn| X 
wW ,¢l/= 4 7» ¥COS (uy| X| /@) + Sin (up| x} /a); exp cosne (Ss 
a n=2,4 T n (l — n~-)?/ a 


where uw, = (1/2) n(n? — 1)'*-k. The first and the third terms in Eq. (6) are neglected, because the coefficients 


C:, and G;, in the second term of Eq. (6) are of the order of (k/”)* as compared to those of the first and the third 


terms. 
The counterpart of Eq. (S) in reference 1 can be obtained, 


: m<R fc 2 ne 
_ Xx 2k k \? ( 23 ae , —f,|x 
WwW.) ,¢)] = >. yCOS (H,|X| /a) + SIM (Ay) X| /a); exp cos Ne J 
a n=2,4 T 1 a 


> 


where ff, = (1/2)n?/k. 
The ratio of the amplitudes of the corresponding 


terms in Eqs. (9) and (8) is 1 to (1 — n~*)~**. The p 
percentage error of the terms in Eq. (9), the solution . qexXh 
- - ~ ~ . 7 r « 
of reference 1, as compared to that of Eq. (8) is 1 — =" ae J 

Pp 


(1 — n-*)9/?, 1.e., 35, 9, 4, and 2 per cent for m = 2, 4, caiaatl a 
@/n=400 











6, 8, respectively. This result is independent of the 


ratio of a/h. It agrees with the one in Table 1 in 5 x 10°} 
reference 2, which was obtained for the special values > 
of a/h = 100. It can also be concluded that the results = one’ 
from the complete Donnell equation will differ from _ 
those from the Fliigge equation by less than one per s 
cent. © ane 
Fig. 2 shows the radial deflection for a/h = 100, 
400, and 1,000. 0°} 
RADIAL DEFORMATION FOR CYLINDERS OF FINITE °o 2 4 6 8 7. * 3 8 


LENGTH 
Fic. 4. Maximum radial deflection of cylinders with various 
The solution, w_, developed in the previous section length-radius ratio (simply supported and built-in end 


for an infinite cylinder can be used as a particular in- 
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where @, 1S the solution of the homogeneous differential 
equation of Eq. (1) 
and w, cancel each other; therefore, Eq. 


In many problems, the leading 


terms 10 
§) should be used for w_ in the evaluation of w(x, ¢ 


I 


However. @,(¥, ¢) is determined by the values of 


y. ¢) and its derivative at the two ends x = L; and 
—L., where w., can be approximated by Eq. (S), 
because in all the other terms the exponential factor 


iminishes very rapidly as x a increases. (It is 


issumed that L; and L. are at least of the order of mag 





kx \ " - k(L, — z P k(Ly 
_¢ = pd F, cos —_ F sin 
a 2,4 ! 1 a 
[ k(L: + 
I COS 
a 
\ “— «6Of n(Li — 2 
Und un(L 
in ‘ > 41 Mh, cos + Hz, sin 
a 2,4 t a 
[ (Le + x 
' IT COs 
| a 
represents the solution for the edge effect... The 
influence of the first term on the right-hand side of Eq. 
11) is confined to the vicinity of the end x = Ly, while 
that of the second term to the vicinity of x = —L 


Since u, S are not large numbers, the effect of wy; 1s not 
localized. The radial deflection of the finite cylinder, 
v= Ww, + wv; + wi, can be represented by w’ = w, + 
1) for the part of the cylinder not too close to the ends. 
rhe four sets of constants //,,, [72,, FZ3,, and J7;, in wy 
can be determined uniquely by two boundary condi 
tions at each end. The question now is which two 
of the four boundary conditions should be prescribed 
on w@y;. This can be settled by a careful examination 
of the relative orders of magnitude of the different terms 
in the boundary conditions, while w; and w;; obey the 


following relationship: 


OW} Ox RO( wy a = 

h n) Ol(1 @ (Ow O¢) | Wy aog (13 
Ow; Ox = O(n*wy ka) = 

nk) O[(1 a) (Ow, O¢)| < Ow, AOe (14) 


The two boundary conditions for w’(x, ¢) are summa 
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nitude of the radius.) From Eq. (S), it 1s clear that 


w(x, ¢) can be represented by a finite trigonometric 
Series, 
W(X, ¥ > w,(x) cos ng (10 


Since n Sm k, the complementary integral w.(x, ¢ 


can be split into two parts w; and wy, (see reference 7 


1.e 
w(x, ¢) = wil(kx a), ¢|] + wul(u,X a), ¢ 
where 
— 2 —k(L, =. 
exp COS NL F 
a 
x . €@G. + x —k(Lo +4 | 
+ F;, sin exp cos ne 11 
a ad | 
»(Ly — Xx = (Ly 2 
exp COS Ne 7 
a a 
—~ Lo+ x —pu, (1 a | 
+ I7;, sin exp cosnge (12 
a J 
wy = —w and Owy; Ox Ow Ox ) 
(2) Simply supported edges, 
Wil = —-W and O°’; ON — (O-w O-x 16 
(3) Free edges, 
O°w)1; Ox" = —(O-w_. ON and O°*wy; OA 
O*w_ On 17 


The application of w(x, ¢) to solve unsymmetrical 


problems in the shells can be found in reference 9 
where boundary conditions for built-in edges and free 
edges equivalent to those in this paper are given. 

If the concentrated loads are applied at the center 
of the evlinder i.e., Ll, = L,-= ZL 


subjected to the same set of boundary conditions, the 


and both ends are 


radial deformation of the cylinder at the point of appli 
cation of the concentrated load can be expressed as: 


w’(0, O) = w.(0, 0) + 2PR*(Ehr x 


(1S 


rized as follows. (1) Built-in edges, where w_(0, 0) = the corresponding radial displace 
ment of the cylinder of infinite length 
7 —2 exp (—2n,L a) [2 + sin (2u,L a) — cos(2u,L a) — exp(—2u,L a 
Jy = ; : Na 
1 + 2 exp(—2u,L a) sin (2u,L a) — exp (—4u,L a 
lor built-in ends, 
—2 exp | —2u,L a) |cos (2u,L a) + sin (2u,L a) + exp | —Zy L/a 
J, = ISb 
l + 2 exp | —2u,L a) cos (2u,L a) + exp (—4u L/a 
for simply supported edges, 
J 2 exp (—2u,L a) [2 + cos (2u,L a) — sin (2u,L a) + exp (—2u,L a ) 
Je => = (1Ne 
1 + 2exp (—2y,L a) sin (2u,L a) + exp (—4u,l a 


lor free edges. 
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It is clear that /, tends to zero when L a tends to 
infinity. This means that, regardless of the kind of 
boundary conditions, the solution for the cylinders of 
finite length approaches the solution for cylinders of 
infinite length when ZL a approaches infinity. 

The radial displacement w’(O, 0) is plotted against 
L a for a eylinder with a kh = 400 in Fig. 3 for free 
ends and in Fig. + for simply supported and fixed ends. 

The Rayleigh solution’ which fulfills the boundary 


conditions for free edges, Eq. (17), gives 


wr(0, 0) = 0.149Pa?/(4DL 
or @p(O, 0) = w,(0, O)Eh p = 0.149k4*(a/L 


When the value of a / is fixed, ®p(0, 0) is in inverse 
proportion to L a and tends to zero as L a tends to 
infinity. It is clear that the Rayleigh solution cannot 
be applied when La > 1. When La = O(1), the 
results obtained in this paper for cylinders with free 


ends agree with the Rayleigh results (see Fig. 3). 


DISCUSSION 

Although in this paper only a special unsymmetrical 
problem has been treated, a general remark can be 
made regarding the accuracy of many approximate 
equations. 

The complete Donnell equation, Eq. (1), 1s accurate 
enough for unsymmetrical problems, for engineering 
purposes. 

If the cylindrical shell deforms into many circumfer- 
ential waves, i.e., 

Op O¢ > p (19) 


where p stands for uw, v, or w, Eq. (1) reduces to the 
equation in references | and 6, 
8 + [12(1 — v*) ‘a*h?] (04w/Ox*t) = (1 D)Vig 
The latter equation was also obtained by Green and 
Zerna."” 
If itis known that 


Op ady > Op Ov (20) 


Eq. (1) degenerates to (the homogeneous part) 


[(02/Oe2) + 1]? (O4p ‘a4d¢*) + 
12(1 — v*) (a7/h?) (O4p/Ox*) = 0 (21) 


Eq. (21) is equivalent to Finsterwalder’s equation for 
AM, (see reference 3), because condition (20) is necessary 
and sufficient to justify Finsterwalder’s assumption 
that 0.\/, Ox and O.\/,, Ox are negligible in the equa- 
tions of equilibrium. Both Eq. (21) and Finster- 
walder’s equation are fourth order in x instead of 
eighth order. The “solutions due to edge effect,”’ 
wy, are eliminated. 

If both conditions (19) and (20) are valid, Eq. (21) is 
further simplified to Schorer’s equation,’! 


(O8p /O¢’) + [12(1 — v*)a’/h?] (04p/Ox*) = 0 
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For ‘short’ cylinders with L a < VY 2rk, and wit 
both ends fixed or simply supported, the leading terms j; 
Wi1(v, ¢) cancel with the corresponding terms in ; 


(x, ¢) and the radial displacement can be written as 


W \X, ¥ = WAX, 9) T Wik, ¢ 
5 “ ¥) COS Ne : 





where p is an integer much greater than unity. | 
therefore solutions in refer 


(22) implies condition (19), 
ence 1 are accurate for “‘short’’ cylinders with simp] 
supported ends. 

Since wy; in Eq. (12) and the predominant terms j; 
w.. also fulfill condition (20), Odqvist’s solution,‘ which 
was based upon Schorer’s equation, is also accurate for 
“short” cylinders with L a < +7 27k |see Fig. 4 

For a thin cylindrical shell with flat heads, the en 
conditions prescribed on the shell due to the existence 
of the flat heads are equivalent to those of the simp); 
supported ends, because the flat plate is much stiffer 
than the shell in resisting radial and circumferential | 
displacements, 1.e., w = O, v = 0. and is much weaker 
than the shell in resisting lateral deflection and rotation 
ec., NV. = GO and M. = 0. 

Finally, it is worthwhile to point out that, for the 
case, 

Op O¢ = O(p) = O(adp Ov 
the differential equation (1) degenerates to the mem 
brane equation.* 


REFERENCES 


1 Yuan, S. W., Thin Cylindrical Shells Subjected to Con 
trated Loads, Quarterly of Applied Mathematics, Vol. IV, No. 1 
pp. 138-26, April, 1946 

2 Yuan, S. W., and Ting, L., On Radial Deflections of a Cylinder 
Subjected to Equal and Opposite Concentrated Radial Load 
Infinitely Long Cylinder and Finite-Length Cylinder with Simpl 
Supported Ends, Journal of Applied Mechanies, Vol. 24, No. 2, py 
278-282, June, 1957 

$ Timoshenko, S. P., Theory of Plates and Shells, (1st Ed | 
pp. 427-440, 446-449, McGraw-Hill Book Company, New York 
1940 

* Odqvist, F. K. G., Action of Forces and Moments Symmet | 
rically Distributed Along a Generatrix of a Thin Cylindrical Sh 
Journal of Applied Mechanics, Trans. ASME, Vol. 68, 
106-108, 1946 

5 Donnell, L. H., A Discussion of Thin Shell Theo y, Proceed- 
ings of the Fifth International Congress of Applied Mechanics 
pp. 66-70, Cambridge, Mass., 1938 

6 Donnell, L. H., Stability of Thin-Walled Tubes under 1 
sion, NACA Report No. 479, 1933 

7 Fliigge, W., Statik und Dynamik der Schalen, p. 118, pp. 12¢ 
127, Julius Springer, Berlin, Germany, 1934 

§ Love, A. E. H., The Mathematical Theory of Elasti 
Ed., pp. 564-583, Dover Publications, New York, 1944 

> Gruber, E., Die Berechnung Zylindrischer, biegungsstetter 
Schalen unter beliebigem Lastangriff, Internat. Vereinig. Briicket 
u. Hochbau Abh. 2, p. 196, 1934 

Green, A. E., and Zerna, W., Theoretical Ela 
p. 422, Oxford University Press, Oxford, England, 1954 

11 Schorer, H., Line Load Action on Thin Cylindrical Sh 

Trans. ASCE, Vol. 61, pp. 281-316, 1935. 


sticity, Ist Ed 








d wi 
terms 
Sing 
ae FE. L. RESLER, JR..,* 
| . 
SUMMARY 
tv. E 
in refer. | The equations describing the flow of an electrically conducting 
ie iuid in the presence of electric and magnetic fields are written 
1 simph wn with the aid of certain simplifications appropriate t 
iautical applications. In order to estimate the probable 
terms 17 significance of magneto-aerodynamic effects, some data on 
whicl nductivity of pure and ‘‘seeded”’ air are first examined. Di 
irate fi ; nsionless quantities representing the ratios of forces and of 
: : rrents are then formed and their values studied for conditions 
‘ flight in the atmosphere 
the end Some examples of magneto-hydrodynamic and magneto 
X1Stence gasdyvnamic effects in simple flows are given These include 
simply two cases Of Poiseuille flow of conducting liquids with applied 
1 stiffer ignetic fields and the case of quasi one-dimensional gas flow 
ie aa pplied electrical and magnetic fields. In the last case, 
Crential | ttractive possibilities are found for controlled acceleration or 
Weaker | deceleration of gas at subsonic and supersonic speeds, even in 
tation nstant-area channels rhe behavior of the flow is characteris 
different in different regimes of Mach Number and flow 
for the elative tc certain “significant speeds’ that are dependent 
ratio of electrical to magnetic field strengths These 
tudied, and a chart is constructed to relate the length to 
speed ratio of a maximum-acceleration constant-area channel 
It is concluded that the advantages that may accrue from 
mem igneto-aerodynamic methods are sufficiently attractive to 
ustifv the considerable research and engineering development 
will be required. Among the unsolved engineering problems 
the reduction of surface resistance of electrodes in contact 
conducting gas, development of techniques for seeding, 
Come 1 provision of the required magnetic fields in flight 
» Ne 
SYMBOLS 
‘ylinc 
oad. half-width of channel in cases of Figs. 4 and 5 
Sim pi = V4p/p in compressible flow 
). 2, py 1 = cross-section area of channel 
| = speed of light 
t Ed | = specific heat for constant pressure 
York = specific heat for constant volume 
D channel cross-section dimension [Eq. (15 
ymin | E = electric field vector (magnitude E 
a F body force 
S, | ) = current density vector (magnitude 
H magnetic field vector (magnitude // 
oceed- Hi = Hartmann Number 
nics h = perturbation magnetic field vector (magnitude / 
nh = stagnation enthalpy 
I f breakdown-potential constant [Eq. (50 
= typical body length 
12t = Mach Number 
Presented at the Aerodynamics Session, National Summer 
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mass flow per unit are 


= constant pressure gradient in Eq. (28 


= static pressure 


q = fluid velocity vector (magnitud 
VU = heat added by dissipation and/or other effects 
excluding Joule heating 
k = resistance per unit length of channel (Fig. 5 
Re = Reynolds Number 
R = magnetic Reynolds Number 
t = time 
1 = temperature 
u = flow speed (average) in channel 
u, 4; = significant speeds in channel flow 
J = typical flow speed; speed of flight 
U = electrical potential 
y, = Cartesian coordinates 
a = fractional ionization 
; - 
é = lin oli*dx 
J x 

re = viscosity coeliicient 
f = mass density 

= conductivity 

Subscripts 

x, = components of vectors, along x, ¥, iXes 
0) constant, especially initial conditi in channel 


flow 


INTRODUCTION 


I \ FLUID is a conductor of electricity the possibility 
in it that will affect the fluid flow pattern in a signifi 
to 
the 
centrifugal, and 


arises that an electric body force may be produced 


cant way. To be sure, we are not unaccustomed 


dealing with fluid flows with body forces 


€.g., 
gravitational force and the inertial 


Coriolis forces in moving coordinate systems —but 


all of these are examples where the body force is more 
The attractive thing 


or less outside of our control. 


about the electric body force (which we will call the 
e.b.f.) is that it can be controlled, insofar as the current 
and the magnetic field can be controlled, and perhaps 
made to serve useful purposes such as acceleration ot 
deceleration of flow, prevention of separation, and 
the like. 

A body force in fluid mechanics is any force that 1s 
applied by an external agency directly to the fluid 
elements, as distinguished from the pressure, for 
example, which is a force between adjacent particles. 
The e.b.f. involved here is the counterpart of the force 
that drives an electric motor; it is given by the vector 
product of the current and the magnetic field strength 
and is therefore perpendicular to both of these vectors. 
In the case of a fluid, we expect electric-current flow 
are 


through a volume rather than along a wire and 
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therefore concerned with the current-density vector 
j, measured in amperes per unit area, for example. 
Similarly, the strength of the magnetic field H is 
lines-of-force per centimeter, 
Then 


measured in square 


usually called ‘‘gauss.”’ the body force per 


unit volume is given by the vector equation 
F=jxXH (1) 


It is assumed here that the fluid is nonmagnetic —i.e., 
that it has the same permeability as empty space; 
this permeability has been taken equal to 1, since we 
propose to formulate our problem in electromagnetic 
units. 

This force can be used to pump a conducting fluid, 
such as mercury, through pipes.' Recently it has 
been used in gasdynamics laboratories to help drive 
shock tubes;” here the conducting fluid is air or another 
gas heated to high temperature by discharge of a 
condenser through it. Astrophysicists also have a 
current interest in magneto-hydrodynamics because 
there are at least two different situations in their field 
where this body force apparently acts: (1) in the sun 
and stars where gases are highly ionized by the extreme 


temperatures, and (2) in nearly empty space where 
the temperatures of gas clouds are not high but the 
ionization is nevertheless sufficient, because of the low 
densities, to make them good conductors. 

The possibility of these effects becoming important 
in aeronautics has naturally arisen, since the tempera- 
tures that occur in flight at high speeds are sufficient 
to ionize air appreciably. If, as predicted, these 
temperatures approximate the temperatures of electric 
ares, it is clear that the air becomes a conductor and 
that the e.b.f. can be important aerodynamically. 

The purpose of the present paper is to make an 
assessment of these aeronautical possibilities. How 
large are the effects likely to be, and under what 
conditions? What are some of the engineering prob- 
lems posed ? 

Before proceeding to this task, however, we should 
notice another basic phenomenon of magneto-aero- 
dynamics, equal in importance to the e.b.f. and even, 
in some cases, overshadowing it—namely, the Joule 
heating. This is one of the most familiar phenomena 
of electricity. We are used to it in metallic con- 
even more significant in 


ductors, and it may be 


magneto-aerodynamics, since the conductivities of 


gases are so much less. The electricians’ familiar 
formula /*R becomes, in this application, 
heat released per unit volume per unit time = 7? o (2) 


where o denotes the conductivity or the reciprocal of 
the specific resistance. The unit of conductance is 
the mho; here o is a material property, the conductance 
of a unit cube. The conductance of such a sample 
is proportional to its cross-sectional area and inversely 
to its length. Thus o has the units mho per unit 


length. 
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THE EQUATIONS OF MOTION 


To summarize what has been written above, th 
presence of electrical currents in a fluid introduces 
essentially two new effects, the e.b.f., j & H: ang 
Joule heating, 7? ¢. The only change required in th 
equations of fluid motion is the introduction of thes 
new terms in the equations of momentum and energy 
The equation of continuity remains unaltered. \; 


have then 
Continuity: 
(Op Ot) + div(pq) i) 


Momentum: 





(Dq Dt) + (1 p)Vp = (uw p){V°q + (1 3)VV-q! 4 


(] p)j x H 
Energy: 


c,(DT Dt) + p(Dp Dt) (77, pc) +Q 5) 


Here p, q, p, wu, 7, and o are the density, velocity 
pressure, viscosity, temperature, and electrical con- 
ductivity of the fluid; c, is its specific heat at constant 
volume; and ¢ is the time. 
all heat added to the fluid except Joule heat, thus () 


The symbol Q denotes 





includes heat conduction, viscous dissipation, ete 


One simplification has been made, for brevity, in 
namely, the viscous term has been 


but this 


writing Eq. (4) 
written for the special case of constant u 
will not affect what is done below. It should be 
remembered however that yu is actually dependent 
upon 7. The symbols j and H represent the current 
density and magnetic field strength defined above 
The derivative D Dt is the convective derivative 
i.e., the rate of change for a fluid particle. 

Eq. (5) is an expression of the fact that as heat is | 
added to, or subtracted from, the fluid medium it 
appears either in internal energy or in mechanical 
work done. This is the simple expression of the 
First Law of Thermodynamics that applies when the 
thermodynamic and mechanical forms of energy are 
the only ones concerned. Eq. (5) does not attempt 
to include electrical forms of energy in the statement 
of the First Law. It is therefore useful only if there | 
is no interchange of energy except what has _ been 
described above—1.e., Joule heating. In more com 
plicated cases, for example, where the magnetic 
permeability of the fluid is a function of temperature or 
of the field strength, a more sophisticated formulation 
of the energy equation from the first law of thermo 
dynamics is required. 

Assuming that this is not the case in most aero- 
nautical applications, the simple statement Eq. (4 
is all that is needed. It takes on an interesting form 
if it is added to Eq. (4) after scalar multiplication of 
the latter with the velocity vector (thus changing 
Eq. (4) to a statement involving work done by forces 
This form is, assuming the perfect-gas equation of state 
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(OL Ot) + (1 2)(Og* Or) + p(Oe~! Or) +q-Vhi 
(1/p)q-j X H+(/? po) + wjiV'qt+(l 3)VV-q}- 
i QO (6) 


where /i) denotes the stagnation enthalpy, c,7° + gq? 2. 
This form is useful for steady-flow problems, since the 
first three left-hand terms are then absent, and the 
left-hand side of Eq. (6) becomes the rate of increase 
of the stagnation enthalpy of a fluid particle. This 
rate involves two magneto-aerodynamic terms, one 
being the mechanical work done by the e.b.f., the other 
being the Joule heating. As we shall see below, it is 
possible for these two to oppose one another. 


Ohm's Law 


In Eqs. (5) to (5) we have introduced two new 
unknowns into the already complicated equations of 
gasdynamics —viz., j and H (assuming o to be a known 
property of the fluid). 
equations concerning these. One of these is Ohm's 


We must provide additional 


Law: 
j o}\E +8 x H} (7) 


This should be thought of as an empirical result, de- 
scribing in a macroscopic way how the electric current 
flows through a volume of conducting matter under 
the influence of an electric field E and the electro- 
motive force q X H that is set up by motion relative 
to a magnetic field. The process of electrical con- 
duction in gases in the presence of a magnetic field 
is really somewhat involved. Basically, a force acts 
on any charged particle that is moved across a magnetic 
field. In our case the charged particles are the elec- 
trons and 1ons of the gas, and, in relatively dense gases, 
they are constrained to move with the gas velocity 
q. The electromotive force on them causes them to 
drift across the flow, their motion being impeded by 
frequent collisions with the other particles, and a net 
flux results that is described by the second term of 
Eq. (7). The first term represents a similar drift of 
charged particles, but here caused by the electric 
feld E. 

The symbol E represents the total electric field 
strength (e.g., volts per unit length) due to external 
effects, such as the field set up by charged boundaries 
of the flow, and internal effects —i.e., separation of 
charges or polarization.’ E is measured in the frame 
relative to which q is measured. In our approximation 
(gy < c¢ = speed of light), H is not a function of q. 
One can think of the term q X H as representing a 
tiny generator or source of e.m.f. at any point in the 
moving fluid; the effect of its e.m.f. is divided between 
driving a current through the gas and separating the 


positive and negative charges. If no current can flow, 


because of the presence of insulators, for example, the 


whole e.m.f. must be opposed by the field E—i.e., by 
charge separation. More generally, the current density 
J cannot always be distributed like q X H, so that a 


net field E will generally result, in magneto-aero- 


dynamics. A calculation shows that the distance 


through which the positive and negative charges are 
moved in a gas to set up any reasonable field strength 
E is extremely small. It is a good approximation, 
therefore, to consider any infinitesimal sample of the 
gas to be neutral in charge, thus the field E appears 
to be produced by doublets—1i.e., the gas is polarized. 
This is the approximation usually made for dielectrics, 
but it is also appropriate here. It does not imply 
that E is zero or negligibly small. 

The importance of the quantity E is illustrated by 
combining Eqs. (6) and (7). The first two right-hand 


terms combine as follows 


qjxH+ (fo) =qjxH+(E+qXH) jl 
E-j 


(NS) 


where the rule of signs of a sealar triple product has 
been used. Thus the net magneto-aerodynamic con 
tribution to the energy is seen to involve directly 


the electric field E. 


Maxwell's Equations 
The introduction of Ohm's Law has related j to q 
and H, but only with the involvement of a new un 


known, E. 
terminate by means of Maxwell's field equations 


The set of equations is rendered de- 


barj curl H 9) 


curl E —(OH Or) 10) 


As before, we have assumed the fluid to have the same 
magnetic permeability as empty space and have used 
electromagnetic units. We have neglected ‘‘displace 
ment currents” in writing Eq. (9).! 


RELATIVE MAGNITUDES OF TERMS 


Having written down the equations of motion, we 
are now able to investigate the relative magnitudes of 
the terms that appear, and in this way make some 
prediction of the importance of magueto-aerodynamiuc 
There is an equally important reason for 
namely, that the set of 


effects. 
making this investigation 
equations is formidable and will require simplification. 
Aerodynamicists have had little success in trying to 
solve the equations of compressible viscous gas flow 
without drastic approximations, and it is even less 
likely that they can be handled with the added com- 
plexity of electro-magnetic terms. It will be necessary 
to approximate in the same rational way —i.e., to 
retain the most significant effects while discarding 
others in appropriate areas and types of flow. One 
may ask in this spirit, will there be any application 
of the magneto-aerodynamics of ideal fluids? Will 
the boundary-layer approximation still apply? Will 
the same approximations that the astrophysicists 
make in magneto-hydrodynamics (which will be 
described below) be permissible in aeronautics? 

To begin this part of the study, we need, first of all, 
some estimate of the probable magnitudes of the 
conductivity o and the magnetic field strength H. 
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Fic. 1. Calculated conductivity of air and of air seeded with 
0.1 per cent potassium (by weight) as functions of temperature 
for two altitudes. Corresponding Mach Numbers J giving 
approximately these stagnation temperatures are indicated. 
(From references 4 and 5.) 


In the light of Ohm's Law, Eq. (7), we might also ask 
about the probable magnitude of E, but we shall 
answer this question by specifically assuming that it 
has the same order of magnitude as q X H. This 
assumption implies that the only electric fields present 
in the air are those induced by magneto-gasdynamic 


phenomena. 


Conductivity 


As has already been pointed out, electrical conduc- 
tion in a gas is accomplished by a drift of charged 
particles. Conductivity is therefore increased by 
increasing the number of charged particles and reduced 
by increasing the probability of their collisions with 
impeding particles—i.e., atoms, molecules, and charged 
particles. 

Let a denote the “‘fractional ionization” of the gas, 
defined as the number of original particles ionized. 
Thus a is nearly zero for air at room temperature, and 
a = | represents a completely ionized gas, in which 
each original particle has given up one electron. Al- 
though @ can exceed 1, the state of complete ionization 
represents a relatively large value for most applications, 
attained only at very high temperature. When a is 
small, so that charged particles in the main collide 
with neutral molecules, one arrives at the result 


a~al-v? (11) 
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This regime is probably typical of aeronautical cases 
but, in general, relation (11) is misleading in form sing 
a itself is strongly dependent on the temperature 7 
Thus even in this regime, o varies rapidly with 7. 

At the other extreme, when collisions between 


charged particles are dominant, the conductivity 


becomes independent of a, and relation (11) is re 
placed by 


o~ TY (12 


There is an important case, however, in which the 
strong dependence of o on 7 may be minimized, 
namely the case where conductivity is achieved by 
seeding. This means that an easily ionizable substance, 
such as potassium or sodium, is introduced into the 
air. The air temperature may be high enough to 
ionize completely this seeding substance, though it is 
insufficient to ionize the air itself appreciably. The 
result is then to remove the temperature dependenc 
of a on 7’ in Eq. (10), making o relatively insensitive 
to temperature. The effects of seeding with various 
substances have been studied by R. Rosa in a Cornell 
University thesis. * 

Some data on conductivity of air with and without 
seeding, from references 4 and 5, are presented in Fig. 1, 
The Mach Numbers spotted along the curves are 
only approximate and are intended to orient. the 
reader as regards the relationship of stagnation tem- 
perature to flight speed. It is seen that o for atmos 
pheric air does not exceed about 0.1 mho/cm. even at 
the temperature corresponding to satellite speeds. 
On the other hand, one may obtain values near 1.0 
mho/em. by seeding, even at Mach Numbers con- 
templated for manned aircraft. 

For purposes of orientation, it might be helpful to 
list the conductivities of other similar substances, in 


the same units. For example, 


Copper 600,000 mho/cm. 
Mercury 10,800 mho/em. 
Salt water (saturated at 25°C.) 0.25 mho/cm. 


Pure water 0.0000002 mho/cm 


On the basis of Fig. 1, we shall adopt the value 
o = 1.0 mho/em. as a standard for certain numerical 
comparisons below. This value has the virtue that 
it is easy to change in any numerical result by appro- 


priate multiplication. 


Magnetic Field Strength 


For the magnetic field strength //, on the other 
hand, there is no completely rational basis of estima- 
tion. Nevertheless, some milestones can be identified. 
A field strength of about 20,000 gauss saturates all 
iron core at room temperature. It is not certail 
that iron cores can be utilized in aeronautics in view 
of their weight; still, we can be sure that if they are 
used the saturation value will be the maximum field. 
This would seem to set only a very high upper limit 
even though in principle an air-core magnet could be 
made stronger. 
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Another indication of practical field strengths may 
be suggested by the value available in a permanent 
magnet. This runs to perhaps 1,500 or 2,000 gauss 
for good modern materials. 

1 s 

It seems probable that in aeronautical applications 
magnetic fields will be produced by air-core electro- 
currents are generated magneto- 
gerodynamically. The 


can be estimated in terms of vehicle dimensions, but 


magnets whose 


dimensions of such solenoids 
since the number of turns will be determined by such 
matters as weight and resistance, one still does not 
have an estimate of /7. In some simple computations 
carried out below we shall adopt rather arbitrarily the 
value [/ 1000 gauss. Once again, its main virtue 
may be that it affords easy correction to other values 


for those who do not concur. 


Dimensionless Ratios 


In fluid mechanics it is useful to compare the magni- 
tudes of dynamic (inertial) and viscous effects by 
forming their that geometrically 
similar patterns of flow are being compared. The 
result, of course, is the Reynolds Number. It often 
has values of hundreds of thousands if the typical length 
is taken to be body diameter or wing chord, and one 
concludes that the viscosity can be ignored through- 
It is only a slight refine- 


ratio, assuming 


out most of the flow field. 
ment of this procedure to introduce the boundary-layer 
thickness as the typical length (or, conversely, to 
inquire what typical length is required to make vis- 
cosity important) and to show thereby that inertial 
and viscous effects are comparable in magnitude in 
a boundary layer. 

Let us undertake to apply this method to compare 


the e.b.f. with the dynamic force—i.e., the first and 


last terms of Eq. (4). In order of magnitude, the 
ratio is 
e.b.f. jH o VIL ¥ 
a = (hp) 
dynamic force (1/2) pV?/L (1/2) pV? 


where |’ is the typical speed, say the flight speed, and 
L the scale-determining length. Here we assume, as 
mentioned above, that the electric field E has the same 
order of magnitude as q X H;; in other words, that the 
current density j is of orderoV//. If L is put equal to 


an aircraft or missile length (or diameter), one is 
comparing the total force per unit cross section of the 
flow, which can be applied over this length, to the 
dynamic pressure of flight. Some values are plotted 
where altitude and Mach Number are the 


It appears 


in Fig. 2, 
variables and the assumptions are stated. 
that the e.b.f. will be big enough to affect the flow 
appreciably at high altitudes and high speeds. 

Naturally, if we form the ratio of body force to 
viscous force we obtain the product of this number 
and the Reynolds Number: 


e.b.f./viscous force = 


(e.b.f./dynamic force) K Re (14) 
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where Ri p\L/u. Thus, this ratio is many times 
as large as the one plotted in Fig. 2, except in boundary 
layers, wakes, etc., where the significant Reynolds 
Number is much smaller. The immediate conclusion 
is that an inviscid theory of magneto-aerodynamics 
is generally useful but must be interpreted in the light 
of boundary-layer theory. In other words, boundary 
layers and similar viscous regions have their common 
significance, in magneto-aerodynamic flows as well. 
The number calculated in Eq. (14) becomes twice 
the square of the Hartmann Number® when the signifi 
cant length is the diameter or breadth, D, of a pipe or 
channel: 
pl*| X 


2HTa°* lo VIED /(1/2 


lp |D MK} 2all*D be 15) 
A third interesting ratio is obtained by comparing 
Let us compute 
the 


In doing so, we 


electrical currents rather than forces. 
the 
necessary to produce the field H. 
shall assume again that the current density 7 is of order 
oVH. The current to produce H, on the other hand, 
depends on the dimension and the number of turns 
Dimensionally, 


ratio of motion-induced current to current 


in the coils of the electromagnet. 
this current ~ ///L, but a large factor, the number 
involved. Ignoring this omission, 


of turns, may be 


one has 


og VL 


a VH/ (T/L) 
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Fic. 3. Sketches of proposed device using magneto-aerodynamic 
effects to decelerate or control re-entry missiles 


Because of its resemblance to a Reynolds Number, 
this ratio has been called the “‘magnetic Reynolds 
Number.’” It can also be interpreted as the ratio 
of motion-induced magnetic field to applied magnetic 
field. 
generated involve their own magnetic fields, in ac- 
This can be thought of as a 
perturbation of an field. The 
the perturbation would follow from 


In this sense, one recognizes that the currents 


cordance with Eq. (9). 
applied magnetic 
magnitude of 
Eq. (9) by a process of integration equivalent to the 
calculation of the velocity field of a vortex; the result 
would certainly be of magnitude oI//L. Thus the 
ratio of induced to applied magnetic field strengths is 
equal to the magnetic Reynolds Number. This may 
be a more satisfactory interpretation of R,, since it 
does not involve the number of coil-turns mentioned 
above. 

In astrophysical and similar problems involving 
either high conductivities and or very large dimensions, 
the magnetic Reynolds Number may be so large as 
to show that the motion-induced magnetic field is 
many times as large as the applied field. In such 
cases, the applied field may sometimes be neglected 
This leads to the concept of a fluid of infinite 
For such a fluid the differential Eqs. 
For example, Eq. (7) 


entirely. 
conductivity. 
(3) to (10) are much simplified. 
requires that E + q X H be zero, since j must remain 
finite. It can be shown that under these conditions 
the magnetic field is convected with the fluid. Since j is 
finite, the Joule heating disappears. There is a con- 
siderable body of literature pertaining to this case, 
especially because of its astrophysical importance, as 
has been mentioned. 

Returning to aeronautics, let us investigate the 
probable values of our magnetic Reynolds Numbers. 
Suppose that o = 1.0 mho,/em. and L 1.0 m. as 
before, then 


KR = OVE 


10~-°I’(m. ‘sec. ) (17) 


where the factor 10~° has been introduced to relate 


the present units to electromagnetic units. Thus 
R,, is less than 1 10 even at satellite speeds, and the 


approximation of infinite conductivity seems not to 
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be appropriate for aeronautics. On the other hand. 
considering the missing factor related to coil-termis, thy 


values of R,, in Eq. (17) may be large enough to mak 





an aeronautical system self-powered—1.e., the currents } 
induced by flow of seeded air may be sufficient to power 


the electromagnets that produce the field. 


EXAMPLES 


In this section we propose to present a few « xamples 
of magneto-aerodynamic flows, chosen to_ illustrate 
some of the effects discussed above and some of the 


difficulties of solving the equations. 


Deceleration of Flow | 
Consider the body force j X H with the current 


density given by Ohm's Law (Eq. (7) 


j XH o\E X H+ (q X H) X H} IS 
The force o(q X H) X H always has a component 
opposing the flow—-1.e., tending to decelerate the flow 

A possible application of this effect occurs in missile 
re-entry where it is sometimes desirable to increase 
accomplished as 


the missile drag. This might be 


sketched in Fig. 3. A missile nose is shown enve'oped 
by a layer of hot, ionized air. 
magnetic field H be provided and that the current 

generated, j, be made to flow through electrodes, as | 
sketched. 
would activate the electromagnet. 
as is mentioned again below, the electric field strength 
due to polarization, E, would be small if substantially 
two-dimensional conditions occur in the ionized layer 


It is proposed that a 


Inside the missile, the current generated 
In this situation, 


Thus the vectors representing the e.b.f. are shown in 
Fig. 3 to be opposite to the flow direction everywhere. 

The result would be a considerable drag on_ the 
missile; the force would be applied to the structure of 
the electromagnets. From an energy viewpoint, one 
recognizes that kinetic energy is being converted to 
electrical form and then being dissipated in Joule 
heating of both the ionized air and the coils of the 


electromagnets. 





y } 


LLLKLLLLLLLLLL LLLLL/ 
UU(y) 





x 





INCREASING Fi, 








sh CM hhhhts 
A, 


Poiseuille flow of a conducting liquid, I 
field applied in plane of flow 


Fic. 4 


th 
ste 


P 


pr 
col 


din 
the 


pre 


wh 


an 


r hand, 
"ms, the 
-O make 
urrents 


o) pr Wer 


‘amples 
uStrate 
of the 


“‘urrent 


(18 


ponent 
e flow 
missile 
icrease 
ed as 
coped 
that a 
urrent 
eS, as 
erated 
ation, 
ength 
itially 
layer. 
wn in 
lere. 
1 the 
ire of 
, one 
ed to 
Joule 
f the 


netic 





| 


PROSPECTS FOR 


We have not made a detailed study of this device. 
\ device based on similar principles was studied in 
detail by R. M. Patrick in a Cornell thesis.° Patrick's 
objective was to produce lift rather than drag, so that 
he considered the device at angles of attack. 

In ipplications of this kind it is clear that the 
distribution of the e.b.f. across the ionized layer is 
nonuniform because the distribution of q is nonuniform. 
In view of the usual boundary-layer profile, it appears 
that the decelerating force must increase in magnitude 
with distance from the body The result 
must then be a tendency to alter the boundary-layer 


surface. 


profile in a sense that would diminish both heat transfer 
These tendencies have been recog- 


mentioned 


and skin friction. 


nized bv Rossow,” whose work will be 


again in connection with the next example. 


Two-Dimensional Channel Flow (1) 


A classical case of magneto-hydrodynamic flow is 
the one sketched in Fig. 4, namely two-dimensional 
steady flow between infinite, parallel, insulating walls 
Poiseuille flow) with a uniform magnetic field H, 
applied normal to these walls. 


by Hartmann! for an incompressible fluid. 


This case was treated 
Under 
the hypothesis stated it can be assumed that all physical 
quantities are functions of y only (excepting, of course, 
the pressure P, which can vary linearly with x). 

Eq. (10) implies div H 0, and, in the present 
problem, this means that the y-component of H is 
Similarly, Eq. (9 


constant and equal to JJ». states 


div j 0, and this means that the y-component of 
jis constant; it must therefore be equal to zero since 
the top and bottom walls are insulators. For incom 


pressible flow, then, Eq. (4) becomes 


r u(d7u/dy*) — 7H 19) 
where P denotes the constant pressure gradient Op Ox 
and 7, is the z component of 4s. 

The equation for the z-component of momentum 
now states that no current flows in the x direction 
1e.,j has only a z-component. The remaining equation 
of momentum then becomes 


Op oy if 3 (20) 


where // Further conse- 


quences of the single-component nature of j follow 
namely, that //., the s-component of 


is the x-component of H. 


from Eq. (9 


H, must be constant and 


(?1) 


inj. = —OII, dy 


According to the hypotheses of the present problem, 
the constant /7, must be equal to zero, but in any case 
it does not affect the dynamics of the flow—i.e., the 
solution is the same if the applied magnetic field has 
a Z-component. 

An analysis of the electric field vector E is accom- 


plished by means of Eq. (10), which states curl E = 0 
in this problem. This leads to the conclusion that 
both the x- and z-components of E are constants. But 
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), now states that the x-component 


Ohm's Law, Eq. (7 
is zero, that the y-component iS equal to u//., and that 


j. is given by 


J ojikE ull t 22 


where /y) is the constant z-component of the electric 


field, determined by boundary conditions at the ex- 


tremities of the flow. Combining Eqs. (19) and (22), 
we have the differential equation for u(y 
u(d*u dy") — ofly*u P+ ok Il 23 
the solution of which is 
P+ chill | cosh \ a llyy | 
ul - ; 24 
oll, { cosh V o/s Iya / 


where a is the channel half-width as shown in Fig. 4 


The familiar parabolic profile of Poiseuille flow is 


obtained from Eq. (24) by carrying out the limiting 
process ¢ — O. 

It is interesting that in this case the velocity profile 
shape does not depend on /), although the current 
density certainly does, as well as the Joule heating. 

The profiles in Fig. 4 are sketched from Hartmann's 
show how the profile of 
flow is flattened by 


The flow of fluid for given pressure 


paper. They parabolic 


Poiseuille the magneto-hydro 
dynamic effect. 
gradient Op Ow is reduced. 

The solution is completed by 
(21) for &.. J 
Also, one can combine Eqs. (20) and (21 


integration of Eq. 


being given by Eqs. (22 


and (24 
and integrate, 
vielding 

Sw)/7,* + function of x | ma 


Sax)H,? + Px 4 


p —(] 


—(1 constant f 


This type of result has led to the terminology 
“magnetic pressure’ for like [7,?/Sr. It is 
clear that in the present problem the pressure gradients 
the interaction of the 


terms 


across the channel arise from 


induced current with the accompanying induced 


magnetic field. 

The alteration of profiles by magneto-hydrodynamic 
effects in this problem suggests that a similar effect 
may occur in the analogous boundary-layer problem. 
This leads to the problem studied by Rossow* (the 
case of two-dimensional incompressible flow over an 
infinite flat plate with a normal magnetic field). 
Rossow’s results show changes of profile analogous to 
those of Hartmann, as predicted above. Rossow 
proceeds to calculate the alteration of temperature 
profiles, accounting for Joule heating, heat conduction, 
that heat 


and viscous concludes 


transfer to the wall may be appreciably reduced by 


dissipation. He 


magneto-hydrodynamic effects. 

Rossow’s paper involves two interesting extensions 
of this problem. One is a calculation of the analogous 
case with magnetic field fixed to the undisturbed fluid 
stream. This represents, for example, the case of a 
plate moving through the atmosphere and the earth's 


magnetic field. The second is an investigation of the 
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Poiseuille flow of a conducting liquid, II Mag- 
Top and bottom 


Fic. 5. 
netic field applied normal to plane of flow. 
walls of channel are conducting. 


effects of variable o. 
that o varies rapidly with temperature; accordingly, 


It has already been pointed out 


Rossow writes 


a = constant |(#.. — u)/u.] (26) 


where u Although Eq. 
(26) is rather arbitrary, it should give an indication 
of the effects 
when conductivity drops off to a negligible value at 


denotes the stream speed. 


reduction of magneto-aerodynamic 


the edge of a boundary layer. This is not quite the 
case sketched in Fig. 3, to be sure, for there the ionized 
layer is bounded by a shock wave, and the approxima- 
tion of constant o does not seem so objectionable. 


Two-Dimensional Channel Flow (II) 


A slightly different modification of Poiseuille flow 
is sketched in Fig. 5; this represents the same case 
as in Fig. 4, but with the uniform magnetic field 
directed along the z-axis—1.e., normal to the plane of 
the flow. This problem has the virtue of illustrating 
how an electric field E is induced by the motion; it 
also shows how relatively complicated a seemingly 
simple flow may become with magneto-hydrodynamic 
effects present. We shall treat it only in the approxi- 
mation of first-order perturbations. Specifically, we 
assume that the magnetic field vector is made up of 
the uniform applied component Ho plus a much smaller 


perturbation h. 


H = Hy + h, h < Hy 


(27) 


The principal effect of the magnetic field is then the 
production of a current-density component j, in the 
y direction, and the x-component of Eq. (4) becomes, 
approximately, 


P = 0p/Ox = pu(07u/dy?) + j,Ho (28) 


Now Eq. (9) requires div j = 0. As in Hartmann’s 
problem above, we postulate that conditions do not 
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vary with the streamwise coordinate x, but we cannot 


also require independence of z. For the existence 0 
a component j, implies, according to Eq. (9), that 


the induced magnetic field varies with z—1.e., 


Ie Oh, Oz ~ 0 At) 


Thus the flow can be only quasi-two-dimensiona] 
Eq. (29) implies that a streamwise component h, js 
that it This 


signifies that the magnetic lines of force are actual 


induced and varies with gz. Variation 


curved whereas in our approximation we consider 
them to lie straight across the flow in the z direction 
The current-density component j, 1s then approxi 
mately constant and the term j,//) in Eq. (28) repre 
sents an additional pressure gradient. The solution 
of Eq. (28) is the parabolic profile of Poiseuille flow, 
only the flow quantity for given pressure gradient 
being changed. 

The approximate evaluation of j, follows from 


Eq. (7): 
aD ( 9 
Fito, ull; (50 


j, *~ constant 
If the top and bottom channel walls are insulators, 
or if the circuit connecting them is open, the onh 
that the 


circuit is closed, as sketched, by a conductor of re 


solution is J, 0. Suppose, however, 
sistance R per unit length and breadth of channel 
Integrating Eq. (30) with respect to y, one has, for 
constant oa, 


2aj, = —ocAV — oll | udy (3 


where the top and bottom walls are at y £a and 


At denotes the voltage difference between these 
walls—.e., 
%a 

| E,dy = —V(a) + U(—a) — 80 (32 
oe =—— 


But the same current must flow through the exterior 
circuit—.e., 


uP AU/R 33 


Thus, combining Eqs. (31) and (32), 


a 


—Hy udy 


. (34 
jt: = 
R + (2a/¢c) 
The distribution of electric field is given by 
Ey (jy/o) + ulo (30 


as sketched in Fig. 5. 

It is interesting to note that the current density 1s 
proportional to the quantity of flow through the 
channel. We are informed that flowmeters operating 
It is also interesting 
to consider the energy balance in this case. For this 
The only alteration 


on this principle are available. 


purpose we use Eqs. (6) and (8). 
of Eq. (6) for the case of an incompressible fluid 1s 
(To be sure, it 


> 


the deletion of the term pdp~'!/Ot. 





net 
the 
mo} 
pre 


con 


Cannot 


Cnce of 
)), that 


WAN) 


sional 
't h, 18 
iiation 
Ctually 
onsider 
ection 
DP proxi 
repre 
ution 
~ flow, 


‘adient 


from 


(30 


lators, 

only 
t =the 
of re 
annel 
s, for 


1 and 
these 


‘erior 


99 
(30 














PROSPECTS FOR M 
is not common to use the concept of stagnation en 
thalpy for an incompressible fluid because the total 
+ pg*/2, is equally convenient and there is 
and 


head, 
ordinarily no 
In magneto-hydrodynamics, however, there 


coupling between temperature 
velocity 
‘s such coupling by virtue of Joule heating.) 

Thus, the rate at which stagnation enthalpy is 
This rate is j-E, 


By integration 


being increased is given by Eq. (8). 
or in our present approximation, j,/,. 
of this quantity from wall to wall, this rate is easily 
found to be 


« 


—j,"R (336) 
In other words, energy is being extracted from the 
fluid at a rate equal to the Joule heating of the external 
Another check on the 
namely, that the total work done 


resistance. energy balance 
can also be made 
by the pressure gradient Op/Ox appears in heating of 
the fluid and the external resistance. 

The induced streamwise component /, of the mag- 
netic field is not negligible, as Eq. (29) indicates. Thus 
the body-force component j,i, in the equation of z- 
momentum is also appreciable, and there must be a 
pressure gradient Op Oz of the same magnitude. To 
complete this solution by calculating to a_ better 
approximation all of the components set up by the 
motion is a complicated process. This problem has 
been treated, though not by the perturbation process 
suggested here, by Shercliffe.* 

Although this sample does illustrate a case where 
the electric field is not originally known and must 
be calculated, it is a particularly simple case in the 
approximation employed here because the directions of 
the vectors j, E, and q X H are determined by the 
hypotheses of the problem. In general, in two- and 
three-dimensional problems, this is not the case; in 
fact, these vectors may have three different directions. 
The technique then employed is to eliminate j from 
Eq. (4) means of Eq. (9), and to obtain an auxiliary 
partial differential equation by combination of Eqs. 
and (10) 


(7), (9), V1Z., 


(1/4) curl? H = o}—(O0H/Of) + curl (q X H)j (37 


One-Dimensional Channel Flow —Acceleration of a Gas 


Let us consider here a problem of more immediate 


aeronautical interest--namely, the acceleration of 
a gas flow in a channel such as a nozzle or tail-pipe. 
If the Reynolds Number is large and the channel 
cross-section .1 changes relatively slowly, one can 
neglect viscosity and treat the problem by the familiar 
methods of one-dimensional channel flow. In other 
words, one can consider the flow speed u to be a function 
of the coordinate x measured along the channel length 
only. Thus it is really an average speed. The same 
approximation can be made for the pressure p, density 


p, etc. The equations of motion, Eqs. (3), (4), and 


(5), become, for steady flow, 
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Continuity: pu-l constant 38 
Momentum: puu'’ + p’ ojk — ull i 39 
Energy: puc, 1” + urn’ oj\E — wHtE 10 


where the prime denotes differentiation with respect 
to the coordinate wx. 

We recall that, in ordinary flow problems without 
magneto-aerodynamical effects, Eq. (40) becomes a 
statement of conservation of stagnation enthalpy. 
Eq. (38), together with this statement, leads to the 
familiar results of one-dimensional gasdynamics such 
as nozzle flow and choking. These appear as special 
cases in what follows. 

Suppose now that not only /(v) but also the field 
strengths E(x) and //(x) are given. Our point of view 
is therefore somewhat different from that adopted in 
our earlier examples, for here we do assume that an 
exterior to the 
this field 


inagneto-aerodynamically by 


electric field is applied by a device 


channel. We imagine, for example, that 
might be generated 
motion of an aircraft and applied across an interior 
duct or tailpipe. Thus, within reasonable limits, we 
assume both / and // to be under control of the de 
signer. 

In general, for given -1(x), E(x), and J/(x), Eqs. 
(40) cannot be integrated, but the acceleration 
uw’ and the derivative of the Mach Number, J/’, can 


(The Mach Number that appears here 


(3S) 


be calculated. 


is the conventional ratio v/a, where a denotes V yp p, 


1 Cafe The results are 
; | v4 oll* | 
ul : : Je (ud — U3)(u My (4] 
M2 —1 lA p J 
and 


1+ [(y —1)/2)1r 


VW’ os 
iM? — 1 
(uM ., oll’ { 
A (4% — U3)(u u (42 
(4 ap { 
where the significant velocities 1“, #2, and u; are defined 
as follows: 
My I(y y\(E/H 13 
u 1 + yi*) [2+ 4 1) M2] ha (44 
u E/H 15 


The terms involving -1’(«) in Eqs. (41) and (42) are 
familiar to gasdynamicists; they state that the velocity 
and Mach Number decrease in an expanding channel 
at subsonic speeds and increase in an expanding 
channel at supersonic speeds. 

The signs of the magneto-gasdynamic terms, on 
the other hand, depend on the magnitude of u relative 
to the significant velocities ™;, m2, and w;. In par 
ticular, u’ can be made positive even in constant-area 
channels at supersonic speeds provided only that wu 
lies in the interval between uw, and u;. Furthermore, 
\/’ can be made positive in a constant-area supersonic 


channel if uw. < u < u;. This is in sharp contrast 











JOURNAL OF THE 





AERONAUTICAL 





CODE: 
\X\\ du“fax > 0 
4 / f/f aM/dx > O 
Yn or A AsympTOTE AS X-+00 
—O “ TUNNEL ” 
{ CHOKING 























ba | 
O 
O ‘| | 2 M 
IID 
Fic. 6. Magneto-gasdynamic acceleration of one-dimensional 


flow in a constant-area channel. Regimes of positive accelera 
tion and increasing Mach Number are designated by cross-hatch 
ing as shown. 4, tM, and uw; are the “significant speeds’’ de- 
fined in this paper. Arrows denote typical behaviors of flows 
in the various regimes; arrowheads show direction of increasing 
x (downstream). 14, 1B, ete. denote areas of different behaviors 
as indicated 


to the effects of pure heat addition, which always 
decelerates a supersonic constant-area flow and drives 
its Mach Number toward 1. 

A summary of the various possibilities for constant- 
area flow is presented in Fig. 6, where the signs of w’ 
and \/’ in a u, AM plane are indicated. There are 
eight distinct regions in the u, J plane, in which 
different behaviors of u(x) and A/(x) can be predicted 


from Eqs. (41) and (42). These are identified in 
Fig. 6 as follows: 
 w>] (II) M<1 
(A) uw <@ (A) uw< 4 
(B) t<u<u (B) m<u< us 
(CC) m< 2 < mw, (C) te <a < 
(D) u<my (D) u< ue 


In general, a ‘barrier’ exists at 1J = 1, where the 


phenomenon of choking occurs: the denominators of 
both Eqs. (41) and (42) disappear. 


sents choking is evident, from the fact that the curve 


That this repre- 


of u versus x turns back itself and becomes 


double valued after reaching the barrier. 


upon 
Thus the 
sonic speed can never actually be attained in such a 
channel, and the phenomenon would have to occur, 
if at all, at the open end of a channel. 

There are, however, two exceptional cases, denoted 
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“tunnels” in Fig. 6. These are two one-way tunnels 
through the barrier at sonic speed. For these special 
(11d IT Db the 


numerators and denominators of both Eqs. (41) and 


combinations of parameters and 
(42) disappear simultaneously, and, as can be verified, 
finite and continuous deceleration and acceleration, 
The to these 
tunnels have been plotted in Fig. 6, as they always 


respectively, occur. curves that lead 


separate choking from nonchoking channels and 
represent the special combinations that permit smooth 


passage through the sonic speed. 


Maximum-Acceleration Constant-Area Duct, Supersonic 
Flow 
According to Eq. (41), the acceleration of a super- 
sonic stream in a constant-area channel can be made a 
maximum by choosing a certain value of FE J] —ic., 
the value that maximizes the product —(u — u;)X 
This value is 


EH 


(a — 2). 


[((2y — 1)/2(y — 1) ]u (46 


The maximum acceleration is obtained by substituting 
this value in Eq. (41); the result is 


Umar’, = [1/(y — 1)](oH?/4m)[M2u/(M2 — 1)] (47 


where m denotes the constant pu (the mass flow per 
unit cross-section area). 

As an indication of the possible velocities to be 
achieved magneto-aerodynamically, we have integrated 
Eq. (47) for constant y. In other words, we assume 
that /:/// is varied with x so as to maintain u’ always 
at its This 
first, integration of ./ as a function of u, which can 


maximum value. calculation involves, 


be carried out using Eqs. (41), (42), and (47). The 
results are 


M? : 
u constant X ‘ (48 
[(2y + 1) y¥i—- @ 
and finally, 
(1 m) | olds = 
imy — bier 1 u 


In — (49) 


) 
2y+1 t 4 Mi 
yMy? — (27 + 1) “) is aT 
(Qy + l)yMy? (" Ty 
Here the quantities with subscript 0 denote values at 
If ofl? is constant, 
Eq. (49) permits an easy calculation of the channel 


Eq. (49) is 
plotted in Fig. 7 for several values of 1/o, for y = 1.4. 


an initial station of the channel. 
length x — x» for velocity ratio u/ up. 
It is interesting to note that the curve of u/u against 
& becomes practically independent of 1/9 for Ao > 5, 
making w/w a function of & only. 

The corresponding Mach Numbers are not plotted, 
but could be determined from Eq. (48). Since 4 
always increases with x in this case, the form of Eq. 
(48) that toward the value 


indicates tends 
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PROSPECTS FOR 


24 1.4, J 1.65. Numerical 


calculations show that this approach is very rapid; 


V(2y + 1) y, or, for ¥ 


\f is nearly equal to 1.65 when w up reaches 3 for any 
initial Mach Number greater than 1. It is clear from 
this result that there is considerable heating of the gas 
in a typical case, and that appreciable further increase 
{ flow speed might be achieved by expanding in a 
supersonic nozzle. 

The present results (Eq. (49) and Fig. 7 
be used when -& — u// is greater than the ‘“‘breakdown’ 
This critical potential (where 


must not 


potential of the gas. 
arcing occurs and therefore the assumptions of the 
present studies are violated) varies directly with density 
in the range of densities anticipated here; i.e., the 
restriction on our results may be written as 


E— ull < Kyp (50) 
where A, has the value of about 2 X 10% volt cm.” gm. 
for air.* Thus the restriction can easily be expressed 
by combining Eq. (46) and inequality (50) to give a 
maximum permissible value for u, above which Eq. 


16) will lead to a violation of (50). It is 
u< V2(y7 — 1)10A,m/H (51) 
where /7 is measured in gauss, m in gm. cm.~* sec. ~! 


ind w in em. sec. This does not mean that speeds 
greater than this cannot be achieved, but only that 
they cannot be achieved under the maximum-accelera- 
tion conditions of Eqs. (46)-(49). For greater velocities, 
acceleration might have to be accomplished under 
conditions of constant potential H — u//, for example. 
We are studying further the consequences of this 


limitation in practical cases. 


CONCLUSIONS 


that electro- 


important in 


The studies reported here suggest 


magnetic effects become aero 
dynamics, especially if seeding can be used to increase 


Here we 


may 


the conductivity of the air. understand 


aerodynamics’ to mean the realm of flight where 


the continuum approximation is _ relatively useful. 


for example, possibilities of decelerating 
skin and 


There are, 
re-entry missiles, controlling the friction 
heat transfer in boundary layers, and accelerating or 
decelerating channel flows magneto-aerodynamically. 

The last mentioned of these possibilities appears 
especially interesting to the present authors, since it 
permits some significant departures from the restric- 
tions of conventional gasdynamics, such as acceleration 
of supersonic flows in straight tubes. 

There are, however, some formidable practical 
obstacles in the path of engineering applications of 
the principles mentioned here. One of these may be 
the complexity of the equations of motion, since this 
will hamper the theoretical studies that might lead to 
Significant and efficient applications. More serious, 


“We can expect the same value to apply for seeded ait 


MAGNE 
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Fic. 7. Magneto-gasdyvnamic channel flow Flow speed 


ratio as function of the channel-length variable € for a maximum- 
acceleration channel of constant area for several initial Mach 
Numbers 


however, is the problem of the development of magnets 
powerful enough to provide the required fields and yet 
light enough for use in flight. Techniques for render 
ing gases conducting by means of seeding must also 
be made practical. 

Another engineering problem of particular signifi 
cance seems to be posed by the need for electrodes in 
all the magneto-aerodynamic devices suggested in the 
present paper. Hartmann! mentioned the high surface 
resistance of electrodes in contact with mercury in his 
early experiments. Brogan'’ found a similar phenom 
enon in a shock tube, where the overall resistance of 
a circuit including cold copper electrodes in contact 
was times the calculated 


ionized many 


The solution of this rather 


with argon 
resistance of the gas alone. 
basic problem is certainly not obvious at present. 

Thus, successful achievement of magneto-aero 
dynamic effects in flight must be preceded by con 
siderable research and engineering development. The 
possible benefits to be obtained seem to justify such 


efforts. 
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Aeroelastic Instability at High Mach Number’ 


JAGANNATH P. CHAWLA* 


Hughes Aircraft Company 


SUMMARY ‘a 


A brief discussion of the limitations of the piston theory is given 


Using Lighthill's approach, the aerodynamic forces on oscillating Ti = 9 f 
airfoils are derived in a slightly different form to conform to the cae ; 
NACA practice and the familiar nondimensional flutter coeffi : 
cients are tabulated for general double and single wedge and 7 
parabolic biconvex profiles. Extensive parametric studies are ] 
included for preliminary design to aid in estimating the effect on 7 
flutter speed and frequency of airfoil thickness and profile, initial 
angle of attack, modal shape, changes in elastic and inertial prop 
erties of the wing, and altitude of flight a 
SYMBOLS a 
a = velocity of sound =O 
Adé = 1, 2.3 = constants defined in Eq. (23) and 
Table 1 a 
r = airfoil chord 
h = subscript for altitude 
h = amplitude of airfoil oscillation in y 
vertical translation 
h(t = vertical oscillation of airfoil (see 
Fig. 2 Vix, t 
(; = 1,23) = integrals defined in Eq. (29 
® = moment of inertia of wing about V(x) + Vix, t 
elastic axis per unit span 
k = reduced frequency = (1/2) (we/l ' 
K,, = wing stiffness in flexure Z 
1 = wing stiffness in torsion 
L = lift force on the airfoil 
L;, M; = nondimensional flutter coefficients ” 
= 1,2, 3,4 defined in Eqs. (15), (16), and (17 a 
Lo = [.(Mk 
L = L,(Mk) . 
Vi, = M, (Mk) 
ml = wing mass per unit span p 
V = Mach Number of undisturbed stream ? 
= U/a 5 
mM = moment (positive nose up) per unit 
span about the axis of pitching 
oscillations (see Fig. 2) - 
0 = subscript for sea level 
p = pressure at point on airfoil surface s 
Por Pore = pressure, density, and velocity of ” 
sound, respectively, in undis- 
turbed stream ‘ 
Ap = pressure on lower surface minus 
pressure at same point on upper » 
surface $ 
POR = constants defined in Eq. (283) and 5 
Tablel — - 
p 
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ratio of radius of gyration in pit 


to semichord 


(40 q/me?)!/* 

constants defined in Eq. (23 
Table 1 

time 





velocity of undisturbed stre 


critical speed at flutter ) 
l (Wa 

“effective piston velocity” 

(OY/ot) + L(OY/ox 

value of w resulting from ap and 


value of w resulting from } 

coordinate parallel to airfoil chord | 

distance of axis of pitching oscill 
tions behind leading edge 

distance from axis of pitching oscil 
lations (or elastic axis) to center 
of mass 

oO airfoil 


coordinate perpendicular 
chord, measured away from airfoi 
on each side of it 
equation of airfoil surface (upper or 
lower ) at time ¢ 
; . { lower . 
equation of ( ) surface of sym 
upper 
metrical oscillating airfoil 
normal force (positive downwar 
per unit span on airfoil 
subscript for undisturbed stream 
initial angle of attack of airfoil 
variable angle of attack (in pitching 
oscillations 
amplitude of airfoil pitching oscilla 
tions 
(M2 — 1)? 
ratio of specific heats (assumed co1 
stant = 1.4 
larger of the two dimensionless 
quantities 7 and e 
maximum displacement in oscillation 
divided by airfoil chord 
w,/@s = @/@ 
ratio of wedge length to airfoil chord 
(see Fig. 3 
spanwise coordinate divided by wing 
semispan 


pl 


wing density ratio = m 

X/C 

Xo/¢ 

Xa/ 

density of air 

ratio of density of air at altitude t 
that at sea level = pp,/po 

density ratio at 35,330 ft. altitude 

thickness-chord ratio of airfoil 

bending and torsion mode shapes, 
respectively, of wing 

frequency of wing oscillation 

frequency of oscillation at flutter 
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LINEARIZED PISTON THEORY 








Lyti/ w'M L.204 
0.70 0.75 080 0.90 0.95 
~ ~ Tt 1 @) 
Possio__4 AR 
THEORY + L a) 
+ -O0.10 
; 
-O15 
a+ eo 1 -O.15 
dC, /da ~ 
’ ie —=k7B(Ly+il, 
: dcz/da),-0 
Fic. | Vectorial representation of two-dimensional dynamic 
lift curve slope with reference to steady-state value 4/8(7 = 0, 
a () 5 
2 = natural frequency of wing in flexure 
= (1/27) (K,/m)!/? 
Wa/ 2 = natural frequency of wing in torsion 
: 1/9 
= (1/27) (Ke/Ie)*/* 
uw _ w Wa 
Ww = Wh/Wa 
r = 1;?/Iy1: 


= differentiation with respect to x 
= differentiation with respect to tim 


(1) INTRODUCTION 


OLLOWING the original work of Hayes,' Golds- 
F worthy,” and Lighthill,* the piston-theory con- 
cept, particularly as applied to problems in aeroelas- 
ticity, has been developed further by Ashley, Heller, 


Landahl, Mollo-Christensen, and Zartarian’~ at the 
Massachusetts Institute of Technology. Lighthill 


made use of a truncated binomial expansion of the 
“simple-wave”’ pressure relationship, Eq. (1), to derive 
aerodynamic coeflicients for oscillating airfoils at high 
Mach Number. 


(1) 


An outstanding advantage of this approach is the simple 
manner in which the effect of airfoil thickness is taken 
For a given profile at zero angle of at- 
It can 


into account. 
tack this effect depends on the parameter l/r. 
be shown that first-order piston theory (linear approxi- 
mation for p/p.) gives results for the load distribution 
which are identical to those obtained by assuming 
1/M?) < 1 in 
supersonic theory for a wing of infinite span. 
the 
gives terms in the load distribution which are linear in 


the solution derived from linearized 
Inclu- 
sion of second-order terms in pressure relation 
thickness, and the third approximation gives terms 
quadratic in thickness. The magnitude of the param- 
eter .\/7 offers an indication of the number of terms 
needed in the expansion of Eq. (1) for determining forces 
on oscillating airfoils at zero angle of attack. 

It is interesting to note that the results of linearized 
piston theory are equivalent to the expressions given 
by Collar’ if 1/2 >> 1. Collar attempted to derive the 
flutter 
speeds from a rather intuitive approach by assuming 


two-dimensional derivatives for supersonic 


that Ackeret’s theory of steady flow might be applied, 
without modification, to the case of nonsteady flow. 
Temple and Jahn’ provided a mathematical basis for 
Collar’s work and showed that his results gave the cor 
rect (linearized) values for the limiting case of Mach 
Number approaching infinity, 1/* > | 

Ashley, Heller, and Zartarian® 
applications of the piston-theory to problems in aero 


have given some 
elasticity such as the flutter of a rectangular wing, the 
flutter of a flat thin panel, and unsteady transients due 
to gusts. They have also done some elementary trend 
studies for preliminary design. In reference 7, in the 
section on three-dimensional flutter problems, they re 
ported results obtained from strip theory and a stand 
ard Rayleigh-Ritz type of analysis making use of un 
coupled modes. Their work is extended here to include 
a mode-shape parameter, and extensive parametric 
studies are carried out for preliminary design work, to 
aid in the estimation of the effect on flutter speed of the 
airfoil thickness and profile, initial angle of attack, and 
changes in elastic and inertial properties of the wing. 
Adding initial angle of attack to Lighthill’s work, ex 
pressions containing elementary integrals are developed 


and .J\/ 


1) which are tabulated for five different 


for the nondimensional flutter coefficients L 
(2 ‘, 2. a 
profiles. 

PISTON THEORY 


2) LIMITATIONS OF THI 


The inaccuracy and inadequacy of linearized theory 
of oscillating airfoils in supersonic flow, and the com 
putational difficulties of existing nonlinear theories, 
led Lighthill to propose a simple nonlinear theory for 
high Mach Number, based on the piston concept. 
Since the publication of Lighthill's work, the Fluid 
Dynamics the 
Institute of Technology have made extensive studies 


Research Group at Massachusetts 
of the linear and nonlinear theories making use of this 
concept and have pointed to its usefulness for problems 
in aeroelasticity. It should be instructive to briefly 
review some of the salient features of these studies. 
Landahl, Mollo-Christensen, Ashley 
that linearized piston theory can be used for arbitrary 


and showed 
small motions of thin two-dimensional airfoils when 
ever any one of the following conditions is met 

kM? > 1, 


(b) 


(a) MM? 


> | 


») 
Ae | 


(c) (RV)PF>D>1 | 
In order to illustrate these requirements, a comparison 
is shown in Fig. 1, for rigid plunging and pitching os 
cillations, between the results of the linearized piston 
theory, the more exact Possio theory,"’ and the theory 
suggested by Collar. All these theories are applicable 
to a very thin airfoil of infinite aspect ratio, at a small 
angle of attack, and oscillating in a nonviscous, un- 
separated flow free from strong shocks. 

Fig. 1 shows a vectorial representation of the dynamic 
two-dimensional lift curve slope at any frequency in 
= 0) value of +8. This 
For the axis of pitch at 


terms of the steady state (k 
ratio equals k?8(L; + 714). 
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midchord, the limiting value of this ratio as k > 


y 
obtained from Possio theory, is §/.\/. Linearized + 
piston theory gives the same value for /* > 1. 
It is seen that for any value of the reduced frequency 
k the errors arising out of the simpler theories become = i 
progressively smaller with increasing Mach Number. Lp — 6,0 —enhuap ibe 
If 1/ is not large, the results of the linearized piston j @ N > x 
theory, namely, 1; 0, Le 1/(R.A\/), ete., are correct a, 7 4 r 
ge - P S— /*b,a(t) 

only for k ~ ©, as predicted by the requirements of ——e Frit * ; 

— , . . ‘ : ENTER OF MASS 
Eq. (2). Collar’s results for 1), 2, ete., are identical dee eames 
to those of piston theory except that .1/ is everywhere , 
replaced by 8; at moderate values of .1/ his expressions z 
are in error for all values of k. If, as indicated in refer- Fic. 2. Coordinate system and notations for an oscillatiy 

. : oa airfoil section 
ence 9, his formulas for 1,’ and .\/,’ are multiplied by a cites 
[1 — (1/g?)], they agree with the initial terms of the 
: parameter on flutter speed and frequency of oscillating 


frequency expansions given by Possio, while the validity aes : : 
: : I 2 Seer é eR, airfoils at high Mach Number. 
for large .\J and moderate k is not affected. Thus 
Collar's corrected results are valid for k — 0, and " ‘ 
. ; am ae! (3) AERODYNAMIC FORCES ON OSCILLATING 
linearized piston theory is valid fork —~ ©. For inter- iden 

° ° . . é “OILS 
mediate values of k, both theories are in error except 
for large values of /, say \/* = 10. Consider a symmetrical airfoil at an angle of attack 
The authors of reference + also worked out a tip a (see Fig. 2). Let the equation of both surfaces 
(with y measured away from the airfoil in each case 


correction for a rectangular wing, disregarding the span- 
(x). Then the steady-state downwash is 


wise bending. The correction is based on an expansion be y ) 
in powers of k, given in reference 11, and can be utilized a ae: + Vx] 
with advantage within the limitations mentioned above. a s\9 Lr a x) | 3 

Landahl extended the above work by deriving cor- 
rections to the linear and nonlinear terms of Lighthill's 


The negative sign applies to the upper surface and the 
positive to the lower surface. Let the changes in y du 


piston theory. The corrections become important at to oscillations be },(x, 7) and the corresponding down 
lower Mach Numbers, say in the neighborhood of .1/ = wash be w;(v, /). Then the total downwash is 
2. The effect of 6/.1/° and 6/.\/° terms in the pressure cy! = la, 4 , 

ane ; : , ; Ww x) F Ueg F wis, 1) 
coefficient is shown in reference 5. Landahl’s solution 
is restricted to frequencies of oscillation for which where w(x, t) = (OY,/dr) + U(OY;/dx) 5 
(k/\J)* < 1 and to thickness or amplitude of oscilla- . ee : = : . — 
i aoe Substituting Eq. (4) in the cubic expansion of Eq. (1 
tions such that (6.1/)* < 1. He showed that the | 

namely, 


second-order terms in Lighthill's calculations are larger 


than the error in the linear term only when 6.1/* > 1. pp 1+ y(wa_) + [y(y + 1) 4] X 
In the sections that follow, simplicity of Lighthill's (w/a)? + [y(yvy + 1)/12] (wa 
approach is retained in order to make a detailed pee ae 
; ace . ; and subtracting the pressure distribution on the upper 
parametric variation and to study the effect of each 7 ‘ : / 
: surface from that on the lower surface, the load dis 


tribution on the airfoil is 


i 


¥+ 1 , y+ 1 aS + | | 
Ap p 24 || ave + - (Vay) WY,'(x) + (Vay) M°Y,’2(x) + 12 arend” | “— 


y+] yt | y+1 
[ = VY (x) + M?2Y,'2(x) + — , Max)? (w/a - 


yar] : ie ae 2 
k (la | (w,/a.)? + - (w,/a_) i 
{ 12 J 


Eq. (7) is seen to have one part which is independent of the downwash due to oscillations. This is the steady-state 
loading on the airfoil due to angle of attack and profile and may be extracted to obtain the corresponding lift and 
5 a I 'd > 

pitching moment of the airfoil. Thus the two-dimensional lift is 


y¥+1 y+! a y+1 
L = Apdx = 2yp Mave | + - (May)? + —— M | V,/(e)dé + - MW } (Ode | s 
J0 vd 0 0 , 


12 ‘ 


Defining OC, 0a = (OL, 0a), (1/2) plc, the steady-state lift coefficient slope is 


OL l 4 l o 4 l 
( ‘) = (4 uw) oe us (Mav)? + - 3 M } Y,’(t)dt + - 7 we Y,’2 (£) dé 


J 


(Y 





J0 | . a 


Oa 9 








SIT 


alt 


lor 


pre 


illating 


attack 
urfaces 
1 Case 
is 

5) 
nd the 


y due 
down- 


Ipper 
| dis 








AEROELASTIC INSTABILITY AT HIGH MACH NUMBER 249 
For an airfoil with sharp leading and trailing edges, }’,(0) }.(1) = 0, so that 
oc, y+] v++1 1 | 
(4 AM) 1+ (lla ~ V/ } t) ds 10) 
Oa s.s8 | } 7 UO J 
Similarly, the pitching moment about an axis distant £,c from the leading edge is 


: i oe ; , { 
] é F fc \di 1) 
( , ) M | Vy, *(&) (& — é 7 


(OM Oa) (1/2) pl °c’, the steady-state moment coefficient slope is 


rn IC ) A (? ; ) uf 
Ci 


Defining (OC, Oa) 


(*r") 4M) 4142 
Oa/ . i| 


Y."(t) (£ — &)dt 4 
l 


)a |. V./*(&) (¢ : dey 12) 


As to the time-dependent loading on the airfoil, it is 


now assumed that the downwash due to oscillations of es 
ince, vA -2pcU*k? } (Li + tLe) X 
the airfoil is small compared to the steady-state down j 2)) ] ] ) : 
Bai . m ‘ t/(c/Z) t 3 T t£4) @ (1) 
wash due to initial angle of attack and profile— that ts, 
w, a.) is small compared to Ja) and J/Y,’(x), so mM — pco?U2k?} (My + iM.) X 
that the last two terms of Eq. (7) may be neglected. [h/(c/2)] + (My + iMya} 16) 


For an airfoil oscillating vertically and pitching about 


the axis ¥ x toc. y 0 For a harmonic oscillation of the airfoil with / he" 
and a ae’ h lwh and & lwa can be substituted 
w,/a_) Mal(r) (ha) + (@a_) (x Vy) and the expressions for the sectional flutter coefficients 


; : . Serre become 
hence the oscillatory loading on the airfoil is 


| Ly (OZ Oh) (1 4pl*k?) 
Ap p_) Py) l + [Cy + 1)/4] (Way)? 4+ L. (OZ Oh) (w tpl 2k?) 
(y + 1)/2JMY,"(x) + [4 +1) 4] xX i —(OZ Oa) (1 2pcl-k*) 
WeY.*(x) i [i a.) + Ma + (a/a,) (vx — x0)| (13) e ~(OZ Oa) (w 2pcl *k?) 17) 
: : | — | ( 2 pl ) 
and the total load and pitching moment are Mi OM oh) | 2pcl*k 
Me —(O.\/ Oh) (w 2pcl 7k?) 
Zz Ap dx M; = —(O.M Oa) (1 /pc?l*k*) 
“0 M, —(dO.1/ Od) (w pe 2[ k 
and 14) 
. Making use of Eqs. (13), (14), and (17), it can easily be 
ies Ap(x — xo)dx shown that 
Following the lead of the National Advisory Committee Ly M, = 0, Mz Ls t 18) 
L (1/k)Leo, M (1 k)M 1/k)Ls! 


for Aeronautics, it is standard American practice to ex- 


press the section force and moment as 


so that the flutter coefficients that need be evaluated 
are Lo, Ly, and .1/;, which may be expressed as 


1 gf y¥+1 |, yt! i ¥+1 
L aee (Ma)?| +2 Vu Vi(e)dt +2 NE v.28) ae! 19) 
Vk tli | 2 Jo } Jo / 
ly y+ 1 os 
l i+ (lla a — 2) + (vy + IM V.’(&) (&€ — &) dé 4 
Wk | J0 


—— = l | 
: we | V./*(&) (é t,) de! (0) 








Y,(€)= 240 €-€9 


of supersonic 


airfoils which are shown in Fig. 3. 


Matriy 


6 
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‘ Writing LZ, = (Wk)L», e (Mk)Ls, and VW, = (\/b)) 
i tere the flutter coefficients may be expressed in a 
= ees ereregs 8 esate, form: 
= x=cé =tc/2 7, <€<(I-,) 
HC ie. )\- 7) ——7,c = i te im, )<E<! [. A, iz r 
- / r / POSITIV YF t x-AXIS - 
y |S POSITIVE AWAY FROM x a As 7 ] Q 
= c( 7) - , 
i - al V A / R 
GENERAL DOUBLE here { | 
waieeeeen where A, 
A» 1 — 2 
A; 1[(1/3) £o + £07] 
y - ‘i P 
i rc 7 (2/7) V.’(é)dé 
! 7 UV 
- * F l 
Tl: (1 ry | V./2(&)dé 
DOUBLE WEDGE ‘ a 
22,205 7" 
T3 (4/7) V./(é) (& — &) dé 
y J VU 
i a 
! T, (2/7?) V./2(£) (& — &) dé 
<= Jo 
= 7C (I-Ie at’ >! 
Los ; on T. (S/r) | V.’(é) (& — &)* dé 
J0 
SENERA INGLE 21 
WEDGE Ye(€)e@c/2mlE OK EK T° (4/r*) | Y,'*(£) (& — &)* dé 
y =rc/2 N<E <I J 0 
' ! 
a = i and P=1+ [(7 + 1)/4] (a)? 
J QO l(y + 1)/4] (7) 
| aaa EDGE R a [(4 4. 1) 1] (Mr)2 
, be Note that A», and A; depend upon location of pitch 
‘ axis, 7; and 7» on airfoil profile, and 7(7 = 3, 4, 5, 
le = os - upon airfoil profile as well as location of pitch axis. 
The coefficients 7°;(7 = 1, 2,..., 6) are tabulated in 
PARABOL IC Bl INVEX< ae 2 
lable 1 for two general and three particular profiles 




















Fic. 3. Airfoil profiles and their equations 
TABLE I 
Flutter Coefficients for Five Profiles 
General Parabolic 
Coefficient General Double Wedge Double Wedge Single Wedge Single Wedge Biconvex 
T, 0 (¢) 1 1 0 
T (nm + m)/(un nD) 1 1/(4n) 2 = 
2 mh * Mets qT 3 
Ss 
T3 - (2 ” (m + np)| me yo 2h, 1- 28, i 3 
4 
1, [my - (my + BE] /(2nyrp) 1 = 2, (n= 28 )/(un) | (1 = 28) | 4 (2 - 26,) 
1 rd 2 a 8 \ 
Ty -u[1 - », - 3(n° - 1) -2(1 = 2&,) LB w - %, +0) uf -&+ ‘- 3(1 - 2€,) 
-[2 ~/}* no) co} 
1 1 £8 1 2| 16 [2 2] 
1 [ry(2 - 1) + 3 mM (ny + 9) uff - 2, +e,4Bn- e428, a i OMe FB 6 & 





= 2nE + (m, + ME“ |/(mm) 
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LEROELASTIC INSTABILITY 
lTb\ . . 
TR) M; 1) SOLUTION OF THE FLUTTER DETERMINANT 
Matrix 
The critical speed and the corresponding frequency 
vibration for aeroelastic instability of an airfoil 
») 
u(e? —1) + (114+ 2 
ZiSa Tt AM, + VW) Mr a 
Making use of equalities (1S) and separating the real 
id imaginary parts yields the following solutions for 
the critical speed and the corresponding frequency of 
vibration 
( aQy1T l — w,) (1 "ie ‘ 
’ | 2) 25) 
2, TT (3 Wy) J 
E 1 + 14a)" 
Wy = dg 4)) 26) 
1+ ms E14 
where 
i pM /L ; Ly Ls, ms Va 1 
Vex M4 Ms Vo ayy 2 Pag 
fil 2§ 1) Ms + ls "= 
27) 
™ iil sa 4) L + a — Seka 
-v - 
1 + fish G Wy Was 
L Was 4 Ww Wy w Ww 
pitch 
5, 6 . , 
f lu(< — 1) + (Ly i iL») | J; (—=7 
aXIs. I ),.¢ ° { . 2] 
: I~ ae T (\J, + 1Ms) | Mla 
ted in 
rofiles | v1 
where iF o17(A)dA 
7 
>| 
I. o»*(A)dA 29) 
JJ 
| a 
[ I 1(A)@o(A)dAr 
J 
ie Separating the real and the imaginary parts as before 
—._—«{|_~=_ and introducing a modal shape parameter I [3° / Tye 
the following solution is obtained for the critical speed: 
U, ; a T (l—o*) U1 - ‘tel 4 . 
(1/2)z (30) 
Wea Qo + (3/7) 
where the corresponding frequency of vibration is ob- 
) tained by solving the following quadratic in @,° 
) lols -T as5(1 ; az) |e, T 
} laslag — 1 + T) — aoag — a; Jor + 
a;,(l1 — I’ — dg) 0 (31) 
The constants a, to a; are defined as 
g 2| 
° | a, = 47E,20/ 
# ” 
ds l(2mt, + 1.7) — fila — ™ 
(3 fil sa," 
(4 1+ N14 a fe OT (32) 
ads al,() — T) 
Ag 1 + 14a," 
} (ty a (é*/9.*) 


For Tr l, ay = 0, ay Q131;, G2 = Aa, G4, = Ca, and Eq. 


31) yields two values for a7, namely, af = ds dy and 
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section in two-dimensional flow and having freedom 
in vertical translation and pitch may be obtained by 
solving the following determinant |see reference 12, 


Eq. (27) ]: 


5) AEROELASTIC INSTABILITY OF A UNIFORM 


CANTILEVER WING 


Analysis of the Rayleigh type for a uniform canti 
lever wing involves the selection of a set of modal func 
tions to approximate the mode of vibration at the criti 
cal speed. The uncoupled first bending and first tor 
sion mode shapes are selected for this analysis and are 
represented as ;(A) and ¢.(A) where A is the spanwise 
coordinate divided by the semispan of the wing. Using 
the two-dimensional air forces and moments, that 1s, 
neglecting the tip effect which is small for high Mach 
Number, the stability determinant of Eq. (24) is modi 


fied to the following form: 


Qué. + (L3 + 7L,) |: 
iE os 1 RE STE) _ 
—(; dx. The first is the same as Eq. (26), while the 
second value is seen to be trivial as it leads to infinite 
critical speed by substitution in Eq. (25). Thus, the 
case [ 1 leads to the sectional solution of Eqs. (25) 


and (26) which may also be interpreted as ¢;(A) 
go(X) B 


(6) PARAMETRIC STUDIES 


It is now evident that the solutions for the critical 
speed and the corresponding frequency of vibration 
given in Sections (4) and (5) lend themselves with ease 
to parametric studies for use in preliminary design. 

First, the simple case of an airfoil in two-dimensional 
flow and free to pitch only is considered. The aero- 
dynamic damping moment is given by the out of phase 
component from Eq. (16), namely, 


M4 —pc*l kh? Ms — pc-a RAM)AMs (3:3) 


Substituting for 17, from Eq. (22), 


kM) |A;3P + 7;04+ 7T¢6R] (24) 


WM, — pc-a 


Considering the particular case of a double wedge pro 
file and substituting the corresponding values of 7% 
Table 1, 


and 7, from 


Sls - — pea .."(kM) (4{(1 3) — & +¢é 2] ; 1 + 
[iy + 1) 4] (Mao)? + [Cy +.1) 4] (M7)*} - 
l(y + 1) 2) (1 - 2¢)) (Mr)) 


(35) 
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Fic. 4. Aerodynamic damping in pitch (double wedge profile 


Fig. + shows the dependence of the aerodynamic damp- 
ing in pitch on the location of the axis of pitch (&), the 
thickness parameter .J/7, and the initial angle of attack 


parameter .J/a). For the zero thickness airfoil at zero 
angle of attack, the damping is maximum for & = 0 
and | and minimum for & = 1/2. Angle of attack in- 


| 


creases it for all values of & by the factor [1 + (y + 
1) (ay)? 4]. The effect of thickness is to decrease the 
damping for locations of the pitch axis forward of the 
mid-chord and to increase it for 0.5 < & < 1.0. The 
curve for .\/7 and / a) each equal to 0.25 shows a blend- 
ing of the two effects. The outstanding feature is that 
the damping is always positive, thus eliminating the 
possibility of instability in pitch at high Mach Number. 

It is seen from Eqs. (25), (26), and (27) that the 
following dimensionless parameters aifect the critical 
speed and the corresponding frequency of vibration 
(independent of altitude) of the airfoil. In these, the 
frequencies have been normalized with respect to w, 














Fic. 5. Variation of flutter speed and frequency with 
moment of inertia in pitch. (Double wedge profile, constant 


wing stiffness in torsion; 7 = ay = @& = 0, & = 0.5, ta = 0.1 


and the lengths, except the radius of gyration in pitch 
with respect to the chord length c. 
(1) r,°, or the moment of inertia of the airfoil about 
the axis of pitch 
(2) &, chordwise location of the axis of pitch 
(3) &, mass unbalance or distance between axis of 
pitch and axis through center of mass 
(4) a, (w, Wy), ratio of uncoupled flexural and 
torsional frequencies of the airfoil 
(5) wl/, altitude parameter 
(6) 7.1/, thickness parameter 
(7) a7, angle of attack parameter 
(S) I’, modal shape parameter 
(9) y, ratio of specific heats (assumed constant 
1.4) 
The edects of the first five parameters may be in 
vestigated for an airfoil of zero thickness at zero angle 
of attack. For this case, P ee, R 0, so that 


Le A, . i Ly = A» =|]— 2Ey, and M; A 


$[(1/3) — £ + &?]. Substitution in Eqs. (25), (26 
and (27) gives 
U; j ft 2a? + ry2(1 — a7) (1 — €?) prs a 
= (1/2)uM - . ee — : if} 
Cie iuAf int, — (1 — 96,) (1 = £9) + (1 — 26)? — 4401/2) — & + 6,7)9 
; 1 + 4[(1/3) — & + £02] (G,2/r2”) ” 
Wy ss > » Pa - 9 9 uy 
1 — [4€.(1 — 2€0)/r.2] + 4[(1/3) — & + &?]/r,’ 
The parameters will now be varied individually. 
Moment of Inertia in Pitch ‘ 
U,/ (cw) = (1/2)uM r,”? X 
Let w, = O and & = 0.5; then 4(t,2/r.2) + (1/3)/re? \' 2 ss 
(Oe 


@? = 7r,?/[ra? + (1/3)] (38) 
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€, , MASS UNBALANCE IN CHORDS 


6. Variation of flutter speed with mass unbalance and 
chordwise location of axis of pitch 


Double wedge profile; 7 


smce 
Wa (x. Ly , 
(K,m)'*/(1/2)e rq 
Eq. (89) becomes 
x m)' 
uw? wl lo 10) 
( 

EU + (1/3)] [2uME, — (1/3)]8 
ind 

wy (2c) (K,/m)" *]/[r.2 + (1/3) ]'* (41) 


For prescribed values of A,, m, &, and wl, ly, and 
are seen to be proportional to [7,? + (1 3)] ie? 
Normalizing Eqs. (46) and (41) with respect to their 
values for 7,° |, the variation of flutter speed and 
frequency with radius of gyration and its square is 
cannot be less than 


shown in Fig. 5. Obviously 7, 


¢,; however, it is interesting to note that the critical 
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€,, MASS UNBALANCE IN CHORDS 


FiG. 7 Variation of flutter frequency with mass unbalance and 
chordwise location of axis of pitch 
ah = 0, re? = 0.25, pM = 250 


speed and frequency can be doubled as the mass dis 
tribution of the concentrated 
masses at the leading and trailing edges inward to the 


airfoil is varied from 
center of mass, holding the torsional rigidity constant. 
Similar design charts may be prepared for other values 
of @, and &». 

In the past, the aeroelastician has been concerned 
more with the critical (flutter) speed than with the fre- 
With 


the ever increasing use of servo and guidance systems 


quency of vibration of the airfoil at this speed. 


in piloted and pilotless aircraft a knowledge of the fre 
quency becomes imperative in order to avoid servo 


elastic instability. 


Location of Pitch Axis and Mass Unbalance 

It is evident from Eqs. (36) and (37) that the effects 
of the chordwise location of the axis of pitch and its 
distance from the center of mass (the mass unbalance) 
on the critical speed and frequency are interconnected 














° 50 100 150 200 250 
ALTITUDE PARAMETER, uM 
With altitude at fixed frequency ratios 


r= aw = 0,8 = 0.5, fa = 0.1, ra? = 1/3 
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with the frequency ratio a, (w, w,). However, it is 


well understood that, for a given mass unbalance, 
forward shift of the axis of pitch increases the flutter 
speed at all values of a. For low values of a, (0 S 
a, < 0.5), the flutter speed generally decreases with 
increase in &,, for practical values of the mass unbalance 


(O < & < 0.25). For moderate values of @,(0.50 


Sa 
a, < 0.70), the flutter speed has a minimum approxi- 
0.10; further increase in mass unbalance 
1.0, the mass unbalance 


mately at &, 
becoming beneficial. For @, 
is beneficial for all positive values of &,. 

Some typical curves for 0 are shown in Fig. 6. 
It is instructive to examine the combined and the indi- 
vidual effects of & and & on the critical speed and 
frequency. 

The full lines show the variation of Ul’, (cw,) with 
£, for a fixed location of the axis of pitch, and the dotted 
lines show the effect of moving the axis of pitch for a 
fixed location of the center of mass. A vertical inter- 
cept gives the effect of moving both axes forward keep- 
ing &, constant. For example, start with point A with 
t = 0.5, é, 0.15 and U, 6.16. De- 


creasing the mass unbalance by 


U, (CW) 
10 per cent (&, 
0.05) raises (", to 9.65, giving an increase of 57 per cent. 
Moving the pitch axis forward by 10 per cent (& 
0.4), leaving the center of mass at 0.65 ¢ raises [", to 
Moving both 


0.15) raises 


7.23 giving an increase of 17 per cent. 
axes forward by 10 per cent (& 0.4, & 
(’, to 10.38, giving an increase of 6S per cent. The 
advantage of decreasing the mass unbalance is thus 
apparent, especially for midchord and forward locations 
of the axis of pitch. 
Fig. 7 shows the corresponding variation of ay 

Ww, w, With & and &. For small values of &, the change 
in w, with & is small, but it becomes significant for 
larger values of &. In general, a, increases as & 
and (& + &,) are decreased. 
0), a, decreases as the axis of pitch is moved 


But for a mass-balanced 
airfoil e. 
away from the midchord location, the minimum value 
being for & = O or 1 and equal to 0.3975 in this case of 


o, = Oands,? = 1/4. 


Torsional Frequency and Its Ratio to Flexural Frequency 


Eqs. (36) and (37) show that the critical speed l 
and the frequency w, increase linearly with torsional fre- 
quency w, if the flexural frequency is also raised pro- 
portionately leaving @, unchanged. However, if «, 
is held constant, the effect on Ul’, of raising w, is differ- 
ent for w, > w, from that when w, < w,. Some typical 
curves of U, versus a, for different values of u.l/ are 
For w.l/ = 250, take the initial value 
Doubling the tor- 
3.) 


shown in Fig. Sa. 
of @, 0.8 for which U, = 3.8. 
sional frequency raises l’, by a factor of 2(6.5 
3.32. On the other hand, if the initial value of a, were 
1.6, doubling w, leaves l’, unchanged. Doubling the 
flexural frequency at w, = 0.S gives U, 7.65, which, 
in this case, doubles Ll’, Thus an obvious conclusion 
that may be drawn from these curves is that, for 0 < 


a, < 1.0, the wing should be stiffened torsionally and, 
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ABLE 2 
Effect on Flutter Frequency of Doubling Airfoil Tors 
Flexural Frequencies at Different Frequency Rati 

] 1 





Wo/ We 

, a 
a. wo doubled @ doubled 

h h awe 

a eae en 
0.25 1.96 1.09 
0.50 1.84 1.27 
1.00 1.58 1.58 
1.50 1.38 eto 
2-00 1.27 1.04 





for a, > 1.0, it should be stiffened flexurally 1m order 
to raise its critical flutter speed. 

Fig. Sb shows a cross-plot of Fig. Sa. <All curves 
are seen to be of a similar nature irrespective of the 
ratio. Another common char 


value of frequency 


acteristic is the steep rise in U, for values of «lJ less 
than 20 
than S or so. 

As for a,, it is independent of w.J/ and its variation 
The effect on w, of doubling 


that is, for values of wing density ratio less 


with a, is shown in Fig. 9. 
®, OF a, at different values of a, is shown in Table 2 

It may be well to point out that it is the square 
appears in the expres 
sions for U, and ay. ¥ A.) the 
parameter (A, A,) could be used instead of (@, w,) for 
A replot of Fig. Sa (for u.1/ 250 only 
10 where the abscissa is 


the frequency ratio (q,") that 


Since @),” r.c/2)*(K, 


these plots. 
and of Fig. 9 is shown in Fig. 


@,” instead of a. 


Altitude Parameter uM 


The variation of Ul’, (A, m)°~ and U’, (cw,) with ul/ 
for representative values of 7,° is shown in Figs. 
and I1b, respectively. These curves as well as those 


of Fig. Sb indicate that l’, varies approximately as 


l r1.* . . 
(uM) This observation leads to rather simple 
results. 

U,,/l [(ulW)>o (uA), ] > - 1? 
Squaring both sides and substituting .1/ U/a, and 
mm ky o gives 
l f py: 10}) dog 
and 
M, M, (701) (o-oo) Dp, p 


Given the flutter speed Uy, and the corresponding 
Mach Number .\/,, at one altitude, their values at an 
other altitude may be obtained from the above relation 
ships which are plotted in Fig. 12. Furthermore 
for 35,330 < h < 105,000 ft. where (a2 a;) 1.0, 
and ./, become inversely proportional to the densit) 
Fig. 12b shows a plot of (¢, o) for 35,330 < 
o at 35,330 ft. altitude. 


ratio. 
h < 105,000 ft. with o, = 
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Fic. 12. Variation of flutter speed and Mach Number with wedge profile; 7 = ay = ah = 0, & = 0.5, fa = 0.1 Fic. lla 
etude Fic. 12a (top 0 <h < 35,330 ft. Fic. 12b (bottom top Constant wing stiffness in torsion. Fic. 1lb (bottom 
Constant wing torsional frequency 


0,33 } 105,000 ft 








) 


256 JOURNAL OF THE AERONAUT 





[CAL SCIENCES 1958 


APRIL, 


a. Double wedge profile, & = 0.4 b. Double and single wedge profiles, & = 0.5 ‘ Double wedge profile ¢ it 
Fic. 13 Effect of thickness on flutter speed fa = 0.1, ra? = 1/3, uM = 259 
Alternate Approach where 
The flutter frequency is independent of altitude as f, 24» 11 — 3) (1 2) /g 
seen from Eq. (26) or (37). Thus for given inertial f, Fag _ (1 ~ 9¢) (1 — ¢2) ( 1 
and elastic properties of the airfoil and location of the f. _ 9¢)2 ee i of ° 
; —- 2&)? — 4 [(1/3) £o + £7] -1 3) 


axis of pitch, Eq. (36) may be written as 


Squaring both sides of Eq. (45), substituting / 


U, (1/2) (cw)@, uM hi ((uAL)fo + f ty 2 (45) : i ; 
aM and solving for \/ yields 
M, (1 A)cPwe?@P7 (fi fo) (ua) — (fs fo) (1 /m) 
[(1/2)ew,|?(fi fo) (uo oa?) — (fs fo) (o/ uo) 
((1/2) (€ w, ay) 7h fr) | wo (p po) | + (1/3f2)(o/po) {7 


(RAM )o7(fi fe) [wo (DP) Po) ] + C1 38f2) (0 uo) 
[uo fo(P) Po) | [fil(R.D)o? + (13) (a ao)? (o/ mo)? 


The second term is much smaller than the first except 
when (&.)/)) and wo are both small and of the same order 
Hence the flutter Mach Number 
versely proportional to ihe pressure ratio, the same rela- 


of magnitude. 's In- 
tionship as derived in Eq. (44). 
Since (~P/ po) 


pressed as 


(oa*) ao”, the flutter speed may be ex- 


el (1/2) poo [ao7 U0" 


9 
2) [po(aomo)?/ (ppo)] [fife (RMD? 2 
l(p/ po) (49) 


that is, the dynamic pressure at flutter speed is inversely 


or Gr = 


proportional to the pressure ratio, or, the product of the 
square of the indicated speed at flutter and the pressure 
ratio remains constant. 

It is evident from Eqs. (47), (48), and (49) that, be- 
tween the altitudes of 35,330 and 105,000 ft., where the 
speed of sound in air is constant, the flutter speed, 
Mach Number, and dynamic pressure vary inversely 
as the density ratio. 


Airfoil Thickness and Profile 


In order to study the effects of airfoil thickness and 
profile, the airfoil is assumed to be at zero angle of 


WO, = 


[ P al, [domo ala‘ ay )f> | y 4 


[fi(RMD) 5? + (1/3) (a/ao)? (o Ho)" | tS 


By virtue of the aforesaid, the flutter speed is invers 
proportional to o(a/ao), the same relationship as shown 
in Eq. (43). 

The dynamic pressure corresponding to the flutter 


S] eed iS 


o*(a, do)*fo"| [fi(RAD)o? + (1/3) (a/ao)? (wo)? |? 


attack. 


and the general single wedge. 


Two profiles are considered; the double wedge 
The former is repre 
sentative of the family of closed profiles for which 
V’.(0) VY.) = 
the family of profiles with blunt trailing edge, that is 
V1) ¥ 0. 

13a, 


0, while the latter is representative ol 


13b, and 13c show the effect of airfoil 


thickness on the reduced flutter speed U’, of a double 


Figs. 


wedge profile for three different locations of the axis 0! 


pitch, namely, £& 0.4, 0.5, and 0.6 respectively, over 


a range of (a,/w,)? from 0 to 4. U, has been normal- 


ized with respect to its value, 7.67, for & = 0.5 and 7 = 


0, and is denoted by U’,’. It is to be noted 
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that for 0 < a, < 1.0 airfoil thickness lowers the flutter 


speed, the severity of the effect decreasing as the axis 
f{ pitch 1s moved aft. However, l’, with &, 0.4 
nd A7 0.25 is higher than the l’, with € 0.5 
ind Wr 0: and the U’, with &, 0.5 and Jr 0.25 
is higher than that with & 0.60 and Jr 0. It is 


easv to Show that for steady-state conditions the center 
if pressure is located at a distance [& + (1 2) (Ly Le) |e 
from the leading edge; that is, at about 0.43c for the 
double wedge profile at zero angle of attack. It would 
seem then that the location of pitch axis at & 0.4 
should yield the double advantage of high flutter speed 
ind low hinge moments for use as a control surface. 

On the other hand, for a, > 1.0, the airfoil thickness 
raises the flutter speed and becomes increasingly bene- 
ficial with rearward movement of the axis of pitch. 
lhe dependence of thickness effects on frequency ratio 
was also shown by Garrick in the 1957 Minta Martin 
Lecture. 

The corresponding curves of airfoil frequency at 
flutter speed are shown in Figs. I4a, I4b, and I4e. 
Again, the flutter frequency ratio is normalized with 
@) 0, 


resp-ct to its value, 0.707, for & = 0.5 and 7 


ind is denoted by @,’.. Thickness decreases the fre 
quency and the effect decreases with aft position of 
the axis of pitch. However, since mass unbalance is 
kept constant in these curves, the frequency increases 
with forward shift of the axis of pitch, as already shown 
in Fig. 7 for @ (0). 

The effect of changing the profile shape is shown in 
For a general single wedge (n 0.5, Fig. 
se) and 7 0 the curve is the same as for the double 
wedge. But for 7.1J > O the flutter characteristics of 
this profile are different from those of the double wedge 
as shown. It is to be noted that, for the single wedge, 
the comparative thickness effects are smaller for 


Fig. 1b. 


() « . 0.5; they are slightly worse for 0.5 < a,” < 


!.), and are substantially lower for @,? > 1.5, that is, 
less beneficial to flutter. 


Decreasing n, that is, decreasing the length of the 
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a & = ().4 


INSTABILITY 
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Effect of thickness on flutter frequency at zero angle of attack 
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wedge, reduces the flutter speed for 0 < a, 1.0 and 
increases it for a,” > 1.0. The effect, however, is rather 


small. 


Initial Angle cf Attack 

For an airfoil of zero thickness, the initial angle of 
attack reduces the flutter speed by a constant factor. 
With the parameters of Fig. 13b and la 0.25, this 
multiplier is 0.982. For 7.1/ 0.25, the effect of the 
initial angle of attack is almost negligible for 0 < a, 
1.0; for a, > 1.0 it is shown in the Figure and is of the 


order of 2 to 3 per cent reduction in flutter speed. 


Modal Shape 

For a double wedge airfoil at zero angle of attack, the 
effect of the modal shape parameter I is shown in Fig. 
15 for two values of 7.\/, namely, zero and 0.25. In 
clusion of mode shapes in the analysis is seen to increase 
the flutter speed. It is interesting to note that the ini 
tial angle of attack and the mode shape have opposite 
effects on flutter speed and, over a wide range of the 
frequency ratio @,, they are of the same order of magni 
tude. Hence it seems that, at least for small angles 
of attack, these two parameters may be disregarded in 
the preliminary calculations of flutter speed without 


much sacrifice of accuracy. 


CONCLUSIONS 


The results of this study may be summed up as fol 
lows. 

Parametric studies of the type shown in Section (6) 
can be of invaluable help in preliminary design. 

Flutter speed can be raised considerably by decreas- 
ing the moment of inertia in pitch while holding the 
torsional stiffness constant. For average values of the 
radius of gyration (0.8 < rz, < 0.6) it varies almost 
linearly with r,. 

For midchord and forward locations of the axis of 


pitch, mass balancing can be used to great advantage. 





O05 Cc & = 0.6 


(Double wedge profile; a = 0.1, ra? = 1/8 
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= T ] speed over a wide range of frequency ratio. Neglectin 


——":025 them in preliminary calculations should not result ; 


a great loss of accuracy. 
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A Turbulent Analog of the Stewartson- 
Hlingworth Transformation 


FRED E. C. CULICK* ann JACQUES A. F. HILL** 
Massachusetts Institute of Lex nology 


SUMMARY 


The Stewartson-I]lingworth transformation is applied to the 
momentum equation for compressible boundary-layer 
y-coordinate transformation unspecified, how 


ving tiie 


It is shown that the transformed equation is the integral 





ientum equation for incompressible flow if (a) the effect of 
ressibilitvy on the boundary-laver shape parameter H/ can be 
nted bi 
esented by 
H+ 1 =(H,; + 1) (T/T, 
( the x-coordinate transformation is chosen to be suitably 
ited to the ratio of skin-friction coefficients in compressible 


ind incompressible flows 


Experimental evidence is presented which shows that condi 
a 


is satisfied for turbulent boundary layers up to .\/ 5 


nsformation is chosen according to (b) and an equation is 
1 


turbulent boundary-layer growth in 


resented which gives the 
mpressible flow in terms of a simple quadrature. The pre 
ictions of this equation are then compared with some measure 
ments on wind-tunnel nozzles 
SYMBOLS 
speed of sound 
local skin-friction coefficient 27,,/ p,1,* 
( average skin-friction coefficient 


D reference length 
boundary-layer shape parameter 6*/@ 
Mach Number 
\ exponent in empirical skin-friction equation (11 
x Reynolds Number based on distance from leading edge 
Reynolds Number based on momentum thickness 
temperature 
‘ velocity component parallel to wall 
coordinate along wall 
coordinate normal to wall 
coefficient in empirical skin-friction equation (11 
ratio of specific heats 
boundary-layer displacement thickness 
variable scale factor in the transformation, defined by 
Eqs. (3 
f boundary-laver momentum thickness 
rv viscosity 
kinematic viscosity p/p 
density 


T shear stress 


Subscript and Superscripts 


evaluated in the free stream just outside the boundary 


lavet 
corresponding quantities in the related incompressible 


flow 
evaluated at the wall 
Received June 24, 1957 Revised and received December 27, 
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= evaluated at an isentropic stagnation point in the free 


stream 


= evaluated at the initial point of the boundary-layer 


growth calculation 


t 


= evaluated at the reference temperature 7* 
DISCUSSION 


Sb PROBLEM of calculating the growth of the 
boundary layer on a body in an airstream, or on 
the internal walls of a duct through which air is flowing, 
has long been of interest to the aeronautical engineer 
Methods for handling this problem are continually be 
ing developed and improved. In recent years there 
has been much interest in methods applicable to com 
pressible flows. One approach, which has been very 
successful in the laminar case, is to seek a suitable trans 
formation which will make applicable also to com 
pressible flows the large body of results available for 
incompressible flows. It is the purpose of this paper 
to examine the relationship between turbulent boundary 
lavers in compressible and incompressible flows and to 
show that a transformation developed for laminar flows 
is approximately valid, with a slight modification, for 
turbulent flows as well. 

The specific problem which will be considered is that 
of the growth of the boundary layer in a two-dimen 
accelerating flow with arbitrary gra 
dient. The restriction to 
the difficult problem of predicting separation and yet is 


sional pressure 


accelerating flows avoids 
satisfied in many supersonic flows of practical interest. 
In incompressible flows of this kind the growth of the 
boundary layer may be calculated very conveniently 
by simple quadratures such as Thwaites’ formula! for 
the laminar and Truckenbrodt’s” for the turbulent case. 
If the flow is not too far from adiabatic, the laminar 
problem in compressible flow may be solved by apply 
ing the Stewartson-Illingworth transformation*: * to 
Thwaites’ formula, as suggested by Rott and Crabtree. 
The question which will now be examined is whether 
this same transformation is also applicable to the turbu 
lent problem, as suggested in references 6, 7, and 8. 
Using the conventional symbols, the integral momen 
tum equation for compressible boundary-layer flow 


may be written 


do dx) + O\(H + 2) (1 /u,) (du, dx) 4 


(i ‘p.) (dp, ‘dx) ] Ter! Poth Cy/2 (1) 


Stewartson originally transformed the boundary-layer 


differential equations using 
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TABLE 1 
Values of H; and of N and a(.V) in Eq. (11) Appropriate to 
Various Ranges Res for the Turbulent Boundary Layer in 
Incompressible Flow 


Reg \ at N 
LO0-5 . 000 } 0.0123 1.4 
500-50. O00 5 0 OO8S5 1.3 
3 000-600 ,000 6 0 0062 1.25 
10 ,Q00-10 az 0 OO48 1.2 
. 
/ } 
) 
\ n p p,)ay 7 Qo/a,ju, 1 / 2 
v7 


as well as transforms for x andv. Here 


q = T./T ott ?/20- [= (T./To)* for y = 1.4] (3) 


and the subscript refers to stagnation conditions in 


the external stream. The transformed integral equa- 


tion obtained by using Eq. (2) and leaving the x-trans 
formation undefined is 
(do’ dx’) + 6'|} 71+ 1) (7.7) + 1) X 

1/u,') (du," dx’) | (dx dx’)n(cy 2) (+) 
with A’ nO 5) 


This result may now be compared with the integral 
momentum equation for an incompressible flow; 


(dO; dx;) + 6; 7; + 2) (1) u,;) (du; dx;) Cy/2 (6) 

with fluid properties py and wy evaluated at 7). If the 
conditions 

(H + 1) (7./To) H,+ 1 ‘) 

Ix" Ix . f f S) 


are satisfied, the transformed flow may be considered 
incompressible and may be represented by the solution 
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Fic. 1. Variation of shape parameter, 7, with Mach Number 
for turbulent boundary layers 
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> 
v 
Fic. 2. Present theory for variation of Cr ¢ ind ¢ 
Mach Number (7* = 7 
of Eq. (6) with 
Q’ 0 i i 


The first condition is clearly one which must be checke 
by examining experimental results or otherwise; where 
the second may simply be used to define the require 
y-transformation. 

Both /7 and c; depend on Reynolds Number and pres 
sure gradient as well as on Mach Number. The pres 
sure-gradient effects in supersonic flow are not know 
well enough to be considered here; it seems that the 
are small in any case for accelerated flows. The Re) 
olds Number dependence, however, must be taken int 
account, and it is necessary first of all to note that t! 
values of momentum-thickness Reynolds Number 1 
the two flows represented by Eq. (1) and Eq. (6 


different; in fact, 


Reg) poll, O; polt,'0" uw 
(u. Mo) (p.u.6 u,) u, bo) Re 


The variation of c,, and /7; with (Re,); in incompres 


sible flows is well known. It is convenient in this cas 
and consistent with the approximations inherent 1 
use of integral methods) to use the approximate for 


mula 


Cr/Z a(N)/(Re,); T,/ Pott 


Very good agreement with the Karman-Schoenher 


formula over wide ranges of (Rey); may be obtained | 


choosing appropriate values of .V and a(.\V), as give! 


in Table 1. Also given are appropriate values of /7/, cot 


responding to each choice of .V. 


tar 


any compressible flow are known to fall in a certa 
range, it is a simple matter to calculate the correspot 
ing range of values of (Re,); from Eq. (10) and het 


to obtain suitable values of .V and /7;,. 


Now that //; has been defined, condition (7) may b 


checked experimentally for flows in which @ and 


have been measured. Fig. | shows the observed var 


ation of /7 with Mach Number for supersonic flows 0! 


a flat plate and on a set of wind-tunnel nozzles. Als 


shown is the variation predicted by Eq. (7) for the tw 





If the values of 6 1 


tw 
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ae... ilues IT, corresponding to the two sets of measure apparent because the latter have been chosen to repre 
nents. For flows with zero heat transfer at any rate, sent the effects of compressibility with Re,, rather than 
ndition (7) seems to be satisfied if lJ <5 Reg, tixed 
[he required x-transformation defined by Eq. (S The ratio c, c,, corresponding to the relation (10 
in onlv be evaluated in terms of the ratio of the skin between the Reynolds Number in the two flows may 
friction coeflicients in compressible and incompressible now be obtained by multiplying Eq. (15) by 
lows. This ratio may be conveniently evaluated if R R \ = 
the effects of compressibility are expressed in terms of ela) EO eal 7 
the concept of an eTective temperature /* as used, for and may be written 
instance, in reference 11. The formula for compressible ' 
( Mi My 1 i Me 
flow 1s written a 
tT. p a \ o*u,0 uw : : . 
= = : The x-transformation (S) for turbulent flow now be 
p, p Cy, 2 12 
comes 
with the same value of a(.V) as in Eq. (11 Various ie Pa = £1)/2 he 
( expressions for the ratio cy c,;, may be written on this ad daa ra ® ned 
asis, each corresponding to a particular relation be or 
tween the Reynolds Numbers of the two flows. For dest de roy os ; sa 
mparison with experiments, it is convenient to use aac ; ei si . 
= at YX (TOTO \ 13 rhe relation completes the definition of a transforma 
ee tion which (at least for small heat-transfer rates) relates 
hecke ind the value of 7*, which is purely empirical concept the integral properties of a turbulent boundary layer in 
whereas will here be chosen by fitting this expression to the meas compressible flow to those in a related incompressible 
require urements of Coles.’ The choice of .V, like that of // flow. The transformation is completely defined, for 
will be based on Table | and the transformed Reynolds 7 1.4, by the three relations 
nd pres Number (Rey); given by Eq. (10); the value of -V Lis - . 
1€ pres ippropriate in this case. Fig. 2 shows a comparison 2 pinta 
knowr i the experimental results with the expression (15 6 1 ./10)* 0 20) 
7 evaluated for : f ] ] ' jy 21) 
~~ [*=1 14) mo a 
en int The only difference between this and Stewartson’s 
hat the This choice seems good for flows with zero heat trans transformation appears in Eq. (21 
iber 11 fer and 0 < JJ < 5. The discrepancy with various Application of this transformation to Trucken 
6) ar ther expressions for /* given in the literature is only brodt’s simple formula’ for the growth of momentum 
thickness in incompressible flow, 
ex D 
D Y (Mee / m)'/* V + 1)/N] a(N)u,,- 8° tee d(x/D) 4 
e D 
mpres 6, D u,;D d d 22) 
casinane vields the result 
‘ent 1! 
so Se D) (T, 1 V+1) NJ a(V)(aD 1 » u-* wo 
* D 
MPTE™ (T/T) 9 d(x/D) + (@/D) (Te/Te)*h* 7" (M.,/ MY" 23 
enherr where D is a suitable reference length. The similar result obtained by applying Stewartson’s transformation is 
1ed by V+1) A , . \ 3) 
. 4 6/D) (T,/T7 y N + 1)/N] a(N) (aD /1 tal Piel ede: 
gm x 'D 
I; cor MPTe™ (T, Ta)4 d(x D) + (6D) (1. 7 dtl! TES DP ae 24) 
ertal! These ses ne : : , . ro “ee : , : “a : 
mee hese equations may be used, in particular, to predict the growth of momentum thickness on a flat plate in com 
pressible flow, where 1/, and 7, 7) are independent of x. The result obtained from Eq. (23) is 
enct 
Ox = Cr 2 =$((N +1) Nja(V)!* “*? (px ‘pe VDT (ey sl r./75" 25 
lay 
nd 
van Where Cr is the av erage skin-friction coefficient. This Cr Cpr;] Rex fixed My bm, = x 
Ws 0! expression mav be compared with the corresponding . \ \ or 
Also | One for incompressible flow to obtain the effect of com- ld <0 


| eae ’ ; a ae a ‘ 2 5 : 
e tw Pressibility on average skin friction in the form for equal values of Re, in the two flows. The agree 





262 JOURNAL OF 


ment with Coles’ measurements, shown in Fig. 2, is 


quite satisfactory. 

Finally, Eq. (25) has been checked by comparisons 
of its predictions with measurements made on the sur- 
faces of M = 2.5 and M 3.5 nozzles of the Naval 
Supersonic Laboratory's wind tunnel. In each case 
the momentum thickness was measured both at the 
throat and in the nozzle exit plane; the measured values 
at the throat were used as initial values in the calcula- 
tions. The measurements were also used to compute 
values of (Re,);, on the basis of which a value of NV = 5 
was chosen. The results are given in Table 2, together 
with percentage errors based on the measured values; 
the predictions of Eq. (24) are also tabulated. In each 
case the new transformation gives a satisfactory result, 
closer to the experimental value than the result of apply- 
ing Stewartson’s transformation. 
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6 (in Q% Error 
Measured 0.061 0) 
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Discussion of the Lecture 


(Continued from page 229) 


on spheres at Mach 10 in the hypersonic wind tunnel 
and fired them at 10,000 ft. sec. in the range to measure 
the drag and real gas effects, it was considered to be of 
only academic interest to a few research people who 
were not really interested in hardware since everybody 
knew, at that time, that hypersonic vehicles must 
have sharp-pointed noses to reduce the drag. Mr. 
Allen’s dissertation shows how fast and unexpectedly 
such results of systematic research can advance the 
solution of apparently insurmountable problems. His 
charts show in a few simple revealing curves the prin- 
cipal characteristics of long-range vehicles including 
those of their re-entry process. It is, for example, a 
startling fact that the maximum deceleration of spheri- 
cal bodies with weights varying within eight orders of 


magnitude, and entering with the same speed, does not 


vary at all. Only the altitude changes where the peak 
occurs. 

Once such general principles are demonstrated, de 
tailed calculations can immediately be pointed toward 
the missing links in the final solution. For example, in 
finding within the indicated limits the correct shape oi 
the body, the dependency of aerodynamicists, ballis- 
ticians, and materials and structures people on each 
other shows up as never before. Their arguments about 
the fraction of kinetic energy which is converted in heat 
and how much of it enters the vehicle are considerably 
simplified by having such charts as we have just seen. 

I wish for the sake of our future work that we had 
more of this kind of communication which puts in 4 
nutshell clear and precise information. 

Mr. Allen, I enjoyed your lecture. 














Bloc 


Tra! 


A Ce 

an 
The 
Rem 


AN 


Prog 





lanes of 


rison of 


19-676 


indar 


eaders 


Line 


ol. 22, 


Smoot} 
t Plat 
Insti 
t No 


rsoni 
Table 
ute ol! 
1954 




















Readers’ Forum 


RIEF REPORTS of investigations in the aeronautical sciences 
lished in the JOURNAL are presented in this special department 


to a maximum of 1,300 words or the equiyalent of 


illustrations 


Committee does not hold itseif responsible for the 


Publication is completed as soon as pos 


and discussions of papers pub- 
Entries must be restricted 
one JOURNAL page including formulas and 


The Editorial 


opinions expressed by the correspondents 


sible after receipt of the material 





Contents 


On Optimum Nose Curves for Superaerodynamic Missiles 
S. TAN 


Blocking in the Supersonic Wind Tunnel 
BAIn DayMAN, JR 


Transient Temperature Distribution in Aircraft Structures 
IRVING FRANK 


An Iteration Method for Solving Linear Problems in the 


Theory of Shallow Shells 
WILLIAM A. NASH AND P. L 


SHENG 
A Note on the Optimum Distribution of Material in a Beam 
for Stiffness B. SAELMAN 


Note on Hydrogen as a Real-Gas Driver For Shock Tubes 
PauL W. HUBER 


Calculation of Supersonic Flow Past an Axially Symmetric 
Cylinder... MartHA W. EvANs AND FrRANcIS H. HARLOW 


Combustion in the Laminar Boundary Layer of Chemically 
Active Sublimating Surfaces 


M. R. DENISON AND D. A. DOOLEY 


Effect of Loads on the Spring Constant of a Particular Type 
of Flexure FREDERICK G. BLOTTNER 


Hypersonic Heat Transfer to Catalytic Surfaces 


SINCLAIRE M. SCALA 


A Comment on the Methods of Calculating Natural Modes 
and Frequencies of Vibration SAMUEL SILVERBERG 


The Effect of External Sound on the Vortex Shedding From 
Cylinders DONALD I. GARBER 


Remarks on Systematizing the Calculations of the Natural 
Modes of a Free Structure... R. Kappus anp D. CLERC 


A Note on Isentropic Compressible Flow of Air With Vari- 


able Specific Heats kK. E. TEMPELMEYER AND L. SELF 
Propagation of Weak Waves in a Dissociated Gas 
F. K. Moor 


263 


?64 


265 


267 
2OS8 


269 


269 


271 


)T,) 


re 
“10 


ore 


-in 


~é 


276 


278 


OT 
209 





263 


On Optimum Nose Curves for 
Superaerodynamic Missiles 


H. S. Tan 
Consultant, Therm-Electric Meters Co., 
November 27, 1957 


I A RECENT 
missile-nose 


under the assumption of regular molecule reflection? is precisely 


N.Y 


Ithaca 


NOTE,' Dennis has pointed out that the optimum 


curve for superaerodynamic regime obtained 
that for the hypersonic regime derived according to Newtonian 
slip flow 

While 


namic and hypersonic missile-nose shape is obviously of impor 


the above-observed coincidence between superaerody 

tant significance, there is no such coincidence for the case of ran 
dom molecular re-emission, which, to certain extent, is actually a 
more realistic model. Thus, it seems that a more thorough sta 
tistical investigation of the actual manner of free-molecular re 
flection should be conducted next. On the other hand, it seems 
not without interest at this stage to compare the two curves cor 
responding to these two different assumptions. The integra 
tions are briefly indicated in what follows, together with tabu 
lated results and curves drawn to scale As already noted,‘ the 
fairing of nose tip beyond y’ = y/& does not merit discussion be 

cause this point can always be brought as close to the axis as we 
like by a simple change of scale 


The differential equation for the missile-nose curve is given by 

y” = —(y’?/y) (1 + y”)/(R y"" (1 

Its first integral with respect to y is obtained by noting y” = 
y'dy'/dy, and is given by: 

y = C1 + y2)*tP2/y 2 


Its first integral with respect to x is obtained by noting that y” = 
dy’/dx and that y is given by Eq. (2 Indeed Eq. (1) then gives, 


upon integration, 





Contributors Please Note 


Readers’ Forum items must be confined to the 
equivalent of one page in the Journal. Contribu- 
tions that exceed this limit will be returned to the 
authors for condensation and rewrite. To avoid incon- 
venience and delay, please adhere to the following 
specifications: 

(1) Five, double-spaced, typewritten, manu- 
script pages (8'/> by 11 in.), including 
wide margins, formulas, and headings 
equal one printed page. 

(2) For every illustration, deduct at least 
one-half or as much as one manuscript 
page, depending on the size of the 
illustration. 














264 JOURNAL OF THE 


AERONAUTICAL 











Fic. 1 
TABLE | 
(xo = 0, vo = 0.010 
k= 3 k 2 
yw’ = V3 Y = v2 y = V2 

V X Vv Y V x V 
V3 0 0 O10 
V2 0 0.010 0) O.010 
l 0.00255 0.0130 0. 0008 0.0109 
0.50 0. 0466 0.0406 
0.25 0. 668 0.235 
0.20 1.519 0.440 
0.15 $. 956 1.005 4.822 (0.973 0.759 0.177 
0.12 11.97 1.935 11.63 1.881 1.480 0.273 
0.10 24.69 3.810 24.00 3.220 2.560 0.391 
0.075 5.940 0.693 
0.05 391.3 26.10 380.5 25.39 20.55 1.540 
0.01 243,600 3,256 237,000 3,166 2,567 38. 50 

Cc 0.00325 Cc 0.00316 Cc 0. OOB85 

c 0.00314 C. 0. 00327 Cc. 0. 00441 

r= —( S (1 + yy?) YP? (Rk — y’?)/y' |dv’ + OG (3 


The resulting curves for these two respective cases are easily ob 
tained by quadrature as follows: 


Regular molecular reflection, k = 3 


Random molecular re-emission, Rk = 2 


It can readily be shown that the asymptotic behavior of Eqs. 
(4) and (5) are indeed given by Eqs. (9) and (11) of the previous 
note.* It is also interesting to note that the optimum nose cor 
responding to random molecular re-emission is substantially more 
slender than that corresponding to regular molecular reflection 
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Blocking in the Supersonic Wind Tunnel* 


Bain Dayman, Jr 

Jet Propulsion Laboratory 
Pasadena, Calif 

November 97, 1957 


A N IMPORTANT PROBLEM in the use of supersonic wind tunn¢ 
for research and development work is to determine, for 


given model configuration, wind tunnel, and Mach Number. ¢ 


Inetitrite at 
Stitule O . 


California 





maximum permissible model size. Tests were run in the Je; 
Propulsion Laboratory (JPL) 12-in. supersonic wind tunnel j 


order to obtain a reliable method of predicting the blockir 


characteristics of models during an attempted start —i. 


supersonic flow could be established in the test section for a par 


t 


ticular model (size and shape) at a given Mach Number Fig | 


* This paper presents the results of one phase of research carried out at tt 
Jet Propulsion Laboratory, California Institute of Technology, under ¢ 
tract No. DA-04-495-Ord 18 d by the De 


Ordnance Corps 


sponsore partment ol 














Fic. 1. Cone-eyvlinder model in a 12-in. supersonic wind tunn 
showing barely blocked flow 
| 
72 it cit T T ? 
Symbol o Test model 
a + 
Oo - Hemisphere-cylinders and spheres 
68 A, | ° 
4 146 Cone- cylinders 
° 19. 7° Cone- cylinders 
. a =| 276° to 366° | Cone- cylinders 
64 — “A = 
O Usk 
$y 2 
7 ° 20 = Cone-apex angie 
- © $i@ = Open symbols ~ start 
© ° Solid symbols ~ block 
5s 
a 9 
3 * 
2 56 
a 
5 © ° % + 
2 4 
= 
°o 
© se 
o 5 
g o 
= 48 | a. 
$s * 
= ~ 
a © | 
” Es 
40 
36 } } @ 
* 
re) 
32 
2 3 4 
Mach number 
Fic. 2. Uncorrected blocking data for indicated models ™! 


JPL 12-in. supersonic wind tunnel 





whet 


it 1s 
from 
tunn 


sepa 





| 
—— 
| 


TABLE 1 
Mach Number \rea Ratio 
p 1.5 0.0854 
se 20 0.1782 
Sf 0. 2808 
. 0.3276 
0. 4000 
tunne 
lor : 19 } 
i 5 for one model, an example of barely blocked fiow (7 = 
Ny : ) Four classes of models were used: spheres, hemisphere 
linders (3 to 4 in. diameter), cone cylinders (2 to 4.5 in. diam- 


30° to 73° cone apex angle), and a disk (3 in. diameter 





- | Previously, a rule-of-thumb value was used to determine the 
si maximum size of a model that could be run at a chosen Mach 

I 1. he value used was 60 per cent of the maximum one 

it nal theoretical area (based upon a norimal shock wave 

wcross the test section upstream of the model and sonic flow about 

eae he model The ratio of the model frontal area to the test 
Arn section area from this simplified theory can be expressed by the 


following equation: 


Model frontal area 1 Pu/ Pe 
Test-section area A /A* 
where 

A/A* ratio of local cross-sectional area of an isentropic 
stream tube to cross-sectional area at the point 
where the Mach Number is unity (one-dimen 
sional theory 

Pu,/t total pressure ratio (upstream to downstream 


across a normal shock wave 


Table 1 gives several values of this theoretical area ratio (from 
|  theequation for y = 1.400) 

Although the 60 per cent value of the theoretical model area 
is generally adequate for models having slender pointed noses, 
it is unconservative for blunter models. Fig. 2 presents the data 
from these blocking studies in the manner normally used in wind- 

| tunnel operation. The start and block conditions are normally 
separated by 0.02 in Mach Number. The ordinate represents 
the frontal areas of the tested models as percentages of the 
theoretical model areas. It is readily apparent, for the data as 
mne i group, that this method of presenting the data is of little or no 
use insofar as the prediction of blocking during an attempted 
start is concerned 
Fig. 3 presents the data of Fig. 2 with the bow-shock-wave 
strength instead of the Mach Number as the abscissa The 
shock-wave strength is defined as the perimeter average of the 
ratio of the static pressure increase across the shock wave to the 
free-stream static pressure upstream of the shock wave. The 
perimeter average is taken along the intersection of the bow 
1 shock wave with the test section walls, using the shock-wave 


| shape which would have existed had there been no wall inter 


ference rhe test section wall is taken to be at the displace 
ment boundary layer on the physical wall. Hence, the test 
section area is less than the physical area by the area of the dis 
placement boundary layer It is this smaller area upon which 
the data are based \n arbitrary straight line was drawn in an 


ittempt to separate the starts (open symbols) from the blocks 
solid symbols Except for the disk model, the data correlate 
to about 2 per cent in model diameter or in Mach Number 
whichever gives the greater latitude from this blocking correla- 
ion line \t this time, no final conclusions have been reached 
i correlation method which will include disk models 

he ability to start the tunnel with a particular model does not 


Tarrant 


gua itee a successful test The wake shape mav be deformed, 
t the bow shock wave may reflect back onto the model 


\lthough this correlation pertains specifically to the JPL 12 


In. supersonic wind tunnel and to the types of models tested, 
similar correlations probably exist for other supersonic wind 
i tunn ind other types of models. The closer the operating 


conditions of the various tunnels (such as supply pressure, dif 
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FORUM 265 
ae ae : az 4) PT | 
4 Symbo ° Test mode 
} ee : , ; 
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68r & a e 
a 46 Lone -cylind 
a ° 197 Cone-cylin 
64} + . 276 to 366° Cone-cylind 
- » Disk 
4 =~ ° 
o 60} -—_$ = te; se. { 2e Cone-apex ongle 
© PN Open symbols ~ stort 
3s | a ) Solid symbols ~ block 
2 56 wt 
® ° hes ° 
S P 
- . @ 
= 52 
S | . . 
S 
+ . a 
Ss 48} ‘ 
; + 
4 
® , , i 1 
é a4, a 
40! 
| =) 
%} 
7 
321 
° 02 04 06 08 ° 2 4 6 8 20 22 
Shock strength 
Fic. 3. Corrected blocking data for indicated models in JPL 


12-in. supersonic wind tunnel 


1 


fuser efficiency, compression ratio), the closer the correlations 
may be. Preliminary checks in the JPL 20-in. supersonic wind 
tunnel support the possibility of extending the correlation to 


other similar wind tunnels 
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Transient Temperature Distribution in Aircraft 
Structures* 


Irving Frank 

staff Mathematician, Allied Research Associate 
Boston, Mass. 

November 25, 1957 


INTRODUCTION 


Sse S APPROXIMATE numerical methods are available for 
calculating temperature distributions in aircraft sections ex 
posed to high heat inputs.’ Such numerical procedures are 
used because of the complex geometry and boundary conditions 
which apply to the ordinary aircraft section rhese methods, 
however, usually require the use of computing equipment, and 
furthermore, lack the generality of analytic solutions. It is, 
determine the temper 


therefore, frequentiy advantageous to 


ature distribution analytically by making appropriate simplify 
ing assumptions. The accuracy of the solutions will in these 


cases depend on the validity of the simplifying assumptions 


\NALYSIS 


Consider an aircraft cross section (Fig. 1) exposed to a heat 
input on one surface or to a symmetrical heat input on both 
surfaces. It is assumed that: the structural members are 
thermally thin so that the one-dimensional heat conduction 
equation is valid in each, the mechanism of joint heat conductance 


can be represented analytically by convective heat relations, 


* This research was performed for the Structures Branch, Aircraft Labora 
tory, Wright Air Development Center, under Contract No. AF 33(61¢ 
5001 
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Returning to the channel section, and assuming tl 














aX X—e x—= ature of the flange at the junction of the flange and t Wwe 
REGION ID {REGION Z REGIONI given by Eq. (6) with y = 0, the temperature distributi in t 
. 3 ce 3 remainder of the section can now be determined 
.—— The heat-balance equations for the four regions (Fig, ] 
| SKIN come 
REGION O*u / Ox" l/a Ou/Ol F(t)/akK Region | 7 
FLANGE 070 /Ox? = (1/a) (O#/Ot F(t)/ak | + 
REGION Vv \ G(t)/aK Region I] R 
‘wes 020/08? = (1/a) (Ou/Odt F(t) /aK Region II] 
O7w/Oz? = (1/a) (Ow/OF) G(t)/bK Region I|\ 
where G(t) is the heat flux from the skin to the flange and wil] } 
| determined by solving the above equations 





The initial and boundary conditions are 





3 Jt | = 





. ; 7 . d it? =O: uF uTHTKHte HKHEw=d 
Fic. 1. Aireraft cross section ' sf 0 g=i: u=%, Bn/dx = du/Oz 
f) att =0, F=0: H=u, O/dF = Ou /OX 
the increase in temperature of the internal structure due to irradi \f dee aa ina = finite, u is finite 
ation from the skin is negligible, and the heat transport proper (4) ats = 0: weve [hak p)/aK* p' qt | 
ties including the joint conductivity are independent of tem ian = dv/dy) = L~'[ghaF(p) /ak* 
perature i) atz =71: (w/dz) = 0 
In addition, it is assumed that the heat flux to the web is inde where L~! is the inverse Laplace Transform and u is the temper 
pendent of the existence of the flange —i.e., the heat capacity of ature in Region I, @ is the temperature in Region II, uw is tl 
the flange is negligible. Thus, the temperature distribution in temperature in Region III, and w is the temperature in Region 
the web can be determined by considering the temperature dis IV. The coordinates for the various regions are shown in Fig 


tribution in the ‘‘T”’ section. The heat balance equations for this Taking the Laplace Transform of Eqs. (7)—(10) yields 


section are 


d2u*/dx? = (p/a)u* — [F(p)/aK Region | 1] 
O71 /Ox? = (1/a) (Ou/Ol) F(t) /aK | ( 
, d*ii*/dx? = (p/a)i* [F(p)/ak 
077 /Oy? = (1/a) (dv/Ot (2 G(p)/akK Region I] 12 
where w and v are the skin and web temperature, respectively, 4 
is the distance from the web on the skin, y is the distance from the d2u*/dk? = (p/a)u* — [F(p)/ak Region II] 
skin on the web, a is the skin thickness, a is the skin diffusivity, : . : . 3 
K is the aircraft material conductivity, / is the time, and F(t) is d*w*/d2* = (p/a)w* — [G(p)/bK Region I\ ld 
the heat input to the skin Solving Eqs. (11) to (14) with the boundary conditions (d 
The boundary and initial conditions are: yields 
(a) att = 0: nu =v =0 u* = T(e7# le-@ + [aF(p)/uKp Region | 15 
(b) atx =O, y = O: 2aK (Ou/Ox) = hb |u 7 |, 7 " 
bK(dv/Ov) = hb [wu — 2] a* = Te-® +e” + [aF(p)/aKp 2r Region II If 
(c) atx = yv= : wis finite, v is finite u* = (e~@ lhe" + [aF(p)/aKp Region III 7 
where /} is the web thickness and h is the joint conductivity. wt = bV/aje"@ + e-@ 21-2 rn 200 /b Region I\ 18 
Taking the Laplace Transform of Eqs. (1) and (2) vields 
where = [habF(p)/a*?K2piq + y) (1 e~<9 
d?2u*/dx? = (p/a)u* — [F(p)/aK] 3 
d*y*/dy? = (p/a)v* (4) Upon setting [1/(1 e-*%)] = 1 + e-70 + 49 


. P ‘ — : ' it is found that the solutions for the temperature distributions 
where F(p) is the transform of F(t), and u* and v* are the trans : : 
J ; oe a , Hs : are all of the form 
forms of «and v. Solving Eqs. (8) and (4) and applying bound 


ary and initial conditions (a)—(c) gives 


e~@ 
— — ‘ ' : >. (p) + laF(p)/aKp 19 
u* = hbaF( p)/2a*K*(q + y)Je~™ + [aF(p)/aKp] (5 a A ara a F(t laF(p) /aKp |e . 
1 L\qG tT 
= |haF(p) aK? p(q + +7)] e~% (6) 
; where € = 1 in the skin and is zero in the web and flange, and tl 
here y = h(b + 2a)/2ak nd gq=WVp . 
where S od) / ad a pa ey e/a A; are independent of pand x 


From reference 5 
l 


L if ¢ ¢ gq” l(q T y)] - la .~ vy)” Va - erfe (x 2 at ) , 7 av ‘at | = b» [—2, Vv afi 1” erfc (4 2 xl } 2 


0 


2/ 
os A [ qxi /g” l(g H 7)] F(p) = [a (—y)"] | z A, Sons + ary” erfe [ (x 2V/ ar) + vv er] 


1 ] 
n— 1 } 
/ r r 1 
z, | —-27V ar | i’ erfe (w/2 ar) ¢ ) F(t tr), dr (2l 
? 0 
If F(t) is any polynomial in ¢, then Eq. (19) can be inverted directly. For let F(t) = ay + ait + aot? +... then F(p) = (ao/p) t 
(a;/p?) + (2a2./p3) +... = (ao/ag?) + (a:/a*q*) + (2a2/a%g*) +... and Eq. (19) can be inverted by using the inversion formula 


in Eq. (20) 
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2. Temperature distribution correlation of a typical panel 


The functions appearing in Eq. (20) are tabulated in reference 


jJand ¢ 
EXPERIMENTAL CORRELATION 


The preceding analysis has been employed to correlate tem 
rature distribution data obtained in experimental tests held at 
e U.C.L.A. thermal facility, reference 6. Fig. 2 compares the 
sperimental and calculated temperatures for one of the test 
panels rhe joint conductivity for this panel was found by 
matching the experimental and calculated temperatures at one 
point. This point is indicated by an “‘X” in Fig. 2. As may 
be seen from the figure, the agreement between the calculated 


ind experimental results is quite good. A full discussion of the 





pt lure followed in performing the temperature calculations 
nd additional experimental correlations are contained in 


reference & 
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An Iteration Method for Solving Linear 
Problems in the Theory of Shallow Shells* 


Nash and P. L. Sheng 
and Research Assistant, Respectively, Department of 
Engineering Mechanics, University of Florida, Gainesville, Fla 





mber 2, 1957 


| HE BEHAVIOR of a thin elastic shallow shell is described by 


» simultaneous differential equations for the deflection « 


ind the Airy function F of the membrane stresses. These equa- 
tions, due to Marguerre! and Reissner? are 

rhese results were obtained in the course of research conducted under 
Grant ( 92 from the National Science Foundation 


FORUM 207 


V°v2F = Eh (2(072/0x0y) -(07w/Ox0y 
07s /Ox*) -(07w/ Oy" 075/04 ) Ox?) ; 
DV2V%w = (02: /Ox?)-(0?F/Oy’ 
022 / Oy?) -(07F/ Ox 2(07s/Ox0y)-(07F/Oxdy) + p (2 


where V2 denotes the Laplacian operator 





rdinates vy) denotes tl 


Here x, y are Cartesian cor 1¢ 
equation of the middle surface of the shell, / thi ell thickne 
and D = Eh'/12 (1 v?) where — is Young’s modulus and » 
represents Poisson’s ratio The applied load f intensity p, 
icts in a direction normal to the x-y plane \ shell is considered 


to be shallow if the ratio of its height to a representative base 
dimension is less than approximately 1/8 
9 


\n iteration-type solution of Eqs. (1) and (2) is possible if the 


equation of the middle surface of the shell can be written in the 


form 

=A Yr, ¥ 
whert x, y) denotes a dimensionless function, and \ represents a 
dimensionless shallowness parameter \lso, & denotes some 
representative base dimension of the shell. If the characteriza 
tion of the shell is such that the maximum value of f (x,y) is 


unity, then AX < 1/8 for shallow shells. This representation of the 
shell was first introduced by Nazaroy 
To perform the iteration, we expand the desired unknowns 
and F in terms of powers of the small dimensionless parameter 


\, namely, 


where w, @2, Fi, Fs, ete. are functions of x,y but are independent 
) 


of X. Substituting these values into Eqs. (1) and (2) and equat 


ing terms containing like powers of A, we obtain the set of equa 


tions 
v2V2F 0) 6 
DV 2V 2% p 7 
AV 20 2F, = Eh} 2(02s/Oxdy) -( 07x /Ox04 
02s /Ox?) + (07) / Oy? 0*s/ Oy 07%; /Ox*); S 
DV 20 we = (07s /Ox*)-(07F,/Oy?) 4 
072 /Oy?)- (07 F, /Ox? 2(072/OxOy)- (07 F)/Ox04 9 
etc 


\dditional equations may be added by equating the squares 
and higher powers of A, and, in each, the right-hand sides will be 
functions of previously determined expressions in the series, 
Eqs. (4) and (5 Thus at every step it is only necessary to solve 
a pair of uncoupled equations for F; and with respect to bound 


ary conditions. In this manner the solution to Eqs. (1) and (2 


may be obtained to any desired degree of accuracy \lthough 


the linear theory of shallow shells is derived on the basis of neglect 
ing, in comparison with unity, terms involving squares and 
higher powers of A, there is no justification for discarding such 
terms in Eqs. (4) and (5 

Further it is of interest to note that in certain cases the bound 
ary conditions may also be treated independently rhis is 
true, for example, if the edge deflection and the edge bending and 
twisting moments together with the normal and shear forces 
along the boundary are prescribed. In this case, the boundary 
conditions regarding deflection and moment may be expressed 
in terms of w only while those regarding normal and shear forces 


may be formulated in terms of F only 
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A Note on the Optimum Distribution of 
Material in a Beam for Stiffness 


B. Saelman 


Lockheed Aircraft Corp., Burbank, Calif 


Design- Weight Engineer 
December 2, 1957 
SYMBOLS 
J torsion constant 
A angular rotation 
7 applied torque 
G modulus of rigidity 
1 enclosed area 
t wall thickness 
HW total weight 
p material density 


perimeter 
length of beatin 

A Lagrangian multiplier 
shear stress 


| volume of material 


T’ HAS BEEN STATED by some authors! that the maximum stiff- 
ness of a beam for a given weight is attained when the strain 
energy is a minimum or when the stress is constant. It can be 
shown that this condition does not, generally, result in maximum 
stiffness; however, it closely approximates the optimum condition 
in some cases 

The optimum distribution of material for torsional stiffness of 
tubular beams is such that the thickness is constant around any 
cross section, and, for positions along the axis of the tube, the 
wall thickness should be proportional to the square root of the 
torsional moment and inversely proportional to the enclosed 
area. For maximum bending stiffness, the effective flange thick- 
ness should be proportional to the square root of the bending 
moment or the square root of the product of the moment and the 
axial length on the beam, depending on whether angular or linear 
deflections are being considered. 

For the case of torsion, stiffness is measured by the amount of 
angular rotation, and, if ¢ is independent of s, is given by: 


¢—= 


Gf 7 Ndx = (1/4G){ (Ts 1°t)da l 
where 7, s, A and ¢ are functions of x 


The weight of the tube can be expressed by 


pf, sds 2 


Considering the weight fixed, the minimizing condition for @ is 


W = 


found by application of The Calculus of Variations: 


(0/dt) [(Ts/A*t) + Ats] = 0 (3 


Performing the operations in Eq. (3), solving for ¢, 


(W/Ap K T i. (VTs Ads | (4 


(p/4GW) [foo Ts Ajax | (5 


The strain energy due to twist can be written 


(1/G f, 7 J)dx = (1 26) |  Ptds = 
(1 2G) J, rav (6 


The minimizing condition for U with W fixed is 


and elimi 


nating by Eq. (1), 


Q=(1/A)vVT/rxX = 


i, = 


ie 


(0/ot) [f%ts — As] = 0 (7 
This leads to the conclusion that the stress, f, is constant 
f=vVx 


and that the strain energy per unit volume is constant 


dU/dvV =f? 


2G (9 
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With the aid of Eq. (2) we 


Substituting Eq. (10) into 
4, = 


p IGH 


It is observed from Eq. | 


t, is proportional to +/ 7 for maximum stif 


/2AVYVA = 


APRIL, 195 


can write 


Eq. (1) we find 


+) and Eq. (10) t 


ha 
Tr 


to 7 for minimum strain energy. 


t the wall thickness 


less and proportior 














Comparing Eq. (5) and Eq. (11) it will be found that, in g 
eral, 0; < 6 except 6; = 6 when the torque is consta 
Fig. 1 illustrates the case for 
Ary = Ap(x?/L?), s = Sr(x/L), T = T; 
inn ii 
—— 
L —r 
SKETCH | 
We have 
= BM, M =4il[R R-1 
and 6.= BN, N = llog(|R/(R l 
where B= pl rSpr*L/4AGWApr? and R = L 
Factors JJ and .\ in Fig. 1 are measures of the relative stiffness 


for the two cases 


become a and L 


The limits of integration for the 


integrals 


If ¢; is a specified strength level and ¢ is the unknown value t 


be determined for stiffness requirements, then if ¢ t T 


7 is obviously a minimum when ¢ 


is a 


latter value is determined by Eq. (4 
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1936 
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Note on Hydrogen as a Real-Gas Driver For 
Shock Tubes 


= Ww 
30 Thermodynamics Section, Langley Aeronautica 
Lat NACA, Langley Field, Va 


y there has been renewed emphasis on the use of non- 


EC! i 
R mbustion shock-tube drivers of low molecular weight 


use of the lower shock attenuation observed than with the 


tically more desirable combustion drivers Nearly as 

k Mach Numbers are obiained when attenuation is 

inted for and more uniform test conditions result from using 
mbustion drivers 

In reference 2 the maximum shock velocities observed along a 

tube were compared with the theoretically computed shock 

wities for two pure gas drivers other than air—namely, helium 


vdrogen Phe computed velocities were based on real air 


the driven gas and on ideal gas drivers. The helium-driven 
locities were significantly higher than theory, while the hydro 
iriven case resulted in values about equal to, or slightly 
er than, theoretical shock velocities 

t ld be expected, intuitively, that these two drivers would 
exhibit more similarity since the velocity of sound in either gas 


much greater than in air, and also in view of the similar shock 


ttenuation characteristics found in reference 1 and the similar 
driver transport characteristics noted in reference 3 for these 
is compared with other gas drivers 

vortant difference between helium and hydrogen, how 





may be found in the validity of the ideal-gas assumption 
is very close to a truly ideal gas both as to low compressi 





nd constancy of specific heat, while hydrogen is far from 





‘ in both of these respects in the temperature and 
ssure range of interest for shock-tube drivers 
mputation of real-gas hydrogen driving real-gas air has 
b made, and the resulting diaphragm pressure ratio versus 
ck Mach Number is plotted in Fig. 1 for two cases of driver 
ressure, / rhe curve labeled p; = 0 is shown to illustrate the 
rge departure from an ideal gas due to the variation of specific 
with temperature The other real-gas curve bearing the 
driver-pressure labels is computed for about the driver 
range used in reference 2. The latter curve was actually 
Db, = 50, 100, 


pressurs 
btained by first plotting three real-gas curves 
id 125 atmospheres), and judiciously fairing through these three 
ves to represent an approximate theoretical case for the driver 
conditions of reference 2. The computation was made using 
Riemann's unsteady-flow characteristics relation where the quan 


tity conserved (see reference 4) in the driver gas is o + u, where 


> 
o> F ap dp 


| this quantity is evaluated along an isentrope for each driver 
pressure using real-gas values for velocity of sound and density 
of hydrogen from reference 5. The corresponding pressure and 
fluid velocity are then matched across the discontinuity to the 


iir shock values, as in reference 2, to obtain diaphragm pres 


It is also seen from Fig. 1 that a pronounced effect of pressure 
upressibility) results when comparing the py = O curve with 
the points indicated for finite driver pressure 
interesting observation is that the effect of compressibility 


\y 
\Tl 


s to partially compensate for the effect of variable specific heat 
fhe net result of the computations, however, is to show that 
iwdrogen is not so unlike helium as a driver gas with respect to 
the shock Mach Numbers obtained in comparison with theory, 

real gas is assumed. The reasons why both gases yield 
significantly higher than theoretical shock velocities is still un 
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Calculation of Supersonic Flow Past an Axially 
Symmetric Cylinder* 


Martha W. Evans and Francis H. Harlow 

University of California, Los Alamos Scientific Laboratory 
Los Alamos, N.M 

December 9, 1957 


W'! HAVE USED a numerical calculation method! to deter 
mine the axially symmetric flow pattern about a flat-nosed 
cylinder moving parallel to its axis at Mach Number 17, = 1.58 
along the axis of a rigid cylindrical chamber. The approach to 
steady state was studied in detail for the case in which the cylin- 
der, starting from rest, instantaneously attained its fixed speed 
The air was represented by a nonviscous polytropic gas with 
specific heat ratio y = 1.4. The calculations were performed 
on an IBM Electronic Data Processing Machine, type 704 

An extensive discussion of the general computation technique 
has been given by the authors in reference 2. The gas is repre 

* This work was performed under the auspices of the United States 


Atomic Energy Commission 
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Fic. 1. Air configuration at steady state, showing stream- 

lines, detached-shock position, and sonic line. Dashed line is 

experimental shock position measured by Marschner.*  Hori- 


zontal projections were added to streamlines to show points of 
departure for determining shock position. Position of shock 
reflected from chamber wall is indicated by flow line deflections. 


sented by Lagrangian mass points, called ‘“‘particles,’’ moving 
through an Eulerian computational mesh. The flow region 
studied is adjacent to the front of the cylinder, as shown in Fig. 
1 For the computation, this region was divided into a mesh of 
14 X 22 toroidal cells of square cross section. The coordinate 
system was at rest relative to the moving cylinder 


ticles representing the air were arranged initially with four per 


The par 


cell. As the air compressed against the cylinder face, the particle- 
number average rose, reaching about six per cell at steady state. 
In order to achieve an incoming uniform density, the particles 
were given masses proportional to their radii. The left boundary 
was far enough upstream that the flow across it was unperturbed. 

The only novel procedure not discussed in previous reports was 
the technique for stabilizing the finite difference equations in the 
stagnation region. It is assumed that every velocity vector 
should point to the right and outward, and whenever any com- 
ponent violated this assumption, it was set to zero, and a corre- 
sponding change made in internal energy to preserve total 
energy. 

The configuration at steady state is shown in Fig. 1. The 
flow lines are drawn from particle positions. The calculated 
position of the detached shock is determined by observing where 
each flow line departs from horizontal. Toward the axis, this 
becomes difficult to observe, and the shock position is determined 
by the point, along a horizontal line, at which 1/? has dropped 
to (1/2) (W,* + M,?), where VW, is the Mach Number just behind 
the shock at the axis, as calculated from normal-shock relations 

Close to the cylinder face, the flow lines are poorly represented, 
and their positions are drawn only approximately. The central 
five flow lines tended to coalesce at the cylinder shoulder and 


thereafter to follow a single flow line. 
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The sonic line position was determined as the locus of pojnts 


at which the interpolated cell-wise Mach Numbers were | 
intersection with the shock wave was placed at the point whe; 
the latter makes an angle of 61.7° with the cylinder axi S pre 
dicted by conical shock theory 

There is some indication of a reflected shock from the chamber 
wall. Its position was determined by the flow-line deflection: 


indicated 

The drag coefficient for the front face of the cylinder was 
calculated. As the flow approached steady state, the drag evef 
cient approached, asymptotically (except for small oscillations 
the value Cp (nose) = 1.48. This value is about 8 per cent 
higher than the mean value given by Hoerner. 4 

The same calculations have been performed with several other 
Mach Numbers 


study was made at VM = §8 


To compare with the work of Oliver,® on 
While the nose dr ig agreed 
closely with his, the calculated detached-shock position stood 
outside of his by an average distance which was about 5 per cent 
of the cylinder diameter. The discrepancy could easily be 
tributed to the mesh coarseness, since the actual shock should 
stand less than four cells away from the cylinder face 

The calculations indicate that the time required for the lo 
flow pattern to reach steady state is given roughly by ¢ = d 
CV M 
sound speed, and .\/ 


where d is the cylinder diameter, C., is the free-stre 
is the flow Mach Number 

The experience gained in these calculations shows that a nun 
ber of improvements could be made 

1) The calculation mesh could be made finer. This would be 
at the expense of additional calculation time, but would increase 
the accuracy and fineness of detail. Approximately two calcula 
tion hours were required to reach steady state in the problem re 
ported here, and this increases approximately as the inverse cube 
of the length of a cell side 

2) The incoming air density and velocity could be varied as 
a function of time to represent passage of the cylinder through a 
variable atmosphere and its slowing due to drag. 

(83) The equation of state of the air could be changed to include 
various real gas effects which occur at high Mach Numbers 

$1) The evlinder could be terminated with a generalized nos 
shape 

The first of these improvements is of considerable impor 
tance, especially with regard to definition of the shock. The cri 
teria used to determine shock position are reasonable and yield 
good results, but the shock line as drawn really represents the 
center of a band of rapid variation which is spread over a thick 
ness of several cells. Likewise, lines of constant Mach Number 
are impossible to locate accurately near the shock, and their 
points of intersection with the shock must be determined by other 
means. As discussed by the authors in reference 2, it would seet 
necessary to require the computational cells to be much smaller 
than any significant dimensions in the problem. It is, therefor 
remarkable that many average functionals calculated in this and 
various other test problems agree as well as they do with know 
results 

Calculations are now being prepared to include several of tl 
above features 
plicated shapes is being undertaken, and the calculation mes 


In particular, the generalization to more con 


is being made finer. 
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i pointe | Combustion in the Laminar Boundary Layer of 

_ Chemically Active Sublimating Surfaces* 

t whe 

1S Dr MM. R. Denison and D. A. Dooley 

tems, Inc., Glendale, Calif 
lamb ombe 3, 1957 
flectior 
DETAILED ANALYSIS has been carried out in reference 1 for 
is al A th it and mass transfer from sublimating surfaces where 
g Coe sub] ted material subsequently reacts with the atmosphere 
lati b ry-laver equations were transformed in the conven 
ner ce nanner and the effects of streamwise pressure gradients 
ected in accordance with the argument of Lees, reference 2 

il othe f the diffusion coefficients of all species present are equal, if 
er® on = | 1, and if pu = constant, then the heat and mass 
‘pree sfet 1 be estimated with very little detailed knowledge of 
ste chemical reactions in the boundary layer. Under these 
er ce nditions the mass fractions T,, of atomie species m, irrespective 
be nolecular configuration, and the total enthalpy /,, including 


the chemical enthalpy, can be expressed as linear functions of the 


ratio u/u, where u, is the velocity at the outer edge of the 








‘ rs Fast combustion reactions of the type 
stre vy’ f + v9'O — vp,"Pi + vp,”P2 l 
: well approximated by a thin reaction surface or flame front 
— vhere both oxygen, O, and fuel, f, are completely consumed by 
mbustion to form products, ?; and PP. On one side, cf. Fig. 1, 
uld be | ‘ , 2 pe : 
kvgen is mixed with products and inerts, while on the other 
— ide, fuel is mixed with products and inerts. The combustion 
alcula surface may be located either within the boundary laver or at 
sepa the solid surface. Dissociation and recombination reactions 
e cube : 6 ‘ 
for oxygen and nitrogen are also included. 
The heat flux, go, conducted into or radiated from the body at 
ted as the surface is found to be given by 
ugh a 
da = (PetteCje/2)( Ise + ToeQpow — Br he 2 
clude where I',, is the mass fraction of oxygen at the edge of the bound- 
laver, Q is the heat of combustion at wall conditions pet 
1 nost unit mass of oxygen, Ay, is the heat of sublimation, and h,, is 
the enthalpy per unit mass of air evaluated at the wall condition 
inpor f temperature and degree of dissociation. It should be noted 
1€ Cri that although the form of Eq. (2) is invariant with the physical 
yield model assumed, the numerical values of the last three terms are 
i's the somewhat model dependent. The quantity B represents a ratio 
thick of mass transfer to skin friction or heat transfer: 
mber 
their : 
rhe authors wish to express their appreciation to Dr. F. E. Marble 
other Professor, California Institute of Technology, and Consultant, Aeronutronic 
seen Systems, Inc., who has followed the development of the analysis closely and 
taller ntributed many valuable suggestions, 
fore, 
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Fic. 1. Model of combustion process in boundary layer 
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Fic. 2 (top Effect of wall temperature, pressure, and Mach 
Number on heat transfer for surface combustion 

FG. 3 (botton Effect of wall temperature, pressure, and flow 

parameter on heat transfer at free-stream Mach Number .1/ 10 

for burning rate assumptions of reference 1 


> _ 9 5 
B -/h Pell else » 


where mm, 1s the mass rate at which the surface material enters 
the boundary layer. When combustion takes place at the body 
surface ‘. B = §* = Do v 'M i VW 

For a given value of the wall temperature the combustion 
reaction will remain at the surface so long as the oxygen ts able 
to diffuse to the wall at a rate in excess of that required to oxidize 
all of the surface material sublimed. If the wall temperature 
becomes sufficiently high, however, the surface material will sub 
lime at such a large rate that the combustion zone will leave the 
surface and both oxygen and fuel will diffuse toward that zone 
If m,. is the sublimation rate at a given surface temperature, 
then the reaction remains at the wall if m,. < my, while the com 
bustion zone is somewhere in the boundary laver away from the 
wall if my ,, = My 

The skin friction can be related to B through the boundary- 
layer solutions of Emmons and Leigh rheir results can be 


fitted empirically by the expression 


(Mu/ Pelle) WV R.| + 
0.824 [Cre Pelle) WV R.| 4) 


Cr V Re = 0.664 — 1.388 


where the effective Reynolds Number for blunt bodies is given by 


fe 
Re = ( i Pele eh y*ds / pe? x) 
vV 





The quantity & is unity for bodies of revolution and zero for two 
dimensional cylindrical bodies. 

As an example, the analysis has been applied to the sublimation 
and combustion of graphite surfaces inasmuch as the physical 
and chemical data required for numerical calculation are readily 
available for graphite. These data are summarized in reference 
1. In the case of a high temperature environment such as that 
generated by hypersonic flight, the graphite oxidizes principally 
to carbon monoxide: 

C+0—-CO 


The manner in which the free-stream Mach Number, wall tem- 
perature, and static pressure affect the heat-transfer rate is indi- 
cated in Fig. 2 for the case of wall combustion. It is to be noted 
that over most of the temperature range considered, the heat 
transfer parameter is actually less than it would be for an ordi 
nary boundary layer at zero wall temperature. Evidently, 
the negative potential due to the combination of wall enthalpy, 
he, and heat of vaporization, BA;y~, more than offsets the increase 
in positive potential due to combustion, TyeQpow. The value of 
B* in the case of graphite is so small, B* = 0.17, that the decrease 
in skin friction and in heat transfer due to the effect of mass 
transfer on the velocity profile is rather negligible for surface 
combustion. Incidentally, g2 would rarely be negative 

Plots similar to Fig. 2 could be constructed for larger values 
of B, corresponding to gas phase combustion. As B increases, 
the whole level of heat-transfer rates shifts downward until 
finaily the entire plot is in the negative region. A graph showing 
the complete range of operation at a free-stream Mach Number 


of 10 is presented in Fig. 3. 
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Effect of Loads on the Spring Constant of a 
Particular Type of Flexure 


Frederick G. Blottner 
Staff Member, Sandia Corporation, Albuquerque, N.M. 


December 16, 1957 


INTRODUCTION 


_gesaconee PIVOTS are frequently used in wind-tunnel instru- 
mentation. One special use is to support models in wind 
tunnels when static and dynamic stability coefficients are being 
determined. In order to determine the static stability coefficient, 
it is necessary to know the spring constant of the flexure when 
the wind tunnel is at the desired air speed. Sometimes it is as- 
sumed that the flexure spring constant is independent of the air 
speed or normal and axial forces; therefore, the spring constant 
at zero air speed (normal and axial forces equal zero) is used. 
However, the spring constant can change appreciably with 
normal and axial forces and serious errors in the static stability 
coefficient can occur. 

The effect of loads on the spring constant of one type of flexure 
has been considered by Wittrick,’ but another type as shown in 
Fig. 1 was not considered. This type of flexure has been used by 
several organizations doing dynamic testing. Therefore, for this 
type of flexure, the analysis of Wittrick is extended, and cor 
roborating experimental data are presented in this paper. 


ANALYSIS 


Wittrick has written the relation between torque and rotation 
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as 


where 


7 


Vg 
K,b 


In reference 2, it was shown that for the typical flexure of the 
type considered in this paper, the coefficient 6 is very small and 
can be neglected (see Fig. 2 which is from reference 2 Ther 


) 


= applied torque about the point at which the und 
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T = El/L (K0 + be? + 


flected strips cross, in.-lbs 
angle of rotation of the moving end of the pivot, rad 
Young’s modulus of the material, Jb. /in.? 


moment of inertia of the cross section of the strips 


in.4 
length of each flexure strip, in 


dimensionless coefficients 


) 


the spring constant in in.Ib./rad. is 
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Fic. 1. Typical crossed flexure pivot. 
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Effect of neglecting second-order term on relation 
between torque and rotation. 
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Fic. 3 Effex t of n¢ glecting higher order terms on Eq. 4 
Fic. 4 (botton Effect of axial force on spring constant 


k = EIK/L 2 
e value of A as given by Wittrick is 
A = B, (coth B Bi) + Bs | 
where for this flexure 
= NL;?/4E\l,, 62% = AL2/4F21 


| is axial force, .V is normal force, and subscripts 1 and 2 refer 


vertical and horizontal strips respectively. Using @. = 


| substituting Eq. 3 in Eq. 2, the following is obtained for the 


Exte72/L2) (cot y2 + v2) (4 


[wo points should be noticed about this equation. Since the 


mal force .V can be either positive or negative, the first part 


I Eq 


$f can change from hyperbolic to trigonometric. The 

second point of interest is the value of the spring constant when 
Vand A are zero, This is 

ko = (E,\l/L1) + (Eele/ Le 5 


$4 can be derived if the cot and 


nore useful relation than Eq 
th are expanded in a series and terms of order 2 .. ° and 


greater are neglected. The result of the above procedure is 
ke=hkh+(1/6)(AL NL, 6 

see the effect of neglecting higher-order terms, sec 

t and 6 are plotted for the typical flexure 


In order to 
Fig. 3 where both Eqs 


see Fig. 2 for dimensions of flexure 


COMPARISON OF THEORY WITH EXPERIMENT 


rhe change in spring constant was determined experimentally 
fora flexure of the following dimensions: 
= () 250 in D = (0.500 in 
0.233 in Le = 1.234 in 
0.015 in to = 0.025 in 
Fy = 30 X 108 psi E, = 17 X 10 psi 


In the « xperiment, se veral bodies of different weights, but of the 
same moment of inertia, were placed on the flexure to simulate 
Phen the periods of oscillation of the different bod- 


axial loads 


determined, and, from this, the change in spring constant 
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Eq 6) and ex 


From this figure, one can 


was calculated \ plot of the theoretical curve 
perimental points is shown in Fig. 4 
see that there are large changes in the spring constant as pre 


dicted by the theory 
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Hypersonic Heat Transfer to Catalytic 
Surfaces* 


Sinclaire M. Scala 

Research Engineer, Missile and Ordnance Systems Dept 
General Electric Co., Philadelphia, Pa 

December 17, 1957 


INTRODUCTION 


‘es DETERMINATION Of hypersonic heat transfer to fully cata 
lytic and noncatalytie surfaces has been considered previ 
ously by Lees! and by Fay and Riddell.? This note will clarify 
these limiting cases and will also present some results of a re 
cently declassified study In Fig. 1 are indicated the various 
physicochemical phenoniena which occur in a dissociated gas 
in the mean free path (shown to exaggerated scale) near a surface 
Since atom recombination is only one of several simultaneous 


processes, the catalytic efficiency, defined by 


atomic collisions effective in recombination 


. total number of atomic collisions 


is a gross quantity which includes the effects of the intermediate 


steps. Thus, a treatment is required which reconciles a con 


tinuum and a molecular approach 


SYMBOLS 


{ mass fraction of atoms 
specific heat 
similarity stream function 


enthalpy, including chemical 


I Lewis Number 
I thermal Lewis Number 
normalized product of density and viscosity 
* Presented at The Ramo-Wooldridge Second Symposium on Ballistic 


Missiles, Los Angeles, Calif., June, 1957 
The author would like to acknowledge the assistance 


numerical computations and in drawing the Figure 








cataly sis 


Fic. 1 Schematic model for surface 
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Fic. 2. Nondimensional velocity, temperature, and concentra- 


tion profiles (altitude = 120,000 ft., Mach Number = 20) 
Vv molecular weight of the ith species 
\ Avogadro's Number 
hb static pressure 
Py Prandtl Number 
0 energy flux 
R universal gas constant 
1 temperature 
u x component of velocity, parallel to surface 
vA net mass rate of production of atomic species per unit volume by 
chemical reaction 
Wu net mass collision flux of atoms per unit area 
ZA (c4/C Ae) referred atom mass fraction 
‘ catalytic efficiency 
” similarity variable 
0 = 7/T, dimensionless temperature 
p density 
“ viscosity 


Subscripts 


1 atom 
é edge of boundary layer 
eq equilibrium 
VV molecule 
wall 
7 differentiation with respect to n 


\NALYSIS 


When dissociated air is treated as a binary mixture of atoms and 
molecules, the stagnation point boundary-layer equations be 


come,’ for steady laminar flow, 


{(1/Pr) [Le cacZay + Le" (0,/0)|\y 4 


i ¢ ZA t ws 2p(du, dx 0) (2 
(Ufa: _ | ae + (1/2 Pe/p) — Fg” = (3 

[(é,/ Pr Cale FT EpfO, + (1/Pr 9, Cpa — Com Le « eZ An tT 
Le" (@)/0)| — [tra(ha — hy)/2pT.(du./dx 0 (4 


which are, respectively, the conservation of atoms, momentum, 
and energy. These equations were programed for digital com 
putation on an IBM 704 computer by Paul Gordon, and a con- 
siderable number of exact solutions have been obtained utilizing 
variable Prandtl and Lewis Numbers and including the effects 
of thermal diffusion. We will restrict our discussion here to atom 


recombination on catalytic surfaces when gas pha reacti 
are frozen 

It has been previously assumed! ? that a fully catalytic 
corresponds to a zero atom concentration at the wall @., Cs 
Q. This surface boundary condition yields the maximum at 
diffusion flux to the surface and simultaneously the largest coy 
terdiffusion flux of molecules. However, this is not entirely sat; 


os 


factory because a microscopic condition is being violat 
face recombination of atoms which occurs ultimately by my 
of collision processes cannot occur if no atoms ire prese 
A more realistic interpretation requires that the surfac« 
concentration be small (to ensure maximum diffusion flux 
finite to satisfy microscopic conditions Indeed, it is 8 
} 


on equating the diffusion flux, obtained from the simultane 


solution of Eqs. (2), (3), and (4) above, 


Wy = V 2ewtue(due/dx) [(Le/Pr) caeZan + 


to the collision flux obtained from kinetic theory 
W, ya Ma/V 20 MART) [pa — (Pacg./Pa 


and taking y4 = 1.0 as the definition of a fully catalvti 
compatibility between the boundary-layer solution and kinet 
theory can usually be achieved when c4,, is on the order of 1 
Thus for computation purposes, one may still take 
being synonymous with a fully catalytic wall, if one holds 
the mental reservation that the precise surface concentratior 
atoms can only be obtained a posteriori by satisfying Eqs 
and (6 

With regard to noneatalytic walls, it has been assumed? t 
the atom concentration is a constant across the boundary laver 
i.@., C4,, = Cae iS Synonymous with Z4, 0 \ctually, if or 
includes the effects of thermal diffusion, the net diffusion flu 
of atoms vanishes (74 = 0) for a surface that can maintain aw 
concentration equal to 0.8 ¢4,.. The physical explanation for t 
is that the thermal diffusion coefficient is negative and atoms tet 
to diffuse thermally from a cold region to a hot region. 1 
means that for a cool wall, thermal diffusion opposes concentr 
tion diffusion and for a sufficiently large value of c4 CAe 
CAy = O.8 ca4-) the two fluxes balance each other completel 


For larger values of c4,, the net atom flux is away from the sur 


face, which implies a net surface dissociation of molecules 
any case, when ¢ { = (Ae, the atom concentration is not const 
See curve 4 of Fig. 2 

Having disposed of the two limiting cases, it is clear that sur 
faces having finite catalytic efliciency must maintain a surface 
atom concentration which lies between the limits 0 
O.SCa- 

\ typical calculation for the full range of catalytic efficiency ! 
shown is Fig. 3. The variation of surface energy transfer 


the net collision flux of atoms corresponds to the boundary-layer 











solutions of Fig. 2 
tz1.0 
a 
200-2. 
a ® re = *5.46 BTU ALTITUDE 120,000 FT 
| = FT.< SEC.5 WALL TEMPERATURE 800° R 
2 160} ~ : , 
°o B,* EFFECTIVE COLLISIONS 
alo TOTAL COLLISIONS 
u | 120: > 
ol. B=10.1X 10 = 
| 
aj 80 Q 393 | 
due 
— 2155 
-* 
-4 due 
2 t,=2.07X10 t:0 —_ 
Q rae 
= or-———- 2.76 -=200—* S 
= due J GUe , 
-40 dx dx J 
2 4 os 8 1.0 


Za, REFERRED ATOM MASS FRACTION AT SURFACE 


Fic. 3. Atom and energy flux vs. catalytic efficiency. 
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CONCLUSIONS 


With t xception of pyrex, (ya = 1.2 10~4), most common 
teri ve a catalytic efficiency for oxygen atom recombina 
t ceeds 50 X& 10 which means that they are vir 
fect catalysts at least at low wall temperature (see Fig 
Further experimental effort seems warranted if one wishes 
init the heat release due to catalytic recombination of 
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A Comment on the Methods of Calculating 
Natural Modes and Frequencies of Vibration 


Mechanical Engineer, Flutter and Vibration Group 
3s Aircraft Company, Inc., El Segundo, Calif 
94,1957 


jee PROCEDURE outlined by N. O. Myklestad,! and that sug 
rested by W. T. Thomson,? for the determination of natural 


des of vibration, contain expressions for certain elastic coefh 
ts, which might be modified to generalize and simplify the 


leulations It is noted that in the evaluation of these coeffi 
ind dp, the solutions become indeterminate for 


nt / 
S 1, dv, 


ease of beams of uniform sections. Each of the expressions 


juantities contain a term of the form 


CSC ¢ 
i 
1 /b? fla_h 
Shad 
limit (1/b")- | f(a, b = 0/0 (indeterminate 
fr » 
In the above expression, a represents the bending stiffness (£/) of 
station » and 6 represents the difference between the stiffness at 


in and that at station (7 + 1 
In the derivation of these expressions the graph of bending stiff 
less (EJ) of the beam plotted against the beam span Was approx! 
ited by a series of straight lines. The purpose of this note is 
) suggest that this approximation be replaced by one in which 
bending flexibility (1/E/ 
, is approximated by the series of straight lines 


graph, again plotted versus the 


t is seen from Fig. 1 that, for the section of the beam between 


ere A represents the flexibility (1/E/) at station , C represents 
flexibility (1/EJ) at station (n + 1), and the quantities 
vw, and d; represent the elastic coefficients of the My kle 
lmethod. <A similar procedure can be used to determine the 


velicients of Thomson’s method 
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n BEAM SPAN 


It might be noted that there exists an error in the relationship 
between the two sets of coefficients as given by Thomson Phe 


relationships should read 


ami = ami, a = a = dy 
du; = au aM ds dm d 


It is seen that these new expressions for the elastic coeflicients 
ire continuous for all finite values of A and ( Further, when 
1 equals C (uniform section of beam) the expressions yield the 
exact solutions for uniform beams. In addition to this, the above 
expressions greatly simplify the calculations by eliminating the 
logarithmic terms and certain products of terms which appear 
in the Mvklestad expression 

his approximation is in every respect as valid as that pre 
sented by Myklestad, since the error introduced by one approxi 
mation is within a small percentage of the error due to the other 
Phe method which will produce the better accuracy 1s dependent 
mtinuity and 


on the shape of the original curve Further, the « 


simplicity of the new expressions should aid greatly in the pro 
graming of the general problem of beam vibrations for auto 


matic computers 
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The Effect of External Sound on the Vortex 
Shedding From Cylinders 


Donald |. Garber* 
Lewis Flight Propulsion Laboratory, NACA, Cleveland, Oh 
December 30, 1958 


i pn INVESTIGATION reported herein was a part ol the develop- 
ment of an instrument utilizing the vortex-shedding tech 
nique! for stream-velocity measurement Phe possibility of 
resonant modes of acoustic oscillation excited by the vortex 
shedding suggested a study of the effect of sound waves on the 
vortex-shedding frequency 

The apparatus is shown in Fig. 1. Sound was introduced 


laterally into the vicinity of the vortices ro eliminate any 


standing wave reinforcement, the apparatus had only two pos 
5 


Since the wave 


sible resonant dimensions: 3/8 in. and 50 ft 


9 


length of the imposed sound was 2.5 in. to 2 ft., no reinforced 


resonance existed in the range studied 


* Undergraduate Student, Carnegie Institute of Technology 
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that the test apparatus gives Strouhal Numbers 
which are higher than those obtained by Roshko,! 


Fig. 2. Roshko, however, used a large duct, so that the cylinder 


It is notable 
as shown in 
aspect ratio high. The low aspect ratio data show 
Strouhal Number decreasing with Reynolds Number,? 
high Strouhal Number increasing with 
Reynolds Number that 


mediate aspect ratio exists where the Strouhal Number might be 


was very 
while the 
have 


aspect ratio data 


This suggests the possibility an inter- 
constant 

Fig. 3 illustrates the variation of vortex-shedding frequency 
with the introduction of sound frequencies for a constant Reyn- 
olds Number. At the 
tensity of sound at the point of entry into the duct is given as a 
Reflection from the opposite wall could 
When the frequency of 


bottom of Fig. 8, the approximate in- 
function of frequency 

possibly add 6 db to the 
the introduced sound is sufficiently lower than the vortex-shedding 
is no effect on 


values shown. 


frequency without sound input point (A), there 


the shedding frequency. However, as the frequency of the 


introduced sound is raised, the vortex frequency soon becomes 


“unstabilized’’ and begins to oscillate (B Gradually as the 


is increased, the vortex-shedding frequency is 


“locked-in” 
Over a limited range the vortex-shedding 


input frequency 


decreased until it is and the two frequencies are iden 


tical and constant (C 
frequency follows the frequency of the introduced signal until a 
frequency is reached that is higher than the original vortex fre- 
an ex 


quency (D). Raising of the input frequency results in 


treme ‘‘unstabilized”’ region where the vortex frequency oscillates 


over a considerable range (FE) \t still higher input frequencies, 








OSCILLOSCOPE 
| | 
ISIGNAL | | lFRE- 
|GENER-|_| | |QUENCY 
ill | |METER | 
= a 
—4 —> S«_ 
AMP-— ee | |HOT-WIRE 
[Liter ~~ — 


— 


HOT WIRE PARALLEL 
TO BLUFF BODY 


° 7 50' TUBING— ———= 
60.4" ys 8 TO VALVE & VACUUM=> 











“BLUFF BODY 
0.1"D STAINLESS- 
STEEL TUBE 


Fic. 1 





28} 


Ss 


PRESENT INVESTIGATION 


—— ROSHKO 


NUMBER, 
tw 
b 





STROUHAL 
n 
° 


0S NUMBER Re 


REYNOL 


Fic. 2 


AERONAU 


T 
































ICAL SCIENCES LPRIL, 1958 
2 L3§ _f —— 
¢ 488 ° 
25a ee 
x nae 
\~ue 
bay” : iz 
ae a 
VA _, Pitta. 
Fic. 3 
the vortex-shedding frequency can be decreased to a value t 


is lower than the vortex-shedding frequency with no sound intr 
duction (F 

As the Reynolds Number ts increased, the vortex signal becomes 
and the vortex-shedding frequency is less subject t 
interference. Whether to the of the vor 
tex signal at the higher Reynolds Numbers or to the increase i: 


less regular, 
this is due irregularity 
the kinetic energy of the system is not known. 

The results of the investigation indicate that ducting resonance 
sources affect the frequency of vorte 


or external sound may 


shedding from cylinders 
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Remarks on Systematizing the Calculations of 
the Natural Modes of a Free Structure 


R. Kappus and D. Clerc 

Office National d’Etudes et de Recherches Aéronautique 
Chatillon-sous-Bagneux (Seine), France 

December 31, 1957 


DISCUSSION 


[ IS WELL KNOWN that the calculation of the natural modes 
of a structure with a finite but rather high number of vibrat 
usually requires iterative methods. In addition, wher 


body equilibrium must bi 


ing masses 
the vibrating body is free in space, fre¢ 
taken into account 

This free-body equilibrium can be dealt with once 
with the aid of 


and for 


in an absolutely sure and systematic fashion 


geometric matrix 7, as the authors pointed out in detail in 1955 


This general scheme, which leads in a very simple manner to th 


classical matrix relation for a characteristic-value problem, tr 
quires only the a priori establishment of certain matrices 
geometric, mass, and elastic information. Thus the engineer 


need not enter explicitly into the numerical details of the free 


body equilibrium relations. The method also avoids reit! 


ducing these relations at each iteration step as is sometimes done 

e, for example, references 2 and 3 
In reference 1 we treated a more general case (the vibrating 
structure being additionally tied to a jig by means of certall 
), and details for handling the flexibility matrices relative 
see chapter III). In 


outlined 


springs 
to the various structural parts were given 


this reference two different calculation routes were 
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Of the two, that which consists in undertaking the mode and 
wracteristic-value calculation in terms of the relative displace- 
ents (columns ? proved to be much superior to that using 

jirectly the absolute mode displacements (columns p"')), which 
re the ones ultimately sought.f We have therefore considered 


+ worthwhile here to extract from reference 1, for the most com 
nonly encountered case of a completely free structure, the essen 


lations entering into the calculation of the relative modes 





through presentation of an outline of the calculation scheme 


together with a simplified example. In addition, we supply cer 


un new det iils 
It might further be remarked that a somewhat similar procedure 
Mack 


based on tensor notation has been developed by C. E 


OUTLINE OF THE CALCULATION SCHEMI 


1) Define the column p* consisting of the m* (<6) absolute 


linear and angular) displacements (with respect to its own rest 


sition) of a conveniently chosen reference system adhering to 
} ; 


t define the column p consisting of the m absolute 


e structure 
placements of the vibrating masses of the rest of the structure 


In the simplified example of Fig. 1 we choose the reference system 


is the trihedron defined by the xyz-axis system with origin at the 


rbitrary point 0, and we endow it with certain mass properties 


see below \ssuming, as usual, that the only important dis 


placements are those perpendicular to the xy-plane, we shall 
parallel 


take into account (a) as p*-displacements, translation 


to the z-axis and the two rotations 62) and 6,9; and (b) as p-dis 
placements, linear displacements w parallel to the z-axis of all 


he other points (i.e., points 1 to 6) and, in addition, the rotations 





points’’ 7 and 8 (which do not necessarily coincide with 


points 3 and 4 endowed with point masses J/; and AV; 





Thus, in the example chosen here, if it be temporarily assumed, 
ler to cover the most general case, that the structure has 


vinmetry, 














- = re. 
b , 
p. 
tT »,* 
| 7 
| ‘ ‘ 
f | p 0 and p= Pe = Ue 
Lf 4 p 4 
at ta 


Naturally, in the case of a real airplane the number m of elements 
in the column p will be much greater than 8 


2) Set up the geometric matrix T (consisting of m rows and 
m* columns) so that the column 7° p* represents the rigid-body 


displacements relative to the rest position of the reference system 


hus, in the example chosen (where m = 8, m* = 3), 

F 4) “4 

ry 9] 
l nN 
l x 

1 l \ \ where vy; = vs = 0 

l x 
] \ v. 
0 l 0 

L 0 ] 0) J 











3) Set up the symmetric square matrices ./*, 1/, and D, de- 


fined for all cases by the relations 


where f* and f (referring respectively to the reference system and 
the rest of the structure) are the columns of inertia forces, and 7 
is the column of the mass displacements relative to the reference 


some other notations have also been changed 


In reference 1 is called gq; 


in the present note 
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qa~x 
YZ) 
rf / 
Jo &fr 
y f 
x 
Oxo 4 





VW 
Vz 


Fic. 1 Perspective figure representing a stylized version of 
in airplane having 7 mass points (lo, 1/1, , We) and 3 mo 
ments of inertia in roll (9,9, 0,7, Ors 


“Ay 


matrix D (matrix of influence co 
efficients) is obtained by means of an elastic analysis which may 


system chosen. The flex 


often be rather complicated It should be noted that D (like 7 
is specified always with respect to the reference system (chosen 
more or less arbitrarily), which is therefore assumed to furnish the 
support” reactions. In what follows D is supposed known, its 


ealeulation not being an object of the present note (see, however, 





reference 1, chapter III The matrix D always has order m, 
and is symmetric and regular rhe same h for the inertt 
matrix MJ (which is even very often diagonal \s for the refer 
ence system inertia matrix \J* it has only order m* (<6) and can 


even be singular as is the case in our example where we assume 
that the reference system possesses no pitching moment of in 


ertia: QO, = 0. Thus, in this example, 


A) (0 t Z Fi 6 i 


VJ 


ul 


Note further that point 0 of the reference system may be any 
point whatever 


4) Carry out in successive steps the following matrix oper 


ations:{ 


MT, T(MT), M* + (TMT), M* + (M7 


H W* + TMI IM j 
S = (MT)H 5 
VV = \/-S i 
The matrix S is evidently symmetric and, moreover, always 


singular. The modified inertia matrix M is therefore also sym 


metric; it is generally re gular 


5) Solve, for the lowest modes desired, tl 


the classical char 


acteristic-value problem ** 


DMr® = (1/w,2)-r, ee | i 


r being the column representing the ith relative natural mode 
6) For each of the modes thus obtained, calculate the column 


~ A bar over a letter indicates transposition of the corresponding matrix 


we recallthat AB BA 


** For iterations concerning modes higher than the first, see for example 


reference 5. The procedure of this reference, called the method of ‘‘aug 
mented matrices matrices allongées), has recently been programed by 


Clere for the IBM 704 calculator 
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— THI» 8 
ind finally the column 

1 + 9 

1 constituting the 7th a ite natural m 
In the mple considered in Fig. 1 the matrices D and J\/ art 
of order s rhere exist, however, only 8 — 1 7 modes 
iracteristic frequencies; for it is clear that I* is here 
degenerate, and it can be shown that this is also true for 

hence for D.\J (see Addendum I of reference 6 
Of course, in practical cases the structural symmetry is always 
taken into account; as a consequence, the symmetric and anti 


symmetric vibrations are uncoupled In this case the following 


columns and matrices will be separately defined (see reference 6 
s, Ts, Ms*, Ms, Ds 
ind D4 7 2 VW VW D; 
It should also be mentioned that it is always possible to choose 


ypother reference system constituted, for example, by three points 


located in a plane reference 6 


See 

Note finally that the generalized masses can be calculated 
in two different wavs (see Addendum II of reference 6 

u Mp + p*OM*p* POM, 1, 2, 3, 
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A Note on Isentropic Compressible Flow of Air 
With Variable Specific Heats 





K. E. Tempelmeyer and L. Self 
AR 8, Inc Arnold Eng neering Deve opment Center 
ime anoma fenr 
Jan ary z 1958 
SYMBOLS 
b, ¢ coefficients 
specific heat at constant pressure, B.t.u./lb.m-°R 
J mechanical equivalent of heat, 778, ft.-lb. f Btu 
7 free-stream Mach Number 
absolute pressure, lb.f/ft 
R gas constant for air, 53.55, ft.- Ib.m-°R 
ibsolute temperature R 
| velocity, ft. sec 
ratio of specific heat 
é density, lb.m_ ft 
\ ratio of value for variable specific heat to value for constant 
specific heat 
Subscripts 
constant 
pressure ratio 
1 temperature ratio 
t Stagnation condition 
f density ratio 
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DISCUSSION 
apnea SSIBLE FLOW RELATIONSHIPS are usuall 
gas medium which has constant specific 
equation of state given by the perfect gas law 
satisfies these conditions accurately over onl 


pressure and temperature, u 


relationships for air flows at h Mach 





significant errors 


ire valid over considerably wider pressure 


se of the usual tab 


1 


11] 


u 


Number 


However, compressible fl 


w eC 


may be derived for air from the perfect gas equat 
the variation of specific heat with temperature is « 
Fig. 1 Equations derived with variable specific 
relatively accurate values when 

0 <p <3 atmos 300°R 1 1,50 
and in this range the additional complication arisi1 


of a more complex equation of state 
man*) may be avoided 

Compressible flow relationships have been 
1 for 


constant specific 


perfect gas with ecific 


variable sy} 


heats; however, because 


relationship employed, the variable specific 


somewhat difficult to use 


1¢ 


such as 


ile’ 


jua 


ind temperature 


10n OF 


he < 
MR 

g fron 
Catt 


aerived int 
its as W 
the con 
heat equat 


Simpler relationships for con 





isentropic flow of air have been derived by using the px 
equation of state and relatively simple but accurate expr 
for the c,-7 relationship 
The c,-7 curve has been separated into three 1 ge 
shown in Fig. 1, and empirical equations accurately fitted t 
section of the curve with a mean c, deviation of 0.2 per cent 
the values of reference 4 rhe resulting equations 
form 
Cp = + bf + cf 
Phe coefficients ¢ ind ¢ for each range are listed in Table 
Euler’s equation for isentropic flow may be writte s 
J Vdi fc, dT = 
If Eq. (1) is utilized in the integration of Eq. (2 
V2/2) + al + (b/2)7 c/3)1 const 
raBLeE | 
Temper 
ature 
Range, Coefficients 
Range R 
] 100-400 0.23893 0 
> 100) 
1.700 O.2318 0.1040 Kk 10 0.7166 & 
3 1,700 
1.500 0.2214 0.3: x 10 0.3776 > 
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_ gnation condition is assumed at one point im the use of equations based on constant is q signit t higl 
I may be simplified to express the following relation- temperatures. Consequently, if the pressurt temperature 
4 static d stagnation temperature stream ire wit the previously specifi i t f equa 
I rf isentropic flow of air with variable specific tions based on variable specific heats is dictats the magnitude 
= f the stag ti temperature s illust by 
al 2)J i } 7 R i+ 
1 } ) RF 7 
S giv ) Eq . be deter ne It ) = 
. R/J ) rhus, in ge l when the stagnation te greater tha 
) 
; WOSR., the use of equ s base \ € spec é s 
ve written when 7 and 7; lie wit single tem- . 
, i ter 
ge I y coml a iges D g the coell : 
. ‘ It s Is p sibie t ¢ ¢ thie 1 n I ae 
S 1 W 11k correspond to the range I their par : , 
1 vel how equati S for ust cross I e 4 K 
ture (see Table 1 Wu 
1 diving t / relations 
S gnation pressure ad densit ratios »« D , 
w 4 1 1 er gases W 1« ( id be use tor g é e fi 
48 valuating the following for f the gener energy) ee 
level 1s wit t presc! s 
f sentropic flow : 
> . 
t = dl ( REFERENCES 
S (dp S:Al ; 
f } t is combined with the perfect gus law a Eq i 
) 
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: 1 " x ~ 
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perature ratios for an isentropic process for air with variable York, 194 
: ts, whic 1 \ v written for flow process as 
r f R 2 + 
: 1ges bh +t / / Q 
a substitution into Eq. (6) is accomplished in slightly 
re 
‘ t manner, the expression for density ratio is obtained as Propagation of Weak Waves in a Dissociated 
* 
k 1 R(bCI £(c/2 ™ Gas 
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Ss wit nd then equating the two equati For illustrative pur N A RECENT PAPER, Resler has provided tu fe the « 
ses, the pressure ratio equation, when 7’ and 7; are in ranges I acteristic, or sound, velocity of a gas out of ther lwnamix 
s respectively, is given below All coefficients and cor quilibrium. This note is concerned wit related topic 
essil stant ve been inserted and simplified | ens of k or acoustk 
Ct g f T00/7 T/400 0.003535 uthors ‘ yzed 
SS 7 ) j 0 0.7932 0 ect ssociati 
Its t with tl 
S | ee tions produce results which are in excelle igree cornti white enc 
€ esults obtained from other equations based on vari | - g 
t fron le specific heats, as in reference 1. The equations here, however, It is interesting to consider the equati vhich, f ssociate 
ft re easier to apply than those of reference 1. Additional com gas, replaces the usual acoustic wave equation \ssumption of 
risons with results produced by equations based on the Beattie particular wave fort e.g ' ic . I 
Bridgeman equation of state with variable specific heats* are purposes of illustrati .ttenti — , fect of 
so favorabl For example, at JJ = 3.5 and 7; = 1,100°R dissociation of a diatomic g so that ‘ g 
Beattie-Bridgeman based equations® gave / 326°R., iverage gas constant may be written asf 
reas the equations derived here yield 7 = 325.5°R., and 
: sed on constant specific heat! indicate 7 = 322.6°R . 6 — a@)/2) Rol 
The errors associated with the use of the usual compressible : ' 
i T r - 
nay be indicated by considering a typical flight 
f JJ = 5.0 at an altitude of 140,000 ft ie altitude For small departures from equilibrium, one rite that the 
the static conditions as p = 4.3 Ib.f/ft.2 a Tl = 485°R., dissociation rate is proportional to the difference between the 
upressible flow relationships with constant specific actual temperature and the temperature at which the existing 
ts are used, the stagnation conditions may be computed as dissociation level would be in equilibriut 
t 2,275 Ib.f /ft and T:)cp = 2,909°R Da /Ot et — ; 
stagnation conditions are computed with variable specific where K is a positive reaction rate dt _ one ee 
ts, Eq. (4) may be used to determine the temperature; and, Pat 
: of the gas 
: range | 3, Eq. (10) may Using Eqs. (1)-(3) in conjunction with the hi rized disturb 
( n variable specific 
ance equations 
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p(Ou/et) + (Op/de) = writing 
Op/Ot) + p(Qu/de T. = T + (O7./0a)pAa + (O7./Op)a Ap (8) 
and (c) insert Eq. (8) into the equation obtained in and its 


p(oE/ot) + p(ou/de) = VU : . S Res : 
time derivative, combining the results in such a way as to give 


pmens ' 0/dt) [(0?7'/dt?) — a,? (0?T/dE*)] 4 

(where bars represent mean properties in the absence of disturb- k[(0?7'/ot?) — a,? (0?T/0¢ 
ance), one may derive the analog of the acoustic equation. The where 2 Q ( : =i 
steps are (a) manipulating Eqs. (1), (2), and (4)-(7), obtain an 


9 


. . . . 6 . . ~ whic f > > the “‘frozen’’ speed of 1 radi 
equation involving only 7 and a; (b) further linearize Eq. (3) by vhich may be shown to be the ‘‘frozen’’ speed of sound (for a dig. 


turbance during which a would remain constant); where, defining 


QO = 2q0/Ro7 


(9 + &) (O0T./0a)p + pi Q — [8/(1 + &)] }(O7./Op)o + TIO+1 
(-@ oT, T (11) 


which may be shown to be the “equilibrium’’ speed of sound (for 
a disturbance during which the dissociation level follows its ’ ; ; 
sheath : : : : ; ihe sound speed is generally the higher, there would be a tendency for 

equilibrium relationship with the properties of state prevailing : ‘ . 
‘ ‘ ; oi : e 3 the high and low frequency components to be sorted out, with 

at each instant, i.e., 7 = 7.), and where & is a positive quantity : : : : 
> the high frequency waves running ahead at the frozen speed. In 
-\ ’ | [T ey =\] ' ‘ : ‘ 

k K;) (OT./Oa)p + [T/(7 — &] (Q— 1); 12 fact, recent shock-tube experiments® conducted at Cornell Aero- 
By reference to Eq. (3), it is clear that & is essentially a reciprocal nautical Laboratory in dissociated bromine and oxygen indicate 
relaxation time, or a dissociation rate. that the head of a steep expansion wave does indeed travel at qa 
Eq. (9), it may be noted, has the same form as the equation sound speed frozen in respect both to molecular vibration and 


derived by Morrison? for the propagation of elastic waves when a dissociation 


time lag exists between stress and strain. While the foregoing conclusions are substantially known,® 4 jt 


From Eq. (9), which, for a dissociated gas, replaces the acoustic is perhaps noteworthy that they follow directly from the differ- 


ential Eq. (9 Since relaxation times are long at low densities, 
the flow behind a shock during hypersonic flight at extreme alti- 
tudes would involve wave processes complicated by the consider- 
ations discussed herein. The present approach, further devel- 


oped, might be helpful in formulating an aerodynamic theory of 


wave equation for temperature, certain conclusions may immedi 
ately be drawn regarding propagation of weak disturbances 

In general, both ‘“frozen”’ and ‘‘equilibrium’’ processes are in 
volved in a disturbance. In fact, the form of Eq. (9) implies 
that the frozen and equilibrium operators are out of phase during 
a periodic disturbance, so that the propagation speed must vary 
markedly during each cycle. Thus, a mechanism for wave-form 


weak disturbances for such a situation 


distortion is apparent even in the linear approximation, in con REFERENCES 
trast to the usual distortion or steepening of an acoustic wave, Resler, E. L., Jr., Characteri Sound Speed in Nonisentropicaae 
which is a nonlinear effect Flows With Nonequilibrium Thermodynamic State Journal of the Aero- 
For an oscillatory process, the presence of the extra time nautical Sciences, Vol. 24, No. 11, pp. 785-790, November, 1957 
derivative in the left-hand term may be replaced, in effect, by a Chu, Boa Teh, Wave Propagation and the Method of Characteristics im 
See MORE: : Se to Reacting Gas Mixtures with Applications to Hypersonic Flow,WADC TN- 

frequency If this frequeney is large compared to the dissoci 57.913. M 1957 
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ation rate (k), then evidently the frozen wave equation dominates ; 

raat steales seep hides juation dominate 3 Lighthill, M. J., Survey in Mecha» Viscosity Effects in Sound Waves 
the motion, and if the frequency is small compared to k, the propa of Finite Amplitude,’’ (G Taylor 70th Anniversary Volume), Cambridge 


gation is governed by the equilibrium wave equation. In fact, University Press, 1956 


harmonic solutions of Eq. (9) may easily be obtained which show ‘ Morrison, J. A., Wave Propagation in R Voigt Material and Visco 
1 . 

this change of sound speed with frequency, or ‘velocity disper Elastic Materials with Three-Parametes dels, Quart. Appl. Math., Vol. 

14, pp. 153-169, 1956 


sion,” in the terminology of the literature of ultrasonic absorption, 
is well as the attendant energy absorption in the frequency rang * Treanor, C. E., Golian, T. C., Hertzberg, A., Measurement of Pig 
as ‘Ilas Ps ond: ne abs ( > > > ‘ e : 

8: I | : - Properties in a Dissociated Gas, Paper presented at the Tenth Meeting of 
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where the change in sound speed occurs. 
Equations for Reactive 


In general, any given wave packet would include both low . Wook. WW. wad Rikwo0d. 1... Charact 
Phys., Vol. 27, No. 2, p. 596, 1957 


and high frequency constituents. Therefore, since the frozen Flow, Jour. Chem 
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This Symposium will be mainly concerned with the present position in regard to research being carried out on Rarefied 
It will include theoretical and experimental studies on (1) aerodynamic and 
(2) fundamental physical values that occurred in the preceding 
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